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HYDROLOGY OF LAKE ONTARIO 


lay Morton, 1 and H. B, Rosenberg2 


SYNOPSIS 


description of the physical and climatological features of the Lake 

rio basin is presented. The effects of climate and topography on the 

ly of water into Lake Ontario are discussed. Short and long term records 
ws and supplies of Lake Ontario are related to the factors causing varia- 
in the water supplies. Detailed descriptions of the source of supply are 

- in relation to seasonal and areal distribution of precipitation, tempera- 

and runoff. Flow records of many of the large tributaries of Lake 

rio are analysed and the possibility of long and short-term forecasting is 

issed, 


INTRODUCTION 


1e International St. Lawrence River Board of Control is currently de- 
ing plans for regulating the levels and outflows of Lake Ontario in con- 
ion with the St. Lawrence Seaway and Power projects. The technical 

ss are being performed jointly by various governmental agencies in the 
d States and Canada. The Water Resources Branch of the Department of 
ern Affairs and National Resources is one of the Canadian agencies as- 
ted with this work. The problem is international in all its ramifications 
the International Boundary traverses not only Lake Ontario and the St. 
ence River to the vicinity of Cornwall, Ontario, but also the upper Great 
3 and their connecting channels with the exception of Lake Michigan. 


“Discussion open until October 1, 1959. To extend the closing date one month, a 
itten request must be filed with the Executive Secretary, ASCE. Paper 2017 is 
rt of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
nerican Society of Civil Engineers, Vol. 85, No. HY 5, May, 1959. 


dr, Engr., Water Resources Branch, Dept. of Northern Affairs and 
tional Resources, Ottawa, Ontario, Canada. 
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In the preliminary regulation studies which have already been made, pre- 
conceived rules have been developed and tested under historical water suppl 
conditions. The water supplies utilized in testing the plans have been deriv 
from historical outflow and water level changes and results beneficial to all 
interests have been achieved without the necessity for a detailed study of th 
sources of supply and the other hydrologic factors. If further refinements 
are to be made, operating experience as well as a more detailed knowledge ¢ 
hydrologic data will be required. 

The question may be asked as to what may be gained from a better under- 
standing of the various factors contributing to the water supply to Lake 
Ontario. Firstly, the knowledge gained from such a study will contribute to 
the confidence of those responsible for the operation of the works when the 
actual water supplies to be regulated will assuredly be different from those 
which have occurred in the past which have been the basis for the develop- 
ment of the regulation plan. The hydrology of Lake Ontario is very complex 
so that a thorough understanding of it will be necessary for efficient, flexible 
operation of the control works. Secondly, such studies may enable the de- — 
velopment of long-term and short-term forecasts of water supplies. This 
would contribute to the improved use of the available storage on Lake Ontari 
Thirdly, this knowledge can make possible an assessment of the possible ad- 
vantages to Lake Ontario of the development of regulation facilities on the 
upper Great Lakes. b 

This study was prepared with these considerations in mind. Its principal 
purpose was to investigate in a general way all the factors influencing the in. 
flows to Lake Ontario and to set them forth in a concise form. The effects ‘ 
precipitation, temperature, lake areas and the physiography of the drainage 
basin on the character of the inflow are analysed. From the better under- _ 
standing of inflow characteristics gained from this analysis an assessment is 
made of the probable effects on Lake Ontario of the regulation of the upper — 
Great Lakes and the possibility of long and short-term forecasts of water _ 
supplies, 


Basin Characteristics 
General 


The Great Lakes constitute one of the finest natural regulatory systems i 
the world, At the outlet of Lake Ontario the total drainage area of the sys- 
tem is about 300,000 square miles of which approximately one-third is wat 
surface. The excess water from the upper Great Lakes, i.e. Superior, 
Michigan-Huron, St. Clair and Erie, passes through Lake Ontario and thence 
to the St. Lawrence River and the sea. This water which is transmitted to ~ 
Lake Ontario from the upper Great Lakes is of great importance to an unde 
standing of the hydrology of Lake Ontario since it comprises approximately 
85% of the total water supply to the lake. Figure 1 is a map of the drainage 


basin of the Great Lakes, Figure 2 presents a summary of the pertinent 
physical data. 


Upper Great Lakes 


The characteristics of the water supply to the upper Great Lakes are s U 
ject to the divergencies in climate, topography and geology that can occur i 
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an area of 260,000 square miles extending 750 miles from east to west, an 
650 miles from north to south, 

The region contributing to the water supply of Lake Superior and the 
northern shore of Lake Huron is generally rough, wooded country typical o 
the Precambrian Shield. To the south and east of Lake Erie lies the Alleg! 
plateau which is part of the Appalachian formation. Between these two for 
mations lie the Great Lakes - St. Lawrence Lowlands which comprise the 
greater part of the area draining into the upper Great Lakes system. The 
glacial periods had a great effect on the entire region, grinding down and 
denuding the Precambrian Shield, flooding the low areas and depositing sil 
and sands, while leaving morainic hills, till plains, drumlins and eskers or 
the more elevated parts and generally disarranging old drainage patterns. 

The upper Great Lakes region lies in the path of the prevailing westerly 
winds and is particularly subject to winter storms. In the lee of the lakes 
areas of more moderate temperatures and extremely high snowfalls. Gen- 
erally, the temperature increases from north to south and the precipitatior 
increases from northwest to southeast. The northern part of the drainage 
basin has a mean annual temperature of 31°F. with an average length of the 
frost-free period of only 75 days. Near the western end of Lake Erie the 
mean annual temperature exceeds 48°F. and the average length of frost-fre 
period exceeds 175 days. The average annual rainfall varies from 25 inch 
west of Lake Superior, to over 40 inches to the south and east of Lake Erie 

Although the inflows to the upper Great Lakes and the evaporation losse 
vary widely from year to year, from season to season and from lake to lak 
the diversified runoff characteristics and the great expanses of lake surface 
contribute to a relatively uniform inflow to Lake Ontario through the Niaga 
River. The monthly mean outflows from Lake Erie have varied from a hig 
of 254,000 cubic feet per second to a low of 121,000 cubic feet per second, 

For comparative purposes the outflows from the various lakes which ar 
shown on Figure 2 and the outflows from Lake Erie just mentioned are thos 
that have been derived from actual records. Certain adjustments would be 
necessary to make these flows consonant with present conditions of diversi 
into and out of the Great Lakes System and any further references to flow 
and levels will include these adjustments. A net diversion into the Great 4 
Lakes system of 1900 cubic feet per second is presently in effect. This in: 
cludes the 3100 cubic feet per second diverted out of Lake Michigan at Chi 
and the 5000 cubic feet per second diverted into Lake Superior from the 
Albany River system. 


Lake Ontario Drainage Basin 


General 


The local drainage basin of Lake Ontario has an area of 34,630 square 
miles including 7,540 square miles of lake surface or over 20 per cent of 
local drainage basin. Of the drainage area, 18,710 square miles are in 
United States with the remaining 15,920 Square miles in Canada. From w 
to east the greatest length of the drainage basin is 280 miles and the gre 
width from north to south is 240 miles. The lake has an average depth of. 
proximately 250 feet with a maximum recorded depth of almost 800 feet. 
average water level has been about elevation 246 and the maximum range 
monthly stages has been about 6.6 feet. 
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siography 
‘he geological features which dominate the local drainage area of Lake 
rio are the Great Lakes - St. Lawrence Lowlands, the Precambrian 
Id, the Adirondack Mountains and the Allegheny plateau. These features 
highlighted on the contour map of the local drainage basin presented on 
ire 7-A. 
‘he Great Lakes - St. Lawrence Lowlands follow the shoreline of Lake 
io, striking inland for various distances. The lower lands close to the 
-became covered with silts and sand. At the edge of this lower area are 
e-cut cliffs, beaches, bars and deltas. To the south and west of the lake 
Niagara escarpment forms another boundary to this lower-lying area. In 
higher areas the relief was formed by the shoving and depositing action 
1e glaciers which left behind the morainic hills, drumlins, eskers and till 
ns. The glaciers also blocked off the established southward drainage of 
Finger Lakes in New York State causing them to flow northward into Lake 
ario. This area includes the headwaters of the small tributaries to the 
tern portion of Lake Ontario. 
The Allegheny plateau forms the northern edge of the Appalachian forma- 
. It extends across the southern end of the drainage basin in an east-west 
.etion from Syracuse to Buffalo and is deeply indented by the Finger Lakes. 
plateau slopes upward to the south and ranges in elevation from 1,000 to 
10 feet above sea level. The terrain is very irregular and hilly and has 
n largely cleared of its original forest growth. The highest point that is 
ieved in this formation in the drainage basin is 2,550 feet above sea level. 
headwaters of the Genessee River and the southern and western tribu- 
es of the Oswego River lie within this area. 
The northern part of the drainage basin is located within the Precambrian 
ion. Rock knobs or hills with shallow soil depths predominate. Little of 
area is cultivated and it is still largely covered with forest. In the ex- 
ne north of the drainage basin elevations of over 1,500 feet above sea 
1] are reached but the general elevation decreases to the south and east to 
Thousand Islands Section of the St. Lawrence River. This region includes 
headwaters of the Trent, Moira and Napanee Rivers. 
[he Adirondack plateau is an elevated extension of the Precambrian Shield 
th of the St. Lawrence River and east of Lake Ontario. The plateau is well 
ded with very rough topography. The highest hill in the drainage basin 
ches an elevation of approximately 3,700 feet although elevations of over 
0 feet are reached farther east. The area includes the headwaters of the 
yegatchie and Black Rivers and the northeastern tributaries of the Oswego 
er. The valley of the Black River indents the plateau in a striking fashion. 


mate 
The Lake Ontario drainage basin lies along the track of many of the low 
ssure areas which sweep across the northern part of North America from 
t to east. This results in stormy changeable weather with wide variations 
emperature but remarkably uniform annual precipitation. For the drain- 
basin as a whole the annual precipitation since 1870 has varied from 25 to 
ches with a mean of 33 inches. The prevailing westerly winds combined 
1 the topography play a major part in the areal distribution of rainfall. The 
rage annual precipitation varies from 27 inches in the area below the 

gara escarpment to over 50 inches in the Adirondacks to the east of the 

2, The snowfall is even more variable with the prevailing winds picking up 
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a large quantity of moisture as they pass over the open lake and depositing 
to the lee of the lake. The average snowfall varies from 40 inches at the wi 
of the lake at St. Catherines to 175 inches on the eastern end, south of Wate 
town. 

The temperatures in the drainage basin are variable with the area below 
the Niagara escarpment having the highest temperatures, and the northern 
part of the drainage basin and the higher elevations of the Adirondacks havi 
the lowest temperatures. The average annual temperature is about 45°F wi 
extremes of temperatures of 105°F and -40°F having been recorded. 


Water Supply to Lake Ontario 


General 


An analysis of the water supply conditions for Lake Ontario will be de- 
veloped along two different but complimentary approaches. An analysis wil 
first be made of the long term records and then supplemented by considerat 
of shorter term records with more comprehensive coverage. Because of 
their importance to navigation, continuous water level records for the Grea 
Lakes and the St. Lawrence River have been kept since 1860. It has been 
possible, therefore, to construct records of the water supplies and outflows 
of the lakes for a period of almost one hundred years. Records of tributary 
flows and detailed climatological records, however, have been available for 
a much shorter period of time. 


Net Total Water Supply to Lake Ontario 


The net total water supply to Lake Ontario is defined as the inflow to the 
lake which is available for outflow. It is equal to the inflow from the tribu- 
tary area plus the precipitation on the lake surface less the evaporation fro 
the lake surface. The records of monthly net total water supplies have bee 
based on the monthly mean discharge of the St. Lawrence River at Iroquois, 
Ontario, with the water levels at the beginning and end of the month measur 
at selected gauges around the lake shore. A change of one foot in the eleva 
tion of Lake Ontario in one month indicates that the difference between the 
net total water supply and the outflow has been 80,000 cubic feet per second 
The recorded net total water supplies have been adjusted to make them con 
sonant with a continuous diversion into the upper Great Lakes of 1900 cubic 
feet per second. . 

The net total water supply to Lake Ontario may be broken down into tw 
principal components, the outflows from Lake Erie and the net local water 
supply. The next total water supply has averaged 246, 000 cubic feet per 
second with 210,000 cubic feet per second contributed from Lake Erie and | 
the balance from the local drainage basin. The average monthly distribut 
of these two components, as well as the total is shown on Figure 3-A. In| 
general the net total water supply reaches a maximum in the month of Apr 
and a minimum in the month of October. Most of the variation in the avers 
monthly distribution is due to effect of the net local water supply. Figure 
illustrates the seasonal distribution of the extremes of the net total water 
supply. The highest and the lowest lines are the maximum and minimum 
each month regardless of the year in which they occurred. Generally, the 
extremes vary in the same manner as the averages. The monthly mean ne 
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vater supplies have varied from 416,000 cubic feet per second to 

6 cubic feet per second. The crosshatched area in Figure 3 represents 
se of one standard deviation on either side of the mean. The standard 
ion varies little throughout the year. 

nth-to-month variations in the net total water supplies are generally 
small in the late summer and early fall but can be very large during the 
~and spring. An extreme value occurred in 1936 when the variation in 
mnthly means from February to March was 178,000 cubic feet per 

i, An example of extreme variation for three months occurred in 1913 
he monthly mean net total water supplies for January, February and 

| were 324,000, 220,000 and 326,000 cubic feet per second. 

> variation in the mean annual net total water supplies is shown on a 
mcy distribution on Figure 4. The extreme variation for the yearly 

; has been 93,000 cubic feet per second or 38 per cent of the average. 


w from Lake Erie 


> greater part of the excess water from the upper Great Lakes passes 
h the Niagara River into Lake Ontario. To obtain the total outflows 
Lake Erie the quantities of water diverted for power and navigation have 
dded to the Niagara River flows. These total outflows have also been 
ed for a continuous net diversion into the upper Great Lakes of 1900 
feet per second. 
> flow in the Niagara River is essentially dependent on the water levels 
e Erie. There are three categories of variation in the outflow from 
rie, i.e. transient, seasonal and long term. 
> transient variations are caused by winds and seiches (oscillation of 
aters due to changes in barometric pressure) which may affect the 
level of Lake Erie at the head of Niagara River by more than ten feet, 
ice jams forming at the outlet of the lake at the head of the Niagara 
A large variation in flow resulting from a recent-storm is shown on 
rly hydrograph in Figure 5. 
sonal variations in the outflows from Lake Erie are small due to the 
ating effect of the storage on Lake Erie and the other upper Great 
_ Generally the outflow reaches its peak during the month of June and 
est value during the months of January or February. The monthly 
1 of the average outflows is shown in Figure 3-A. 
diversity of hydrologic factors affecting the water supplies to the up- 
eat Lakes and the retarding effect of the great volume of storage avail- 
esult in long-term variations in the outflows from Lake Erie. Short- 
ariations in water supply are absorbed by the storage with little effect 
outflows of Lake Erie but long-term variations gradually change the 
3 so that the seasonal and transient variations are superimposed upon 
r or lower base. 
*xample of the seasonal and long-term variation of the outflows from 
rie from 1926 to 1935 inclusive is shown in Figure 6. Also shown is 
graph of average outflows and a plot of the accumulated deviations 
verage of the precipitation on the drainage basin of the upper Great 
‘The flow gradually built up from much below average in 1926 to near 
ximum recorded in 1929 and then decreased to the minimum recorded 
. These long-term variations were obviously in response to long term 
ons in the precipitation on the drainage basin. The figure very clearly 
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shows the seasonal variations which were superimposed on a base that was 
much above or below average. Although the precipitation for 1932 and the 

early part of 1933 was average, the deficiency of water supply to the upper 
Great Lakes which occurred in the extremely dry period in 1931 had so de: 
pleted the storage on the lakes that the outflow from Lake Erie continued t 
be much below average. 

The seasonal distribution of extreme high and low monthly mean outflow 
from Lake Erie, regardless of when they occurred, is shown in Figure 3-C 
The trend of these extremes follows the trend of average seasonal variatio: 
quite closely. The monthly mean outflows have varied from a maximum of 
258,000 cubic feet per second to a minimum of 131,000 cubic feet per secor 
or a range of about 60% relative to the mean. The minimum occurred in 
February 1936 and a study of available records indicates that it was causet 
in part by a combination of easterly winds and an ice jam at the head of the 
Niagara River. 

A range of one standard deviation on either side of the average is also 
shown. The standard deviation varies very little during the year. 

The variations in the annual mean outflows from Lake Erie are shown 0 
Figure 4 in the form of a frequency distribution, The maximum variation ¢ 
the annual means has been 68,000 cubic feet per second or about 32% of the 
mean, 


Net Local Water Supply 


The net local water supply to Lake Ontario is the difference between the 
net total water supply and the outflow from Lake Erie. Since it is based on 
the difference between two large values it is possible that the percentage q 
error in its determination may be high. The net local water supply is com 
posed of the land runoff from the local drainage basin of Lake Ontario plus 
the precipitation on the lake surface less the evaporation from the lake sur 
face. The seasonal distribution of the average net local water supply is ~ 
shown on Figure 3-A. The maximum monthly average was 90,000 cubic fe 
per second in April and the minimum was 6,000 cubic feet per second in j 
September, 3 

The seasonal distribution of the maximum and minimum net local water 
supply for each month is illustrated on Figure 3-D. It may be noted that di 
ing the late summer, fall, and early winter, the minimum net local water _ 
supplies are negative quantities. These negative quantities result when he 
evaporation from the lake surface exceeds the sum of the precipitation on 
the lake surface and the land runoff. The net local water supply varies w 
with a maximum monthly mean of 188,000 cubic feet per second and a mini 
mum monthly mean of -32,000 cubic feet per second or a range of about 6] 
of the mean. The seasonal distribution of the extremes bears a close rela 
tion to the seasonal distribution of the averages. A range of one standard 
deviation on either side of the mean for the net local water supply is also 
shown. The standard deviation varies from a maximum in April to a mi ni 
mum in August. The standard deviations for the net local water supply avi 
age 70 per cent of the standard deviations for the net total water supply. — 

The net local water supply varies widely from month to month, notably 
the winter and spring seasons. The maximum variation occurred in 1936. 
the variation from February to March was 138,000 cubic feet per second, 
During the year 1913 the mean monthly net local water supplies varied fro 
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000 c.f.s. in January, to 5,000 c.f.s. in February, to 110,000 cubic feet 
second in March. 

igure 4 shows the variation in the annual mean net local water supplies 
frequency distribution. The maximum variation of the annual means has 
39,000 c.f.s. or 109% of the average net local water supply. This varia- 
is 62%of the maximum variation in the annual outflows from Lake Erie. 
slope of the frequency distribution for the net total water Supply is greater 
the slopes for either the outflows from Lake Erie or the net local water 
ly which indicates that there is a tendency for both components to be high 
ww together. However, the extreme values of both components have not 
cided. 

more thorough understanding of the variations in the net local water 

ly requires a more detailed study of the land runoff in the local drainage 
1 and the precipitation-evaporation balance on the lake surface. 


| Runoff 

and runoff is the inflow to Lake Ontario from the local drainage area ex- 

ve of the lake surface. This area is 27,090 square miles. The major 

taries are the Oswegatchie, Black, Oswego, Genessee, Trent and Moira 

rs which contribute more than half of the net local water supply to Lake 

rio and form the basis of the analysis of land runoff in this study. 

ne Oswegatchie and Black Rivers drain a large part of the Adirondack 

au within the Lake Ontario drainage basin. The flows of both tributaries 

ffected somewhat by the operation of storage reservoirs in the head- 

's. The Oswego River is the largest tributary of Lake Ontario exclusive 

> Niagara River. It drains the southeastern part of the drainage basin 

headwaters in both the Adirondack and Allegheny plateaus and serves as 

tlet for the Finger Lakes. Its flow is regulated appreciably in conjunc- 

vith the operation of the New York State Barge Canal. The Genessee 

"rises in the Allegheny plateau of Pennsylvania and drains the south- 

al part of the local drainage basin, Its flows are regulated slightly by 

eam storage reservoirs. The Trent River rises in the extreme northern 

of the drainage basin in the Precambrian Shield and drains the north- 

al portion of the local drainage basin. It has numerous lakes along its 

e and the flows are regulated for power and for operation of the Trent 

system. The Moira River has its headwaters in the Precambrian 

1 to the east of the Trent River and flows through a region with very 

soil cover. The flows are not regulated. The period of record for the 

stream gauging stations on these rivers varies from 26 to 45 years. 

ne of the major tributaries arise in the Great Lakes - St. Lawrence 

nds to the west of the lake or in the northeastern part of the drainage 
Two smaller rivers, the Credit and the Napanee, are considered to 

resentative of these areas. The Credit River rises above the Niagara 

pment and flows into Lake Ontario west of Toronto. The Napanee drain- 

1Sin is similar to that of the Moira. The flows of the Credit and Napanee 

S are substantially unregulated. The periods of record at the down- 

n gauging stations are 10 years for the Credit River and 21 years for 

ipanee River. 2 

e density of stream measurement stations on the local drainage basin is 

ate for the determination of the areal distribution of runoff. Figure 7-D 

lap of the drainage basin showing the distribution of mean annual runoff 

1es. The area of maximum runoff is east of Lake Ontario in the 


| 
a 
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headwaters of the Black and Oneida Rivers. Here the average runoff is as 
high as 35 inches per year. The area of lowest runoff follows the shoreline 
of the western half of the lake where the average runoff is as low as 10 incl 
Generally the runoff increases with distance from the lake. The reason for 
this distribution is apparent from inspection of Figure 7-B showing isohyet 
of mean annual precipitation. From a comparison of Figures 7-B and 7-D, 
may be deduced that the annual average water loss is generally about 20 ine 
although it is somewhat less to the east of the lake. It may be seen that the 
areas of greatest precipitation and runoff are coincident with the highest la 
In addition, the prevailing westerly winds, after crossing Lake Ontario, are 
affected orographically by the Adirondacks. This is especially true in the 
winter time where the winds sweeping over the open water pick up large 
quantities of moisture and drop it in the form of snow in the lee of the lake, 
The distribution of the mean annual snowfall is shown in Figure 7-C. 

The seasonal distribution of average precipitation on the drainage areai 
shown on Figures 8A, 8B, 8C and 8D. The areal distribution of precipitatio 
during any one season generally follows the same pattern as that for the ye 
A remarkable feature illustrated by these figures is that the precipitation — 
does not vary widely from season to season, Figures 9A, 9B, 9C and 9D sh 
the areal variation of average temperatures for the months, January, April 
July and October respectively. The decrease in temperature with both lati: 
tude and altitude is apparent. The range of areal variation of temperature 
greatest in the winter. The seasonal and areal variations in runoff are sho 
on Figures 10A, 10B, 10C and 10D. During the winter season the effectpf 
areal variations of temperature on the runoff is very marked. In the:northr 
sections of the basin and in the Adirondacks, where the temperatures are © 
lowest, the precipitation is stored in the form of snow so that the runoff is" 
much less than the precipitation. In the lower and more southern part of tt 
drainage basin a combination of higher temperatures, small evaporation an 
transpiration losses contribute to a runoff that is almost as great as the pi 
cipitation. During the spring months the situation is reversed. In the oa 
where the precipitation was stored in the form of snow the runoff is either 
greater or slightly less than, the precipitation, whereas in the other areas 
is considerably less. During the summers the large evaporation and trans 
piration losses result in a generally low runoff. In the fall months the wat 
losses decrease so that the runoff becomes generally higher. The areal 
ations of average water losses during any season do not seem to vary sig ni 
cantly. Thus, although the precipitation does not vary much from season t 
season the effects of temperature and topography cause large seasonal flue 
tuations in runoff. x 

Extreme high and low monthly mean precipitation and temperatures at | 
stations are shown on Figure 11. A range of one standard deviation on eitt 
side of the average precipitation is also shown. The stations at McKeever 
Alfred, N.Y., and St. Catharines, Apsley, and Kingston, Ontario, are con- 
sidered to be representative of the Adirondack plateau, the Allegheny plate 
the western end of Lake Ontario, the Precambrian shield and the eastern € 
of Lake Ontario respectively. . 

It can be seen that, with the exception of Alfred, the average precipi’ i 
does not vary much from month to month. The higher summer average pr 
cipitation at Alfred indicates a somewhat different climatic condition at thi 
extreme southern end of the local drainage basin. The greatest variations 
precipitation generally appear during late summer or autumn, possibly du 
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nderstorm activity or the occasional tropical storm which may cross the 
in. Although only the maximum and minimum monthly mean temperatures 
shown, the degree of variation of these quantities is an indication of the 
jation that can take place within the month. The winter period exhibits the 
est variations in mean temperatures. The minimum and maximum month- 
emperatures for all the stations follow the same trend. The similarity in 
perature variation between two widely separated stations i.e. McKeever 
Apsley may be noted. Apparently the effect of increased elevation at 
Keever balances the effect of the northerly location of Apsley. 
For a more detailed study of the seasonal distribution and variation of the 
d runoff it is necessary to examine more closely the behaviour of the in- 
idual tributaries. On Figure 12 is presented a hydrograph of the mean 
nthly flow in inches for various tributaries and for the net local water 
ply. Also shown is the percentage of the mean annual flow for each month. 
> much higher unit runoff from the Adirondack plateau into the Oswegatchie 
Black Rivers is immediately obvious. The lack of a deep soil cover and 
rage reservoirs in the Moira and Napanee drainage basins has a marked 
»ct on the distribution of runoff. Much higher percentages of the annual 
y occur in the spring and much lower percentages in the late summer than 
h the other tributaries. Moreover, on these tributaries there is a smaller 
centage of the annual runoff during the winter which would indicate a 
ater accumulation of snowcover. This would tend to increase the spring 
off. Throughout the greater part of the local drainage basin of Lake 
ario the peak flow generally occurs in April although on the Genessee and 
dit Rivers in the western part of the basin the spring season seems to be 
re advanced. The Trent and Oswego Rivers have the most uniform dis- 
ution of average monthly flows because of their large size and of the ef- 
ts of regulation. With the exception of the Genessee and Credit Rivers the 
rage flow for January for all the tributaries is higher than for either 
‘ember or February. This would give some credence to the local tradi- 
. of a “January thaw.” The runoff pattern for all the tributaries is, in 
te of these differences, generally similar to that for the net local water 
ply. The difference can generally be attributed to the effect of 
cipitation-evaporation balance on the lake. The maximum and minimum 
vs for each month are compared to the mean flows for each month on 
ure 13. The monthly distribution of the standard deviation is also shown. 
. figure illustrates that the range of variation is smallest in the summer 
aths and greatest in the winter and spring. The reason for the small 
ge in variation of runoff is that the evaporation and transpiration losses 
rease with the availability of water. The decrease in the effect of these 
ses in the fall allows a much greater range of runoff. The effect of large 
perature variations on the snow that has fallen and the small water los- 
tend to make the variability of runoff during the winter extremely high. 
. variability during the spring period results from the variation of the 
w cover at the end of winter as well as the spring precipitation. Generally 
standard deviation follows the same pattern as average monthly runoff. 
tt may be concluded from the foregoing that the variability of runoff is de- 
dent not only on the variation of precipitation but also on the variation of 
perature. A specific example in which the precipitation and temperature 
related to runoff is shown on Figure 14. Deviations in precipitation and 
perature at Belleville for the period 1941 to 1950 are compared graphical- 
yith the deviations in flow of the Moira River. Belleville is at the mouth of 
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the river and it is believed that the deviations of the temperature and pre¢ 
tation from average reflect conditions throughout the watershed with suffic 
accuracy for this study. 

Although there are many interesting occurrences shown on this plate on 
three periods will be examined closely, These periods have been chosen t 
illustrate the factors influencing the variability of runoff which have been ¢ 
cussed previously. The results of the examination can be summarized as 
follows: 


Year Month Devigtion from Average Remarks 
Precip: Temp: Runoff 
in: oF in: 
1947 May 40.98 00.80 Increases in the evaporation an¢ 
June 43.98 41.75 transpiration losses reduce the 
July 44.57 20,56 effects of very much above averg 
precipitation, 

st -1.49 40.18 Decreased evaporation and trans- 
eis at anal 1.41 20.18 piration losses allow the incres 
October @1.42 @1.52 rainfall of September and Octobe 
to be reflected in the runoff fo 
October, g 
1949= November #0.50 “4 -0.59 Accumulation of snow, bg 
1950 December -0.01 +4 00.58 Some snow melt. 4 
January #1.46 08 42.36 All snow melted, . 
February 41.47 =-2 20.48 Accumulation of snow. 4 
March -9.96 -4 -0.79 Accumulation of snow. ‘ 
April 20.20 41.40 Melting of snow accumulated in © 
February and March. > 


> 


It may be concluded from the foregoing analysis of runoff variability 
the local drainage basin of Lake Ontario lies in a transition zone with rest 
to winter and spring runoff. To the north of the basin the snow cover is co 
tinuous throughout the winter so that the winter runoff is usually low and u 
form and the spring runoff high and irregular. South of the basin, continud 
accumulation of snow is rare so that the runoff varies more directly with 
variations in precipitation. In the drainage basin of Lake Ontario itself, a 
either of these conditions can occur depending on temperature variations — 
during the winter months, It may be concluded that, because of the varia 
bility of temperature during the winter, the use of snow surveys for the pi 
diction of runoff may have a limited use only. : 

Although instantaneous and daily peak flows are of little importance in 
study of the hydrology of Lake Ontario, an indication of the flood potential 
the drainage basin is of interest. Figure 16 illustrates a plot of the peak 
flows of the gauged tributaries against the drainage area, The envelope 0 
peak flows and its equation are also shown on the plate. It is interesting 
note that the peak quarter-monthly net local water supply lies slightly abc 
the envelope of extreme flows for the tributaries. This may be explained 
the fact that the period of record represented by net local water supply is 
years which is much longer than the period for any of the tributaries of cI 
basin. The maximum quarter -monthly net local water supply occurred in 
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Precipitation-Evaporation Balance 
The contribution of the surface area of Lake Ontario to the net local wat 
supply is derived from the precipitation falling on the surface as modified | 
the evaporation from the surface. Since there is no direct method of meast 
ing either precipitation or evaporation on such large lake surfaces the 
precipitation-evaporation balance can only be found from a water budget stu 
Such a study has been made for the period from 1934 to 1952 inclusive. Thi 
total land runoff for each month has been derived by adding together all the 
measured inflows, with the exception of the Niagara River, and estimating t 
remainder of the land runoff on the basis of the known unit water yields. Th 
estimated land runoff was then subtracted from the net local water supply t 
obtain the precipitation-evaporation balance of the lake. Since it represent: 
the difference between two large quantities, the precipitation-evaporation — 
balance for any one month is subject to large percentage errors. However, 
the average for each month over 19 years shows the average distribution of 
the precipitation-evaporation balance with sufficient accuracy for the pur- 
poses of this study. 
The average distribution of the precipitation-evaporation balance is sho¥ 
on Figure 15. This figure would indicate that the average contribution of th 
surface area of Lake Ontario to the net local water supply is negligible and 
that the only effect of the lake area is to increase the net local water supply 
during the seasons of high land runoff and decrease it during the seasons of 
low land runoff. The average cumulative effects of the seasonal increases 
and decreases are substantial, being of the order of 60,000 cubic feet per 
second for one month or three quarters of a foot on the lake. } 
There is no method of accurately analysing the two components of the 7 
precipitation-evaporation balance. However, by using the precipitation : 
records of stations near the shore of the lake the apparent precipitation on 
the lake has been determined. By subtracting this from the precipitation-— 
evaporation balance, the apparent evaporation was deduced. 
It may be noted from Figure 15 that the average apparent precipitation _ 
varies very little throughout the year. The average apparent evaporation, — 
however, varies considerably reaching its maximum value in October and 


minimum value in June. It is relatively uniform from August through Deceé! 
ber. 


Outflows and Water Levels of Lake Ontario 


The hydrologic factors influencing the water supply to Lake Ontario hav 
been given careful consideration. In the past the changes in water supply 
have resulted in changes in water levels and outflows from Lake Ontario bi 
because of the large storage capacity of the lake and the constricted outlet 
conditions the changes have been modified and delayed. Thus, although a 
range of one foot in the elevation of Lake Ontario can absorb 80,000 c.f.s.- 
for one month a one foot change in the water level only changes the outflov 
by about 20,000 c.f.s. 

The control section of Lake Ontario is in the St. Lawrence River down: 
stream from Prescott, Ontario. Previous to the construction of the St. 


the river. During the past century, however, slight changes in this cont 0 
section have resulted in overall changes in the levels of Lake Ontario, Th 
water levels which will be used in further discussion are the recorded val 
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ited for these changes and for the diversions. During the winter, ice is 

ea at and above the control section restricting the outflow and causing 

‘ases in the water levels of the lake. Temporal changes in the outflows 

vater levels of Lake Ontario are caused by winds and seiches. Because 

; location and greater depth, the effects of winds and seiches on Lake 

“io are not as great as on Lake Erie. 

) obtain a more detailed knowledge of the interaction of the outflows from 
Erie and the net local water supplies and their effects on the water 

s of Lake Ontario in the past, a further investigation was made. It was 

assumed that the outflows from Lake Erie remained constant at their 

age value but that the net local water supplies varied as they did in the 
These hypothetical net total water supplies were then routed through 

Ontario using standard routing procedures and taking into account the 

1 effects of ice formation that had occurred. The same procedure was 

with another set of hypothetical supplies which were prepared by as- 

ng that the net local water supplies remained constant at their average 
but that the outflows from Lake Erie varied as they had in the past. The 

ts obtained were thus directly comparable with the recorded adjusted 

~ levels of the lake and are summarized in the following tabulation: 


Monthly Mean Water Levels and Ranges of Stage 


Under Outflows from Lake Net Local Water Supplies 
Recorded Erie constant Net constant Outflows from 
Adjusted local Water Supplies Lake Erie vary as in 

, ___.___ Conditions vary as in the past the_past 

um Water 

. of Lake 

‘io 248.493 248.67 247.86 

‘um Water 

of Lake 

Lo 242.84 244.11 244.03 

| ee ee 

um Range : 

age 6.09 4056 3.83 

ge Annual 

of Stage 1.75 1,81 0.62 


e marked decrease in the average annual range of stage, compared with 
ded adjusted conditions which would have occurred had the net local 
‘supplies been kept constant, indicates that the variations in the net 
water supply cause most of the annual variation in water levels. Con- 
ly, the slight increase in the average annual range of stage that would 
»ecurred had the outflows from Lake Erie been kept constant shows that 
‘agonal variation of the outflows from Lake Erie is not in phase with 

al variations in the net local water supply. The effect of the long term 
ions in the outflow from Lake Erie is shown by the decrease in the 

um range of stage with respect to recorded adjusted conditions that 
have occurred had these been kept constant. The relatively greater de- 
in maximum range with constant net local water supplies shows that 
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the variations of the net local water supplies have had a greater effect on th 
variation of water levels than have the variations in the outflow from Lake 
Erie. Since the variation of annual average water supplies is smaller for 
net local water supplies than for the outflows from Lake Erie, it can be con 
cluded that the seasonal variations of the former are of considerable impor 
tance in producing the maximum range of stage. 

It is difficult to evaluate the effect of the variations of the different sours 
of water supply on the regulation of Lake Ontario. Critical water supply co 
ditions can vary in length from three months to three years. In the studies 
which have taken place it has been found that the critical recorded adjusted 
levels can be used as an indication of critical conditions for regulation. Wit 
this criterion and from the results of the foregoing routing study, it may be 
concluded that the variations in the net local water supply will cause most 0 
the difficulties in regulating the outflows and water levels of the lake. It ca 
also be concluded that regulation of the upper Great Lakes which would mer 
ly smooth out the seasonal fluctuations in the outflows from Lake Erie woul 
not assist in the regulation of Lake Ontario although if they could be change 
radically to offset the seasonal variations of the net local water supply som 
improvement might result. A system of regulation which could smooth out 
the long term variations in the outflows from Lake Erie would benefit con- 
siderably the interests concerned with the regulation of Lake Ontario. Con: 
sideration of the vast volume of water involved in these long term variation 
indicates that an appreciable improvement would be very difficult to achieve 


Forecasts of Net Total Water Supply 


The regular seasonal variations in the outflows from Lake Erie and the 
delayed response of the long term variations in outflow to the variations in 
precipitation on the drainage basin of the upper Great Lakes indicate that it 
is possible to forecast the outflows from Lake Erie for a long period in ad: 
vance with a considerable degree of accuracy. However, with the extreme 
variability of net local water supply and its dependence on variations of tem 
perature as well as of precipitation, accurate long term forecasts of this 
component could not be expected. With the net local water supply such an 
important factor in the regulation of Lake Ontario, it would seem impracti¢ 
ble to devise a method of regulation that was based entirely on forecasts. _ 
Forecasts may be of use, however, as a supplement to regulation. Such 
forecasts may be of short or long term. The accuracy of short-term fore- 
casts of water supply would depend upon the accuracy of temperature and 
precipitation forecasts. Long term indications of future supplies may be di 
termined on a Statistical basis utilizing the observed, but not invariable, té 
dency of the annual mean outflows of Lake Erie and the annual mean net lo’ 
water supplies to be high or low at the same time. Because forecasts of <} 
ficient accuracy would be of too short a period and the long term indicatio 
would not be sufficiently accurate, the use of these techniques in the regul: 
tion of Lake Ontario would require limitations on the use of storage estab- 
lished by testing over the period of record, 


SUMMARY 


From the study of the hydrology of Lake Ontario which has been presen 
it has been determined that the water supply to Lake Ontario has two prin 
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onents. The outflows from Lake Erie provide a very large and stable 
flow to Lake Ontario which increases or decreases predictably in re- 

e to long term changes in hydrologic conditions on the upper Great 

3. Although the ratio of the variations is small in terms of mean out- 
the long term variations do represent an extremely large volume of 

. The net local water supply, although small in proportion to the out- 
rom Lake Erie contributes more to the variation of both net total water 
y and water levels of Lake Ontario, Although the average seasonal dis- 
ion of precipitation for the drainage basin is very uniform other factors 
the net local water supplies to vary widely. These factors are the vari- 
in the areal distribution of precipitation on the local drainage basin, the 
mal variation of water losses, the effect of the winter temperature vari- 

; on the accumulation of snow and the effect of the seasonal variations of 
‘ecipitation-evaporation balance on the lake which intensify the effect of 
nal variations of the land runoff. In addition to yielding an understanding 
‘water supplies to Lake Ontario the results of the study have indicated 
ne possibility of making accurate long range forecasts of water supply is 
* remote. It has also been indicated in the study that the regulation of 
per Great Lakes in such a manner as to obtain an appreciable benefit to 
ists on Lake Ontario would be very difficult. 
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RESISTANCE EXPERIMENTS IN A TRIANGULAR CHANNEL? 
Ralph W. Powell,/ M. ASCE and Chesley J. Posey,2 M. ASCE 


SYNOPSIS 


, 400-ft adjustable-slope triangular flume was tested smooth and also 
thened with small rectangular battens at various spacings. A comparison 
sade of Manning’s formula with some more recent suggestions, and the 
city distribution, secondary currents, and other phenomena are reported 
riefly. 


INTRODUCTION 


‘his paper reports the first phases of a comprehensive study of resistance 
ow in open channels being conducted at the Rocky Mountain Hydraulic 
oratory, Allenspark, Colorado since 1954. The major part of the support 
come from grants from the National Science Foundation, the first two of 
th were received through the Colorado State University with Maurice L. 
rtson,3 M. ASCE as co-director. Support has also come from the Water 
purces Division of the United States Geological Survey, the State Univer- 
of Iowa, the Ohio State University, the University of Kansas, and from 
funds of the Laboratory itself. 

the triangular shape of flume was adopted to maintain geometrical simi- 
'y at all depths. Several others, such as Rogers, (5) varwich, (9) Owen, (13) 
Moeller-Hartmann(15) have tested the same shape, but their work was 

: indoors where the length was necessarily limited. The authors feel that 


‘a short flume it is very difficult to tell when the flow is at normal depth, 
e the upstream and downstream depths may be the same as nearly as can 


neasured, for quite a range of depths. (This point is discussed more fully 
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later.) Therefore this flume was built outdoors and made 400 feet long. It is 
believed to be the longest channel systematically tested since Darcy and 
Bazin(1) and quite certainly the longest one of variable slope. As described 
later, the slope of this flume could be adjusted from nearly zero to slightly 
over 1.5%. 


Design of the Tests 


The purpose of an open-channel friction formula is to express the relation- 
ship between the channel shape and roughness, the discharge or velocity, the 
depth of flow, and the bottom slope for the condition known as steady uniform 
flow. Since the yearly cost of open channel construction in the United States 
undoubtedly totals many millions of dollars it behooves us to make sure that 
we evaluate this relationship as accurately as possible, introducing minor 
variables such as the temperature, if necessary, and even increasing the 
complexity of the computation if that would yield a better result. 

Despite the wide use of open-channel friction formulas, it should be noted — 
that the condition of steady uniform flow is seldom encountered, either in the 
field or in the laboratory. The depth of flow for steady uniform flow, known 
as the normal depth, is an abstraction, an extremely useful tool without which 
flow profiles in open channels could not be predicted with any satisfactory 
degree of reliability. ‘ 

In determining normal depths in the laboratory, however, it is essential _ 
that steady uniform flow be attained just as closely as possible, Unsteady 
flow may differ imperceptibly from steady flow in all other respects but the ‘ 
depth of flow. In the present experiments, not only was the discharge intro- — 
duced into the flume held as constant as possible, but the water-surface pro- 
files were repeated as many times as necessary until they reached equilibri- 
um, as far as that could be determined by the precision of setting and readi 
the point gages. For some conditions of flow, this would take twenty minutes” 
or longer. 

After the flow is steady, it is even more difficult to be sure that it is uni-— 
form. Backwater curves of the M1 and M2 type, corresponding to slightly dif- 
ferent tailgate settings, may be very nearly parallel to the bottom, and yet — 
have an appreciably different depth of flow. That the flow is uniform cannot 
safely be determined by reading gages at each end of the flume, for the error 
in setting and reading only two gages, especially if the flume is short, may 
seem to indicate uniform flow whereas actually the flow may be well along an 
M1 or M2 curve, and the apparently uniform depth far from the desired nor- 
mal depth. 

In lesser degree, the same thing is true for flow on steep slopes. Althoug 
the S2 and S3 backwater curves are not nearly as long as the M1 and M2 
curves, there are apt to be standing waves on the surface that make accurate 
setting of point gages more difficult. y 

_ In determining the normal depths, therefore, profiles were obtained both 
below and above the normal depth, each profile for a steady flow condition. 
Fig. 1 shows a typical set of profiles for flow on a mild slope, and Fig. 2a — 
typical set of profiles for flow on a steep slope. The different curves in Fig. 
are obtained from different settings of the tail gate at the lower end of the - 
flume, and those of Fig. 2 from different settings at the entrance. Notice the 
in the 9:52 a.m. series of August 1, 1957 the normal depth of 2.475 feetis 
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Normal depth not obtainable 
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irly well determined by an M1 run, an M2 run, and two runs very near the 
9rmal depth. This is about the maximum depth that could be obtained in the 
ume at the setting of that date, for no normal depth could be determined for 
¢ Slightly greater discharge of 9:26 a. m., July 15, 1957. 

The 1:54 - 3:06 run of June 26, 1957, shown in Fig. 2, gives a normal depth 
1,45 feet which is about as great as could be determined at this roughness 
id slope. Of course, for the smaller depths on a steep slope, the S2 and S3 
ipves are comparatively short, the normal depth is soon reached, and it was 
sually sufficient to measure only the S2 curves starting from the critical 
pth. 

When a test is being started in the field, with the depth for uniform flow 
it determined, the experimenter does not know whether the flow is tranquil 
‘shooting. To be sure, if it is well within either range the characteristics 
the flow will immediately show which one it is. Flow on the critical boun- 
Try range, however, does not have any well-defined characteristic except 
riace instability. Hence, whenever there was any question, profiles were 
ken with different headgate settings and also with different tailgate settings. 
lis was not done with intent to explore the interesting backwater curves in 


iS vicinity, but rather to establish the depth of uniform flow as reliably as 
ssible. 


Description of the Apparatus 


Water is diverted from the North St. Vrain Creek into a forebay, through 
netal flume, and into a weir box 42 feet long, 8 feet wide and 7.5 feet deep. 
e Fig. 3. The flow can be controlled at the point of diversion and there is 
emergency spillway for protection against high water in the forebay at the 
le of the metal flume, which latter can be barricaded against extreme high 
ter. Flow is stopped by inserting the needles at the point of diversion and 
sting leakage over the emergency spillway. It is started again by closing 
> emergency spillway and removing enough needles to obtain approximately 
‘desired discharge. Adjustment to hold the flow steady is made by moving 
needles in a wasteway in the side of the weir box. Only a small amount 
discharge was permitted at this station, the operator of which could see a 
. Of temporary hook gages fastened inside the weir box, and adjusted to 
yw deviations from the desired water level of only a few thousandths of a 
t. 
Discharge from the weir spilled into a deep tumble bay from which it 
ered the variable-slope flume over a curved galvanized metal entrance 
‘tion. The radius of curvature was about 2 feet and the bottom of the 
ible bay was 6 feet below the bottom of the flume. 

The variable-slope flume had a triangular cross-section with central angle 
t did not vary more than 0.0025 radian from ninety degrees. The bisector 
his angle was set vertical; in other winds, each side was set at 45° with 
horizontal. The walls of the flume were made of four-foot by eight-foot 
els of one-half inch fir plywood, exterior grade, stiffened at one-foot cen- 
s by fir two-by-fours. The two-by-fours along the bottom edge were off- 

Slightly so that the reinforced panels nested and formed a tight joint, 
ch was held together by 1 /4 inch bolts at one-foot centers. When the tops 
‘e fastened the correct distance apart by struts and heavy wire ties spaced 
‘wo foot eight inch centers, a rigid eight-foot section of flume was formed 
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hich could be joined to its neighbors by splice plates. Fabrication was with 
hijlips head wood screws. 

The flume was supported on heavy timber bents at twenty-four foot centers. 
he upper surfaces of the supporting timbers were set to a common plane at 
1 inclination of 45° with the vertical. Bearing blocks under the reaction stif- 
ners of the “bent-plate” girder formed by the flume sections rode up and 
ywn the 45° inclines, as controlled by jacks resting on the inclines. The 
ume was, therefore, approximately to a true grade whenever it was in align- 
ent. Final adjustment perpendicular to the incline was made by driving 
lingles under the bearing blocks. Before water was introduced additional 
saring blocks were inserted under the flume at the third-points between the 
Clined supports, and wedged tightly above vertical logs carrying the load 
ywn to wooden plates set directly on the ground. To prevent tipping, the ends 
‘the stiffeners lying over the inclined supports were wired down. Since the 
yper flanges of the bent-plate girder did not have sufficient stiffness perpen- 
Cular to the plane of the side, they had to be braced at the third points. This 
sO prevented the flume from tipping. 

Setting the flume accurately to a new slope required three days at first, but 
ith improved methods the time was cut down to little more than two days. 
ne flume was quite stiff, and raising one jack only half an inch would lift it 
om supports 24 to 48 or even 72 feet on each side without causing distress 
‘noticeable distortions at the joints lifted. It is hoped that when the hydrau- 
2 tests have been completed, the flume will still be in good enough condition 

permit structural tests of the bent-plate girder to be made. 

The plywood panels were given a coat of penetrating sealer and one coat of 
rest green paint when fabricated at the hydraulic laboratory of Colorado 
ate University in Fort Collins early in the summer of 1955. When it was 
‘sired to start testing in 1956, it was found that the inner surface of the 
ime, which had been exposed to the weather, was beginning to check, lose 
int, and even suffer some ply separation. Since the bearing blocks under 
e intermediate stiffeners were not of oak, as had been specified but found 

be unavailable, and the intermediate splice plates not completely fastened 
e to a shortage of screws, the flume had developed a “sway-backed” con- 
tion with sags at the intermediate points ranging from 0.01 to 0.03 foot. 
1ere were also slight local bulges (less than 0.02 foot) over the main sup- 
rts, caused by placing wedges under the stiffeners instead of under the 
aring blocks. 

While the flume was being realigned and the missing screws inserted, it 
came apparent that the painted surface was suffering rapid deterioration. 
though the altitude is 8000 feet above sea level, day-time temperatures in 
2 intense sunlight in the flume went as high as 130°F, while frost formed 
the surface at night. Several solutions to the problem were tested. That 
ich was finally adopted was to paint the entire interior surface with two 
ats of aluminum asphalt paint, using rollers and brushes. The resulting 
rface was fairly smooth, daytime temperatures were much lower, and the 
pid deterioration was stopped. Cracks in the surface and between panels 
re sealed with asphalt-asbestos caulking compound, and, after painting, 
vered with black vinyl chloride adhesive tape. The tape alone would have 
en sufficient to seal the cracks, and perhaps better by itself as the asphalt 
i not have as much flexibility, and weakened the adhesion of the tape. 


After a preliminary test with a maximum load of water, the bolts fastening 


2 bottom together and the wires holding the top struts tight were tightened. 
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ing subsequent tests no leaks greater than a few slow drips were observed, 

n July of 1957 a new method of bracing the flume that was being tried in 

ffort to cut down the time required for changing slope overstressed the 

, and during a run with the flume nearly full, 160 feet of it collapsed. The 

its of the design became apparent when it took only a little over two weeks 

ebuild, using mostly salvaged material. During the rebuilding greatly im- 

red methods of fabrication were developed, and a simple and fast method 

racing and aligning the flume was discovered which solved that hitherto 
ersome problem. 

pecially-built point gages were installed at twenty-four foot spacing along 

centerline of the flume. These gages were held vertical by adjustable 

orts fitted with grooved guides and a spring so that the gage could easily 

et with one hand. 

1 order to read the point gages which were set to indicate the water sur- 
and special level rods which were used to measure the bottom elevations, 

2 levels were kept set up alongside the flume. Intermediate bench marks 

> used to establish the instrument heights accurately, and in addition the 

- gages midway between two instruments were usually read from both. 

he “headgate” of the variable-slope flume, which had to be used when the 

city was supercritical, consisted of a vertical bulkhead with rounded bot- 

Surface made from a log and sometimes from two or three additional logs 

2d so as to decrease further the entrance depth. It was fastened in at the 

r end of the flume, just at the end of the curved entrance section. Due to 

issymmetry of the tumble bay the surface waves below the headgate were 
unsymmetrical. Usually they could be reduced by fastening a baffle 

1 at a location below the headgate that had to be determined by trial. The 

mmetrical surface waves seemed to disappear within a hundred feet of the 

ince in even the worst cases. 

1e “tailgate” of the variable-slope flume was simply a set of needles 

1 could be fixed in place as desired. 

ijustment of the headgate was required when the velocity was supercriti- 

nd adjustment of the tailgate when it was subcritical. In general, it was 

ractice to wait until the flow became steady, as indicated by the point 

3, before recording any readings, and to adjust both head and tail condi- 
if there was any doubt as to which should govern. Occasionaily it was 

, upon repeating a set of profile readings, that changes had taken place, 

ating that the flow had still not become quite steady. The readings were 

ited until the slow changes disappeared. 


Calculation of Flows 


general flows of over 4 cfs were measured with a full-width sharp- 

2d weir and those smaller with a 90° V-notch weir. 

e full-width weir was of stainless steel 24 in wide and 7.948 ft long. It 
bout 0.25 in thick but the upper edge was beveled at 30° to the vertical 

t for a flat top machined to 3/32 in width the same as Schoder and 
r’s(3) Series D to P. At first every time the zero of the gage was de- 
ned, levels were taken at each foot along the crest. It was found that 
rth end was 0.004 ft higher than the south end, and the average of all 
adings was used to compute the discharges. But it was found that this 
ze was apt to be the same as the height at the middle, so that sometimes 
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only this was taken. The upstream face was as nearly vertical as it could b 
set with a carpenter’s level, and the weir plate was held in place by counter 
sunk screws, so that the entire surface was smooth, The end of the weir bo: 
below the stainless steel weir was of concrete and quite smooth and vertical 
but in the plane of the downstream face of the weir plate. 

The 90° V-notch was cut in a sheet of 1/4 inch aluminum-painted temper 
masonite the same length as the full-width weir and 4 ft wide, which was 
placed with its lower part resting against the stainless steel weir plate and 
its upper part supported by a wooden frame. The bottom of the notch was 
about 0.3 in above the crest of the steel weir. The edge of the V was formet 
to the same shape as the full-width weir. It was found to keep its sharp edg 
very well. The ends and bottom were sealed with tape so as to make leakag 
negligible. 

The head on the weir was read on a hook gage attached to a plank which 
was spiked to two trees beside the weir box. The hook gage well was a lucil 
cylinder about 5 in, in diameter which was connected by about 4 ft of 3/4-in 
transparent hose to a short piece of half-inch pipe which projected from the 
side of the weir box at about 12 ft upstream from the weir. The inner endo 
this short pipe was filed to a plane and was flush with the inside of the weir 
box. The hook gage well was drained and refilled from the weir box wheney 
it was found that the water in the well was as much as 1.5°C warmer than th 
water flowing over the weir. The “zero” of the gage was determined by set 
the hook to the cross hair of an engineer’s level, and then reading a rod hell 
on the weir crest as mentioned above. For the V-notch a special target was 
constructed which could be placed in the notch with contact along both edges 
rather than just at the vertex. ; 

When a run was being made the hook gage was read every few minutes a! 
recorded with the time. These readings were plotted against time, and fror 
this graph the average head for each run was determined. The median vari 
tion between maximum and minimum gage reading during a run was about q 
0.006 ft so that the error in the mean gage reading probably averaged about 
0.001 ft. But there may have been another 0.001 ft error in the “zero readi 

The flows from the V-notch were computed by the formula given by Leng 
which reduces to 3 
2.994 


fe % 
Q=Kh with K = 2.395 + R165 w 170 


but not less than 2.485, Ris g0-5y1.5/y and Wis pgh2/g where a is the 
surface tension in lb/ft. Since the laboratory is located at an elevation of 
8000 ft above sea level, the value of g is less than at most laboratories. — 
Upon the advice of Professor Weikko Aleksantri Heisknen of the Ohio Stai 
University, it was taken as 32.136 ft/sec. 

It was originally planned to compute the flows over the full-width weir 
Schoder and Turner’s formula (D).(3) But when a velocity traverse was 
in the weir box opposite the hook gage piezometer to determine Vg and Vt 
was found that the Vp was quite a bit larger than Va; that is, the highest a 
proach velocity was near the bottom of the weir box. Schoder and Turne1 
data show that when Vp is larger than Vg, the flow is about the same as 
by Francis’ simple formula, or sometimes even less, while formula (D) — 
always makes it more. As described elsewhere, (17) a restudy of Schoder 
and Turner’s data was made, leading to the formula 4 
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Q = (3.08 + 0.44 h/P + 0.15 r)(b - 0.003)(h + 0. 003)!+5 


For our weir, b was 7.948 ft and P was 5.5 ft at the weir in 1956 and 1957, 
| this value was used for runs 1 to 108. However, the floor was not level 
| the average depth of the floor opposite the piezometer connection below 
weir crest was only 5.24 ft. Perhaps this should have been used for P for 
se runs but the change in Q would have been a very small fraction of one 
cent. Before work was started in 1958 a new floor was installed with 

5.2 ft at the piezometer connection and 5.0 ft at the weir face. P was 
en as 5.2 for runs 109 on. 
In the formula, r is the ratio of the velocity of approach above the level of 
crest to that below. Of course this was not measured for each run. From 
verses at five different heads and from Schoder and Turner’s results, it 
3 decided that r could be taken as 0.435 h + 0.870. (This was for 1957 and 
8 after additional baffling had been installed at the upper end of the weir 
. For the data of 1956 r was taken as 0.435 h + 0.199. This is different 
a the formula published in (16), so that the full-width discharges given in 
- publication have been revised.) All this leads to the following formulas 
discharge from the full-width weir: 


1956 Q = 7.945 (3.110 + 0.145h) (h + 0.003)1.5 
1957 Q = 7.945 (3.210 + 0.145h) (h + 0.003)1-5 
1958 Q = 7.945 (3.210 + 0.150h) (h + 0.003)1-5 


Mr. R. E. Glover of the U.S.G.S. made some diffusion measurements by 
salt velocity method during the summer of 1956 and kindly supplied us 
1 two values of the discharge which had also been measured by the full- 
th weir. The results are as follows: 


Test No, Qas measured by Difference 
USS; Gao; RoM, HJ. Ro MSH, =U. SiGe, 
1 4.85+.02 4. 88 0.03 cfs = 0.6% 
4 10.08+.11 9.96 0.12 cfs = -1,2% 


“arrent meter gagings taken in the weir box to find r (as noted above) 
2 the following rough checks of the weir formula: 


Date Value of Qin cfs Difference 
by weir formula by gaging 
ig. 24, 1956 17.19 17.17 Peo pape ioe 
ig. 7, 1957 10.22 11.49 41,27 cfs 12.4% 
#8, 1957 16,57 17. 64 +1.07 cfs 6.5% 
fi. 3, 1957 13.50 14.06 40.56 cis 4.2% 
pt. , 1957 19053 19.6 9e +0.06 cfs 0.3% 


pt. , 1957 24.71 24,11 -0.60 cfs 2.4% 
‘" 
> 
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Pitot tube traverses taken in the flume gave the following rough checks: 


Date Value of & Value of Qin cfs Difference 
by weir formula by traverse 


Aug. 24, 1956 1.05 6, 68 6, 65 +0.03 cfs 0.45! 
Aug. 10, 1957 09 14,51 14, 88 +0.37 cfs 2,55!) 
Aug. 15, 1957 LO 14, 44 14, 38 -0.06 cfs 0.41 
Aug. 23, 1957 1,03 17.40 V7579 +0.39 cfs 2.24 
Aug. 27, 1957 1.06 LHS90 L776 -0.14 cfs 0.78 
Aug. 30, 1957 1.04 23.72 PAG YI he -0.33 cfs 1.39 


Determination of Bottom Slopes 


For each slope and roughness a number of runs were made covering a Wi 
range of discharges and a similarly wide range of depths of water in the 
flume. The weight of water in the flume varied from a few pounds per foot t 
as much as five hundred. Since the supporting system was necessarily sub- 
ject to deformations from the varying load, and the heights of the main and ~ 
intermediate supports differed along the length of the flume, it was deemed 
necessary to determine the bottom-profile of the flume under the conditions 
of the test runs, with the water flowing as near the normal depth as possible 
To permit this, special steel level rods were made which could be held in 
position even during the high velocity runs and read from the three levels sé 
up along the length of the flume. Readings were taken at eight foot intervals 
along the bottom. The rod probably interfered slightly with the water -surfa 
profiles, but water surface-profiles were not taken while the rod was in pla 
and the rod was in the flume only part of the time. For some positions of tt 
flume and levels, the supporting bents or other obstacles interfered with the 
line of sight, in which case the readings were simply omitted unless it was_ 
possible to include them by making an additional level setup without too oa 
delay. :. 

Because of the time required, bottom-profile data were not taken for eve 
water-surface profile run. Usually bottom-surface profile data were taken 
corresponding to the highest and lowest runs of a given series. The differ - 
ence between the mean bottom slopes and elevations determined at these ex 
tremes were usually negligible or very slight, so that the values for inter- 
mediate runs could be assumed to be the same as that of the nearer run, OF 
interpolated between. 

The profile data were plotted in pencil on matte-surfaced acetate sheets. 
using a scale of 20 feet to the inch horizontally and an adequate vertical s¢ 
depending upon the slope. After each of the observed elevations had been ii 
in, using small circles centered on the plotted points, the straight line of b 
fit was determined by shifting a straight line drawn in ink on a separate sh 
around under the acetate until its fit was judged to be as good as possible, 
Then a fine pencil line was drawn on the acetate above this inked line, us in 
a steel straightedge, and its elevation at each end of the flume determined, 
permitting computation of the slope. . 

In order to gain some idea of the accuracy of this graphical method of ¢ 
termining bottom slopes Mr. Chintu Lai(19) made a study of the results — 
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ined by 50 individuals from identical prints of the graph of data of 3:00 

no, August 31, 1956 and compared it with the least Squares value computed 

m the numerical data. He found the mean of the graphically-determined 

pes to be 0.008000, with standard deviation of 0.000030. Least squares 
ution gave a slope of 0.007998, with standard deviation (corresponding to 
Scatter of the data) of 0.0000144. It seems reasonable to conclude that 
Slopes were measured to within one-half of one per cent, 


Determination of Normal Depths 


Water-surface elevations obtained from level readings of the point gages 
unted at twenty-four foot intervals along the flume were plotted on the same 
tate sheets as the bottom-profiles, and after the plotted points were inked 
line best representing the normal depth was drawn by a similar technique, 
ing care to make it parallel to the appropriate bottom-profile line. For the 
6 tests, this was done in the drawing room after the runs had been com- 
ced, but for the 1957 and 1958 tests the pencilled plots were made at the 

ne time the observations were taken, the slope and normal depth computa- 
iS being made later. All of these drawings for the 1956 tests are repro- 
ed at reduced size in the R.M.H.L.’s Report No. 21, and it is expected that 
of the 1957 and 1958 profile drawings will be reproduced in a forthcoming 
ort. The results are summarized in Table I. 


Brief History of Search for a Formula for Channel Resistance 


The first formula for the resistance to flow in open channels that is still 
ognized is V=CVRS , proposed by Antione Chezy in 1775. He evidently 
90Sed that C would be constant for any given channel, but this proved not 

e the case. In 1869 Ganguillet and Kutter presented a resume of practical - 
ll of the measurements of open channel resistance that had been made up 
iat time and from them derived the formula which bears their name. It 
very widely used for many years, but its cumbersomeness and the fact 

it is not dimensionally homogeneous has led to its replacement by other 
qulas. 


1 1865 Bazin had proposed that the Chezy C could be expressed as 


a+b/R but in 1897 he suggested instead C = tnege ale - Neither of 


e formulas seem to have been used in this country. In 1868 Philippe 
yard Gauckler, a French engineer, proposed the formula 


ged oink bowsgag 


or some reason he limited it to slopes less than 0.0007. This formula 

is to have attracted no attention at the time. eee) proposed the same 

ula in 1881, without the limitation, and Manning(2) repeated it in 1889 

ut acknowledgement. But he was not satisfied with the formula and. 

ad recommended a much more complicated one. In 1894 Crimp and 

es suggested V = 124 R2/3 1/2 for sewers and water mains, a formula 
1is still used in England. In 1899 C. H. Tutton proposed to the Boston ~ 


neering Societies the equation V = 138 R2/3 51/2 with n to be the same 
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as in Kutter’s formula. Although this is obviously just as satisfactory as whe 
we today call the Manning formula, it seems to have never attracted attention 
It was not until books published by R. B. Buckley in 1911 and P. A, M. Parker 
in 1913 gave the formula 


0. 67 <0. 50 
y --1.486R-°'S 


and for some strange reason called it Manning’s formula, that it began to dis 
place all other formulas for estimating channel resistance. 

After the work of Prandtl, von Karman and Nikuradse had put a new light 
on the resistance to flow in pipes, and had shown the important difference be: 
tween smooth pipe and rough pipe flow, it was natural to ask if channels coult 
not be treated similarly. Probably the first to propose formulas on this basi 
was Keulegan.(4) He was soon followed by Johnson, (7) Powell, (8, 10, 11) 
Kirschmer, (9) Robertson and Albertson(12) and others. However, none of 
these developed formulas that are satisfactory for universal use. In fact it 
became obvious that when a channel acts as smooth (that is any roughness 
which is present is entirely imbedded in the laminar boundary layer) the law 
of resistance must be different than it is when the boundary roughness 
governs. 

There has been a tendency the last few years among the theoriticians to 
treat open channels as equivalent to a pipe with diameter four times the hy- 
draulic radius and then use the Weisbach formula and the Moody diagram for 
both pipes and channels. It is, of course, allowable to use the formula, but it 
must not be expected that f will be the same at the same Reynolds number, ~ 
and the same wall roughness. The free surface and secondary currents woul 
be expected to increase the resistance. 5 

A noteworthy treatment of the resistance to flow in rough channels is that 
of Morris.(14) Although the formulas he presents seem of little practical us 
and some of his generalizations too sweeping, he seems to have shown that — 
there are at least three types of rough channel (and pipe) resistance, name 
isolated-roughness, wake-interference, and quasi-smooth; and thatitis 
foolish to suppose that an “equivalent sand grain size” with the Nikuradse q 
rough pipe formula can be used to represent them all. 

a 


Discussion of Results 


The results of the present tests are given in Table 1. It must not, of 
course, be assumed that values are correct to the last significant figure. ; 
Suppose that about 21 cfs were flowing in the flume at a depth of 2.40 ft and 
slope of 0.00200. This would be measured by the full-width weir with a hea 
of about 0.850 ft. The measured head might be in error by 0.002 ft which © 
makes a percentage error of 0. 235%. Since Q varies as the 1.5 power of 
head this indicates a percentage error of 0.35% in Q, But since the form 
which is used in computing Q is probably not correct, it seems that the 
probable percentage error in Q might be as much as 1%. The slope may 
assumed to have a probable error of 0.00001 or 0.5%. The error in deter- 
mining yp might be as much as 0.005 ft or a probable percentage error of 
0.21%. Since f = 8S Rg/V2 = 90.9 S y,°/Q2, the probable percentage er 
in f would be 
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0.5% + (5% 0.21)% + (2x 1.0)" = 2.3%, 


But now suppose that 0.043 cfs were flowing at a depth of 0.250 ft with a 
2€ of 0.00045. With the same assumed error as before, the percentage 
or in the slope would be 2.22% and for the depth 2.00%. The flow would 
measured by the V-notch weir with a head of about 0.200 ft and probable 
centage error of 1.0%. Since the discharge varies as the 2.5 power of the 
d, this makes a probable percentage error in Q of 2.50%, or allowing for 
ertainty in the discharge formula, a total of say 3.0%. This makes the 
bable percentage error in 


fs 2.22% + (5 x 2.0)" + (2 x 3.0) = 11.9% 


3 indicates that we might expect percentage errors in f in the range from 
12%. When the flow is critical or shooting the error in measuring yn 
ht be even more than assumed above. Since C and n vary inversely as 


Square root of f, the expected percentage errors in them would be only 
it half as much, 


Results from Smooth Channel 


Vhen the values of f as found from the runs on the smooth flume (called 
n what follows) were compared with the values of f for a smooth pipe with 
same Reynolds number computed by Prandtl formula 


“= = 2 log (RY f) - 0.80 


Vt 


ed f, in what follows) the ratio W = {,/f) varied from 0.778 to 2.248. 
3 1-16 and 102-108 where the flow was definitely tranquil, with a maxi- 


epee eae Ww of 
1 Froude number (defined as ) of 0.85, gave an average W o 
= ( vg Yn/ 2 ; 


0. In other words the f was only 93% as much as for a smooth pipe at 
3ame Reynolds number. This is very surprising. This reduces the for- 
L 


Ts 2.074 log (IRVf ) - 0.797 
ranquil flow in smooth triangular channels. The average discrepancy of 
28 of f computed by this formula from the observed values was 7.1%. 

2 the square root of 1.071 is 1.035, this corresponds to an error of 3.5% 
orn, The average n was 0.0092 and the average discrepancy from this 
0.00032 or 3.5%. Therefore, the logarithmic formula has no advantage 
the Manning formula for representing these data. ’ 

irning now to the smooth channel data in which the flow was critical or 
ting (Froude number more than 0.85), the best fit was obtained by making 
15. This reduces to the equation 


1 = 1,866 log (IR Vf) - 0.803 
Vi 


HY 
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for shooting flow in smooth triangular channels. The average percentage 
discrepancy of values of f computed by this formula from the observed 
values was 11.2% which corresponds to an error of 5.5% in Vorn. The 
average n was 0.0094 and the average discrepancy from this was 6.8%, so 
that the Manning formula does not fit these data as well as the logarithmic. 
But the difference is not great, and neither formula gives satisfactory agree- 
ment. The explanation seems to be that when the flow is critical or shooting 
the surface is so irregular that there is too much error in measuring the 
depth. 


Results from Rough Channel 


The roughness used in the experiments described in this paper consisted 
of rectangular strips of redwood 3/16 in. thick by 3/8 in. wide tacked to the 
sides of the flume perpendicular to the flow. They stopped about 1/2 in. short 
of the bottom, so that the flume could drain dry after the day’s runs were 
finished, but it was felt that this would have very little effect on the resistance 
The strips had been sanded on the three exposed surfaces and the 3/16 in. 
dimension did not vary more than a few thousandths of an inch, The strips 
were used without any surface treatment except that they were soaked in wate 
before being applied to avoid their buckling out between brads when the water 
was turned on, Four longitudinal spacings were tested, 24 in. center to cente 
12 in., 6 in., and 3 in. This last required 3200 battens. The longitudinal 
at center to center in feet is represented by the symbol A as in Morris’ 
paper. 

In searching for a formula to represent the rough channel results it was © 
found necessary to divide the tests into two groups - those in which the flow 
was definitely tranquil (the Froude number less than 0.70), and those in whic¢l 
the flow was critical or shooting (Froude number 0.70), and those in which ; 
the flow was critical or shooting (Froude number 0.70 or more). It Should 
noted that these Froude numbers do not include the velocity distribution < 
factor a. ; 

For the tranquil case three formulas were tested. The first was Manni 
It gave the following values of n: 


Batten spacing Value of n Average discrepancy 
24 in. 0.0128 1.6% 
12 in, 0.0150 2.1% 
6 in. 0.0180 2.5% 
3 in, 0.0210 1.8% 


The weighted average discrepancy in values of n for these 52 runs was 2.0 
The same runs were computed on the basis of Nikuradse’s formula for 
rough pipe: 
1 
2 fF = log (r,/eE ) + 0.87 
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Substituting 2 R for To this reduces to 

1 

VE 


rom this formula € was computed for each run and the average value for 
uch spacing and the average percentage discrepancy in f computed by using 
is € was found as follows: 


= 2 log R- 2 log € + 2.342 


Batten spacing Average € Average % discrepancy 
24 in, 0.00459 ft nm 
12 in. 0.0126 ft 3,42 
6 in. 0. 0343 ft SON, 
5 in, 0. 0694 ft 6. 46 


ie weighted average percentage error in f was 4.40% which corresponds to 
2.18% average error inn or V. 

A third formula tried is based on Morris’ proposal for isolated-roughness 
w which for our case reduces to 


y = 1+ 67.2C) n/a 
with h = 3/16 in. = 1/64 ft, 


W =1+1.05 Cp/x 
lues of fp and W were tabulated for each run and values of W plotted against 
\ as in Fig. 5. It was found that the equation 


Y = 4.43- 0.145 (4.88-41/, )* 


resented the data fairly well. (This indicates that the maximum value of 

would occur when A: = 0.2049 ft = 2.46 in. or 13,1 times the batten height. 

is agrees quite well with Johnson’s results. See Fig. 14 on page 558 of 
erence (8). The formula is not expected to hold beyond the maximum.) 
lues of f were compute for all 52 runs by this formula and the percentage 
crepancy found, The mean was 4.53% which corresponds to 2.24% in n or 
Values of W and Cp by this formula are tabulated below: 


Batten spacing y Cy 
24 in. 1565 1,286 
12 in. Zita 1.190 
6 in. 323 1.062 


3 in. 4.32 0.790 
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The fact that Cp decreases with the spacing seems to indicate that even at 
these spacings the roughnesses are not entirely isolated, and that especially 
at the 3 in. spacing there is some wake interference. Also, the fact that none 
of the Cp’s approach the 1.9 value given by Morris, throws some doubt on his 
67.2 value. 

As the mean discrepancies in velocity resulting from a given Q and S by 
the three methods are 2.05%, 2.18%, and 2.24%, none is outstanding. Because 
of its simplicity and familiarity Manning’s would naturally be chosen, but it 
must be realized that the formula suggested here, based on Morris’ formula, 
is the only one of the three which gives a definite value in terms of batten 
height and spacing. 

Turning now to the runs with Froude number 0.70 or more, there are only 
three spacings, as no runs were made with a 3 in, spacing at a steep enough 
slope to give critical or shooting flow. The average values of Manning’s n 
and the average percentage discrepancy from the average for each spacing 
are as follows: 


Batten spacing Value of n Average discrepancy 
24 in. 0.0142 3.17% 
12 in. 0.0173 6.24% 
6 in. 0.0195 2.88% 


The weighted average discrepancy was 4.72%. : 
The same runs were computed on the basis of Nikuradse’s rough pipe for- 
mula with the following results: ’ 


4 
Batten spacing Average € Average discrepancy in f 2 
24 in, 0.0092 ft 5.52% , 
P 
12 in, 0.0279 ft 10.00% 
6 in, 0.0524 ft 7.70% 
The weighted average percentage discrepancy in f was 8.12% which a 


ponds to 3.98% error in n or V. 
An attempt was made to develop a formula connecting / and ) as was do 
for tranquil flow, but the scatter was so great that the result seemed worth- 
less. The Nikuradse type formula has a slight edge over the Manning formu 
for shooting flow on the basis of these tests, but the discrepancies from bott 
formulas are so large that the whole question remains in an unsatisfactory ~ 
state. In spite of the experimental errors present in these tests, it is be- — 
lieved that they show the resistance to shooting flow in triangular channels 
be definitely greater than for tranquil flow for corresponding conditions (i.¢ 
at the same relative roughness for rough channel flow and the same Reynolc 
number for smooth channel flow). This is contrary to the results found for 
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rectangular channels by Powell.(8) (As has been pointed out elsewhere the 
terminology in that paper was incorrect as Table 1 (a) included both tranquil 
and shooting flow and (b) was ultra-rapid flow. Ultra-rapid flow was not ob- 
tained in the tests reported here. To obtain it in a triangular channel requires 
a much higher Froude number than in a rectangular channel.) 

Increased resistance for flow at the higher Froude numbers can be ex- 
plained on the following logical basis. It was recorded by Montana(18) that 
the “capillary rise” of water along the sides of a rectangular flume increased 
greatly when the velocity in the flume became supercritical. The “fillet” at 
the water’s edge, normally very small, was apparently invaded by turbulence 
which so modified the balance of forces that the “fillet” enlarged to 3/4 in. on 
a side. 

This phenomenon undoubtedly occurred in the present tests and in addition, 
it was noted that water from standing and transitory waves caused by rough- 
ness and turbulence “rode up” the sloping side and fell back with considerably 
decreased momentum. In spite of the dryness of the climate at the laboratory, 
the flume was kept wet for a distance above the normal shore line, sometimes 
2 in. or more above it. High speed photographs show this condition and the 
retarded water in the waves. Clearly, this would add to the effective wetted 
perimeter, and cause an apparent increase in roughness if no additional 


wetted perimeter was taken into account. Possibly this effect could be evalu-— 
ated, 


Velocity Distribution ‘ 
; 
Pitot- tube velocity traverses were made for the smooth flume and for the : 
different batten roughnesses. Despite the fact that more than one reading ’ 
was seldom taken at a given point in the cross-section, the process was timely 
consuming. When two or three readings were taken at a fixed point, the vari- 
ability of the results obtained showed that in order to establish a reliable 
value for the average, perhaps as many as twenty readings would have to be 
taken. For this reason the traverses, a sample of which is shown in Fig. 6,9 
cannot be considered to show an accurate average velocity at each point, ; 
They do, however, provide a fairly reliable check on the discharge and value 
of a since the average of many observations is used in arriving at these 
values, 

Of more general interest are pitot measurements made to obtain the vari- 
ation of velocity perpendicular to the mid-point of the side of the flume. Close 
to the flume wall and in the middle of a long straight portion of the wetted 
perimeter, the variation of velocity ought to be free of any effect from the 
relatively distant water surface, or other segments of the cross-section 
boundary. If this relationship could be established in terms of roughness 
flow parameters, the mutual effect of the different boundary segments could 
be studied by means of interaction diagrams or other suitable techniques, thu 
yielding a more general means of evaluating the effect of cross-sectional 
shape. Use of the hydraulic radius implicitly assumes either constant or no 
mutual effect of velocity gradients generated at the different boundary seg- 
ments. While this causes little difficulty in comparing cross-sections that 
are not too dissimilar, results are known to be grossly in error whenthe 
cross-section is unusual, and such sections are customarily divided up into 


sub-areas and the hydraulic radius, average velocity, and discharge of each 
computed separately. 
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The velocity determinations along the perpendicular to the mid-point of t 
wall were made at closely-spaced points near the boundary, though it was 
realized that the pitot indications might not be reliable when the tube was he 
too close to the side. Four or more readings were made at each location, b 
as can be seen from Fig. 7, the variability is so great as to suggest that 
twenty or more measurements would be required to establish the average 
value accurately. Logarithmic and semi-logarithmic plots of the data were 
made, and empirical relations of either type could be fitted fairly well. The 
logarithmic type being simpler, it was used for fitting the plots for the diffe 
ent roughnesses, as shown in Fig. 8. In each case the fit is fairly good neai 
the wall, and diverges near the middle of the cross-section where the effect 
of the other boundaries becomes appreciable. 


Secondary Currents 


During the 1956 tests, it was noticed that floating debris thrown into the 
channel tended to become concentrated near the centerline as it moved dow! 
stream. This seemed to indicate the existence of secondary flow. Moeller 
Hartmann(15) gives a complicated pattern of secondary currents deduced fr 
the velocity distribution in a cross-section. No instruments were available 
that would permit the relatively weak secondary currents to be measured il 
such turbulent flow, but the relative strength of the current near the surfac 
was roughly measured by the following scheme. A large number of small 
pieces of wood, approximately 3/8 in. in diameter by 3/4 in, long were soal 
in water long enough so that they barely floated, Then, for runs with differ 
roughnesses and very nearly the same depth, particles were released from 
rest, one by one, at the shore and their distance from shoreline noted after 
they had travelled 16, 32, 48 and 64 ft downstream from the point of releas 
The results are shown in Table II and are plotted in Fig. 9. : 


Table II 7 : 
: 


Depth Number Average off-shore distance in fe 


Batten of of travelled while floating downstreé 
Date 5 % 
Spacing Flow Values : 
Averaged 16s 32 ft 48 ft 64 ft 
Aug. 19, 1957 6 in, 2.50 30-32 {OM IC SBee2e0sa 22s 
Aug. 23, 1957912 in; 2e5aetod=400 leek ly al R67 212 S07 ameo> 
Aug. 26, 1957 24 in, 20354 eer =28 O44. ehn02) a3 Poel sor 
Aug. 30, 1957 No battens 2.52 31 Bef We a ere eae CTs 


On the basis of this evidence, admittedly scanty, it would seem that fo 
turbulent flow in the triangular cross-section there is a large weak secon 
spiral that is directed inward at the top, from the shoreline to the cente rl 
(Although the floating-particle measurements were all made from the noi 
bank, there would presumably be another spiral, symmetrical about the 
centerline, as suggested by the 1956 observation.) There is a temptation 
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Vezocity IN FEET PER SECOND 


log-Loe Por oF VeLociTy 
VARIATION PERPENDICULAR TO 
Flume WALL aT 173 wiiD-Point 
(S,= 0.00/0/-0.00/03) 


conclude that the spiral is strongest for the 12 in, batten spacing, but more 
data should be gathered before such a conclusion can be regarded as estab= 
lished, | 


Backwater Curves 


-. 


Although the long variable-slope flume could well be used to investigate 
some of the cases of the backwater curves, the objective of the present stuc 
was the evaluation of losses in uniform flow, and the only data that might be 
useful to test backwater curve theory were those incidentally obtained while 
establishing uniform flow. After the practice of plotting the profiles at the 
time observations were made was started very little data of any possible ~ 
value for backwater comparisons was taken, but before that some interestil 
profiles were obtained which seemed to justify analysis as backwater curvé 


of Iowa. Fig. 10, reproduced from Mr. Singh’s thesis, shows a typical con 
parison of observed data and theoretical backwater curves. In general, Mi 
Singh found that the assumption ordinarily made in computing backwater 
curves, that the rate of friction loss is the same as it would be for uniforn 
flow at a given depth, gives very close agreement if velocity-head changes 
are correctly evaluated. He compared a number of recently-proposed bac 
water curve methods, and prepared a new table of functions that permits 


+) TRIANGULAR CHANNELS 59 


Batten Spacing in inches . 
24 12 


ged particles 
in feet 


-svubmer. 
line, 


totally 
shore 


Average distance from shore travelled 
by almos?- 
released of 


1.0 
A = Number of battens per foot 


Fig. 9. Floating-Particle Test of Strength of Secondary Current 


i computation of backwater curves in channels of triangular cross- 
on, 


CONCLUSIONS 


1is paper has reported on 156 runs on a long adjustable-slope 90° tri- 
lar flume. The length of the flume was found to be sufficient for 
letermination of normal depths up to approximately 2.5 ft. Including 
uns classed as critical flow along with those classed as shooting, 
‘esults may be summarized as follows: 


60 


Elevation in feet 


May, 1959 


Fig. 10 BACKWATER CURVES IN SMOOTH FLUME 
Test points compared with theoretical 
curves computed by K. P. Singh 


1956 Tests 
I 


Distance from Flume Entrance in feet 
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Item Smooth 
er of runs, tranquil 23 
er of runs, shooting 33 
»er of runs, total 56 
ing'sn, tranquil 0092 
ing's n, shooting 0094 
of se to fp» tranquil 0.93 
of f, to fp shooting LAS 


ness € in ft, tranquil 


ness € in ft, shooting 


Rough, batten spacing: - 


24 in, 
9 
16 


ZS 


. 0046 


0092 


127in.: 
9 
24 


33 


.0126 


0279 


6 in, 


16 


. 0343 


- 0524 


61 


Total 
Sin: 


18 (E) 
0 81 
18 156 


ese tests indicate that for this shape of channel the resistance is greater 
the flow is critical or shooting than when it is tranquil, For the case of 
uil flow with battens 3/16 in. high by 3/8 in. wide spaced A ft apart cen- 
center, the velocity was given with a mean discrepancy of 2.24% by the 


ula 


2 
“ - 0.145 (4.88 - 1 
cae [4.43 0.145 (4.88 - 1/) ] 
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where f, is the f for flow in a smooth pipe with D = 4 R and the same Rey- 
nolds number, as given by Prandtl’s formula. This also applies (but with a 
somewhat larger mean discrepancy) to flow in the smooth channel (by makir 
d infinite), The authors were unable to derive a similar formula for critica 
or shooting flow. 

A surprising result is that on the average these tests showed a smaller 
loss for tranquil flow in the smooth channel than for the corresponding smot 
pipe. Also contrary to the presuppositions of the authors, it was found that 
Manning’s formula gave about as good results as a logarithmic one, both fox 
the smooth and rough channels. 

Pitot-tube traverses indicated discharges that agreed fairly well with the 
weir measurements. The values of a varied from 1.03 to 1.09. 

Study of the velocity distribution for different roughnesses showed that 
the variation of velocity near the wall could be represented by exponential 
type formulas. 

The existence of two secondary currents with surface direction from sho 
to centerline was inferred from measurements of the motion of nearly- 
submerged floats, and the strength of these currents compared for the dif- 
ferent roughnesses. 

The ordinary assumptions made in computing backwater curves were fou 
to be justified, to the extent permitted by comparisons with the backwater 
curves incidentally obtained in getting the normal depths. 
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Notation 


a constant in Bazin’s formula; total area of cross-section in ft2, 


- a constant in Bazin’s formula. 


length of weir crest: a constant in Bazin’s formula. 
Chezy coefficient = V/VR S. 


coefficient of drag. 


- diameter of pipe in ft. 


ST oe eed 
1 


- partial area of cross-section. 
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Fig. 12. Rapid flow, battens at 12-inch centers. ete ~ second 
exposure shows ride-up on shore, Height of wetting not visible 
‘on far shore, but can be seen on near shore at bottom of picture. 
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f constant in Weisbach formula = 8 g R S/V2. 

fy f as computed from observations on flume. 

fp f for smooth pipe under comparable conditions. 

g acceleration of gravity in ft/sec2, 

h head on weir in ft. 

K coefficient in weir formula. 

n coefficient in Manning formula. 

Pp height of weir crest above floor. 

Q discharge in cfs. 

R hydraulic radius in ft. 

R Reynolds number for weir = g9-5 h 1.5 /y, 

IR Reynolds number for channel flow = 4 R V/v. 

r V2/Vp 

Ss slope of energy gradient. 

Vv mean velocity in cross-section, 

Va mean velocity of approach above weir crest. 

Vp mean velocity of approach below weir crest. 

v velocity in area dA. . 

WwW Weber number for weir = pgh2/o. 

Yn Depth for uniform flow. ; 

a velocity distribution factor = h v8 da j 
vA ; 

Y a constant in Bazin’s formula, 

€ equivalent Nikuradse roughness in ft. i 

r a constant in Gaukler’s formula; batten spacing in ft. 

v kinematic viscosity in ft2/ sec. 

p density in slugs/ft2. 

o surface tension in lb/ft. 

w i5/i. 


am Darcy, H, and Henry Emile Bazin, “Recherches Hydrauliques,” Mémoi 
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MODEL-PROTOTYPE STUDY OF A PLUMBING DRAINAGE SYSTEM®* 


D. A. Gyorog, ! F. M. Dawson2 M. ASCE and E. C. Lundquist? 


SUMMARY 


[his paper describes how dimensionless groups of variables may be used 
xpress the performance of a plumbing drainage system. It is also demon- 
ted that data obtained from tests on a scale model can be correlated with 
,¢ from the prototype in terms of the dimensionless parameters developed. 


SYNOPSIS 


"he results of the study described in this paper can be grouped into two 
gories, First, dimensionless analysis was used to reduce the number of 
umeters required to describe the performance of a plumbing drainage 
em. The data obtained from experiments run on test drainage systems 

3 put into the form of dimensionless parameters and plots were made to 
rmine the functional relationship existing between these parameters, 

nd, studies were made to determine whether the performance of a full- 
drainage system could be predicted from the performance curves plotted 
2 geometrically similar scale model. 

o geometrically similar drainage systems were constructed; one was a 
size drainage system; the other was a half-size model of the same sys- 
The performance curves obtained from the tests on these two systems 
rated that the performance of the full-size system could have been pre- 
2d on the basis of observations made on the model. However, it should be 
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pointed out that only steady flow conditions were considered in this study. 

A numerical illustration is presented at the end of this study which shows 
how the pressure at a typical point in a proposed system can be computed 
from data obtained from a scale model of the system. 


Analysis of the Problem 


A plumbing drainage system is composed of piping which carries water al 
entrained foreign matter away from plumbing fixtures in a building. All 
plumbing fixtures are equipped with traps which prevent odors and filth in th 
drainage system from flowing back into the open portion of the fixture. Most 
traps are only a few inches in depth, and pressures in the drainage system 
must not deviate positively or negatively from atmospheric pressure by 
amounts greater than the pressure head in the traps. If the pressure devia- 
tions are greater than this, the traps may be siphoned or blown out, thus des 
troying their effectiveness as barriers against transmission of odors. Vent 
pipes opening to the atmosphere at the roof of the building and connected to 
the drainage system are used to prevent large pressure deviations from oc- 
curring in drainage systems. If the vent pipes are of the proper size and ar 
connected to the drainage system at the proper points, the pressure througl 
out the drainage system will remain nearly atmospheric even when water ~ 
from the fixtures is being discharged. 

Laboratory experiments on full-size drainage systems are difficult and 
expensive to conduct because of their large size. It would be advantageous, 
therefore, if measurements made ona scale model could be used to predict 
flow and pressure conditions in a full-size system. 5 

In the research to be described two simple drainage systems were con- b 
structed; One a full-size system, or prototype, and the other a half-size ; 
geometrically similar scale model of the full-size system. Figure 1 shows 3 
the geometrical design of both the prototype and the model drainage a 


employed in this study. 

The prototype system employed a 4 inch diameter vent stack, and the 
model a 2 inch diameter stack. Transparent plastic pipe was used so charac 
teristics of the flow could be observed. For each experimental run, pres- 
sures at the various points indicated in Fig. 1 were measured by manometer 
connected to pressure taps on the vent stack. A constant pressure was main 
tained at A during a given test run through control of movement of the inve ti 
air tank. Observations were simultaneously made to obtain sufficient data t 
calculate the air and water velocities through the system. All measurement 
and observations were made only after steady flow conditions had been at-_ 
tained. j 

The drainage system was defined as that portion of the apparatus bounde 
by points A, B and C as shown in Fig. 1. Point A was designated as the en- 
trance to the vent stack of the system. In a real system this point would be 
open to the atmosphere, but for reasons to be explained later, the flow of ail 
into or out of the experimental system, and therefore the pressure at point 
was controlled by the inverted air tank shown in Fig. 1. The horizontal pipe 
through point B was used to admit water into the drainage system to simulat 
the discharge from a typical plumbing fixture. ; 

During the experiments the pressures existing at typical points in the s 7 
tem (such as D, E or F) for various air and water flow rates were studied 
closely. Under dynamically similar flow conditions in the model and protot 
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systems, it was found that measurements made on both systems could be co! 
related using dimensionless parameters obtained from a combination of the 
variables affecting the systems. 

The dependent variables in this problem are: 

P,-P, the pressure difference between the entrance to the stack and a 
point in the stack at which it is wished to determine the pressure (such as af 
D, E, or F in Fig. 1), lb/ft2 

Va, the air velocity in the plumbing stack, ft/sec 

The independent variables that must be considered are: 

Vy, the water velocity in the plumbing stack, ft/sec 


D, the diameter of the stack, ft 
g, the acceleration of gravity, ft/sec 
Ya» density of the air, lbs/ft3 


Y yw, density of the water, Ibs/ft3, and 


Pa-Pc, the pressure difference between the entrance and exit of the 
plumbing stack (at points A and C in Fig. 1), lb/ft2. 

It will be noted that two dependent and six independent variables are in- 
volved. The theory of dimensionless analysis was employed to combine thes 
eight variables into five dimensionless parameters. Furthermore, if no var. 
ables other than the eight defined above are pertinent in the definition of the 
problem, the physical relationship existing between variables can be ex- 
pressed as a function of the five independent dimensionless parameters that 
can be obtained employing this theory. In general, the theory of dimension- 
less analysis is useful for reducing the number of parameters required to 
express the function defining the physical phenomenon. The actual form of 
the function must be obtained from experiment or from other known facts as 
sociated with the phenomenon. % 

For this experiment Eq. (1) shows the five dimensionless parameters tha 
can be formed from the eight variables. The expression, fj, implies that a_ 
functional relationship exists between these five parameters. 


? 


a(PcP) [Va Vw Te 9a Pe) 
YaVia Vw YOO Yw wVw ; 
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Dynamically similar flow conditions will exist in the prototype and the - 
model if all the parameters in the above equation have identical values for 
the prototype and the model. Therefore, if the form of the function is dete 
mined experimentally with a model system, the results should be applicabl 
to the prototype. . 

The function expressed by Eq. (1) may be further simplified. If water ai 
air are the fluids to be employed in the model, as they are in the prototyp 
the ratio y ,/y,, will be virtually constant over the range of conditions li 


to be encountered in either the prototype or the model. This is true becaus 
the pressure and temperature variations that are apt to occur in either are 


not of sufficient magnitude to produce appreciable change in the density of 
the water or the air. = 
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g-P) g [Ma Vw 9(RR) 


ae ot A (2) 
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the viscosities of the two fluids had been admitted as independent varia- 
forms of the Reynolds’ number would have appeared in the functional 
onship. However, the effects of variations of the Reynolds’ numbers 
neglected in these experiments. Experience in other fluid flow problems 
utes that variations of pressure coefficients caused by small variations 
ynolds’ number are usually not great unless a critical region is en- 
ered. In this study it is argued that the viscosity effect and hence the 
lds’ number is of little importance unless the scale of the model is very 


sdure and Apparatus 


e€ model and prototype systems were arranged as shown in Fig. 1. Point 
; designated as the entrance to the plumbing stack. The pressure at this 
was controlled by the magnitude of the counterweight on the inverted air 
The rate of air flow into the vertical stack past point A was obtained by 
irement of the rate of change of displacement of the inverted air tank. 
tack exit at C was open to the atmosphere, and differences in pressure 
en the entrance and exit (P,-P() were created by controlling the rate of 
ment of the inverted air 
e pressure at the entrance to a plumbing vent is ordinarily atmospheric, 
.e pressure at the drainage exit may be something other than atmos- 
3, as determined by the pressure conditions existing in the sewer to 
‘the drainage system is connected, In this investigation the situation 
eversed for simplicity. The drainage end of the system (C in Fig. 1) 
en to the atmosphere, and the pressure was varied at the stack en- 
2 (A in Fig. 2). In this manner, the same effective pressure difference 
°c) was obtained. In addition, it was more convenient to measure the 
[ air flow through the system by measuring the rate of displacement of 
verted air tank placed over the pipe entrance as previously described. 
iter flow into the system was controlled by a valve upstream from B, 
e rate of water flow was measured with a weighing tank and a stop 
. Pressures at various points along the vertical stack were measured 
anometers, as shown in Fig. 1. 
i velocities Va and Vw were taken as the air and water flow rates 
tively in cubic feet per second divided by the cross sectional area of 
ck in square feet. These velocities, then, were only proportional to 
e velocities since the drainpipe was never running completely filled 


ither fluid at any one time. 


| (2) contains four variables, and it is necessary to employ two families 
s as shown in Figs. 2 and 3 to express the relationship between these 


les. In these figures the quantities g(P a-Po)/YwV% and g(Pa-P)/yaVa 
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SOLID POINTS FROM PROTOTYPE DATA 
OPEN POINTS FROM MODEL DATA 


Va /Vw : 
FIG. 2 EXPERIMENTAL RESULTS ~ 


: DRAINAGE SYSTEMS 73 


jotted against Va/Vw for constant values of V,,/VDg. . These curves 
plicable for only one location on the vent stack (point D, Fig. 1) at which 
ressure is P. Curves are not shown for other points along the stack, but 
were found to be of similar form. 

igs. 2 and 3 show that test points obtained from both the prototype and 
10del fall on common curves within the accuracy of the measurements. 
provides experimental verification of similitude in performance between 
rototype and the model. 

ig. 3 may be used to predict the pressure in a similar drainage system at 
nt corresponding to that for which Fig. 3 is drawn (point D in Fig. 1). In 
ion, the required venting air flow rate may be computed if the water flow 
the pressure difference between the stack entrance and exit, and the 
eter of the stack are known, This presumes, of course, that the system 
hich the calculation is being carried out is geometrically similar to that 
to obtain test data. 

1 example will serve to clarify the use of these data. Suppose a prototype 
bing system with a stack diameter of 6!'' and geometrically similar to 
ised in this study were to be set up. The problem could be that of deter- 
ig the pressure at point D (the point for which the curves of Fig. 3 are 

n) and the air velocity required in the stack for a given water rate of 

and for a given pressure differential between the entrance and exit of the 
m (Pa-Pc). If the pressure at the stack entrance, Pa, is atmospheric, 


ressure difference would be determined by those factors governing the 
ure Po. The flow conditions and discharge connections to the system 
stream from point C would govern the pressure Pc. 


i this example, suppose the water flow rate is such that 
Vw = 1.87 Ft/Sec. 

sat 
Py-R, = 0.433 Lb/Ft.? 
it will follow that 

 g(R-R) _ 32.2(0.433) 


: — = 0.064 
YwVy 62.4 (1.87) 


| 
| ace 
| VG2~ VO. StS2 2) 


shows that for these values 
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ily, the following values are obtained 


= 12.0 


Va = 2.0 Vw= 2(1.87)= 3.74 Ft./Sec. 


12.0 %4Va _ 12.0(.075)(3.74)? 
g 32.2 


= 0.392 Lb./Ft? 


CONCLUSIONS 


his investigation was carried out to determine the feasibility of employing 
snsionless groups of variables to reduce the number of parameters needed 
spress the performance of a plumbing system. The technique was applied 
simple system, and the performance curves obtained are not intended to 

[ value as specific design data. However, it has been demonstrated that 
from different size systems (so long as they are geometrically similar) 
oe reduced to such a form that common curves will define the performance 
ith systems. Specifically, it has been shown that four parameters are suf- 
nt to define such performance so long as the difference in the size of the 
2=ms is not so large that possible variations due to Reynolds’ number need 
onsidered. 


| 
| 
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RESISTANCE PROPERTIES OF SEDIMENT-LADEN STREAMS 


Vito A, Vanoni, 1 M. ASCE, and George N. Nomicos 


SYNOPSIS 


‘xperiments were made to determine the relative effects on the friction 
or of a flow of the suspended sediment load and the configurations which 
1 on sand beds of streams. Results showed that the friction factor of a 
am carrying suspended sediment is less than a comparable one without 
ment. They also showed that the reduction of friction factor due to 

ies in bed configuration is much larger than that due to suspended sedi- 
ts 


INTRODUCTION 


fhe resistance to flow of clear water in a channel with fixed boundaries 
been studied extensively and can be predicted with a satisfactory degree 
ertainty. Values of friction factors to be used in the flow equations such 
1e Manning or Darcy-Weisbach formulas have been determined and are 
in reference works on hydraulics as a function of the texture or rough- 
of the channel walls and the Reynolds number. 
‘owever, when a stream has a movable bed, and sediment is being trans- 
ed, the problem of determining the resistance is much more complicated 
‘in the simple case of clear water flowing in a channel with fixed walls. 
-rvations have shown that the friction factors for such sediment-bearing 
lluvial streams vary over appreciable ranges. However, the problem has 
een studied sufficiently to enable hydraulic engineers to predict a friction 
yr for given conditions. 


| 
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The fact that the friction factors for alluvial streams vary has been kno' 
for along time. Despite this, the problem is not discussed at any length in 
the literature, and textbooks are particularly silent on this point. This lea 
all but the expert in this field with the erroneous impression that the fricti 
factor of streams is constant and can be determined once for all time. The 
lack of discussion in the literature of the resistance of alluvial streams ma 
be explained by the fact that the problem is not well understood, and inform 
tion on it is often conflicting. For instance, some evidence presented show 
that the friction factor of a stream increases as the sediment load increase 
while other evidence supports the opposite conclusion. 

The variation of roughness of sediment-laden streams is caused by two 
distinct processes: (1) appearance of dunes and bars on the bed which may 
increase the roughness several fold because of the additional form resistan 
and (2) the damping effect of the suspended sediment on the turbulence in th 
stream tending to reduce the hydraulic roughness. Because of these two og 
posing processes, there has often been confusion about the effect of sedime 
load on roughness. Indeed, it is very difficult to separate the two effects | 
cause, for streams in which a significant amount of bed material is carried 
in suspension, the two processes are always acting together, but with differ 
relative effects. For example, in a sediment-laden stream flowing at high 
velocity over a flat bed of sand, the friction factor, or roughness coefficien’ 
will probably be lower than that for a clear flow over a fixed bed of compar 
ble grain roughness because of the damping effect of the suspended load. O 
the other hand, at a lower velocity with light suspended load, dunes may for 
and the net result will be an increase in the friction factor over that for a 
fixed flat bed of the same sand roughness. 

This paper presents results of experiments in which an effort was made 
separate the two processes in order to determine qualitatively their relativ 
order of magnitude. To accomplish this, runs were made with suspended 
and a movable bed, and then with clear water but with no change in the bed. 
configuration. This necessitated the development of a procedure for soli¢ i 
ing a natural bed configuration in place so that the clear water could be ma 
to flow over exactly the same bed as the sediment-laden flow. A comparist 
of friction factors then showed directly the damping effect of the sediment q 
separately from the effect of changing bed configuration. In the following 
pages a brief review of available information on the problem is given, and © 
then the experiments are described and the results and conclusions deriv aC 
from them are given. j 


Previous Studies of 
Resistance of Sediment-Laden Streams 


Early Studies of the Resistance Problem 


One serious obstacle to the understanding of the resistance of alluvial 
streams by the hydraulician of the latter part of the last century was the 
theory developed to explain the phenomenon, The reasoning back of the th 
started out with the logical idea that the loss in energy head between any t 
stations on a stream is merely the difference in elevation between the sta 
tions. It was then argued that this drop, or head loss, was used up in two 
ways: (1) in overcoming the hydraulic friction of the flow, and (2) in trans 
porting sediment. According to this concept, the more energy used in 
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pansporting sediment, the less was available to over come friction, and vice 
ersa. According to the Chezy formula, the velocity of a stream is propor- 
onal to the square root of the energy slope where the energy slope is merely 
1e loss of energy to overcome friction per unit length of stream. Now, if a 
irbid stream uses part of the total available energy to transport sediment, it 
ill have a smaller energy slope and hence a smaller velocity than it would 
ave if it carried no sediment. Thus, it was concluded that a clear stream 
ould have a larger energy slope and hence, according to the Chezy formula, 
ould flow faster than a comparable one carrying sediment. 

In view of present knowledge of flow problems, it is clear that some of the 
leas upon which the above conclusion was based were incorrect. To begin 
ith, much of the energy to transport sediment comes from the turbulence 
hich itself results from the expenditure of energy in overcoming the friction. 
herefore, this energy of turbulence is no longer available to overcome fric-~ 
on, and the fact that some of it is used to transport sediment should not nec- 
ssarily affect the dissipation of the main energy head. Also, apparently it 
as not recognized that the friction, or roughness factor, of sediment-bearing 
treams varies widely with load and flow rate, and that because of this the 
1ean velocity can vary even without a change in energy gradient or depth. 

The energy concept held by the hydraulicians of the last century, and the 
eductions resulting from it, were so firmly established that they influenced 
1e thinking of many workers of that era, and some vestiges of the ideas are 
) be found even in relatively recent literature. It appears from his early 
ublications that G. K. Gilbert(1,2) adhered firmly to this concept. This, no 
oubt, had an important influence on the thinking of his time because of the 
xcellent reputation he had amply earned in his active and productive career. 

Latham, in 1886, (3) reported that observations of a storm sewer four feet 
1 diameter in Croydon, England, over a period of years showed that the 
ean velocity for a given water depth was always greater when the water was 
lear than when it was turbid. He concluded that the retardation of the turbid 
ater was due to the load of sediment it was carrying. Hooker(2) p. 288, in 
o9mmenting on Latham’s conclusions, noted that Gilbert also reached the 
ame conclusion. Hooker also went on to state, “The law of conservation of 
1ergy will not admit any other decision in this matter, though the statement 
48 sometimes been made that such a retardation of velocity does not exist.” 
ne statement contradicting the belief that a stream at a given stage would 
ow faster when clear than when turbid, was made by McMath, and was in- 
uded in the paper by Hooker(2) footnote p. 289. McMath(4) was apparently 
ting some observations of the Mississippi River which contradicted the ac- 
apted views of the time. 

In his classical flume experiments in sediment transportation, Gilbert(5) 

229, expected to confirm the idea that a clear stream flows with less re- 
stance than a comparable loaded one. He selected pairs of runs from his 
‘periments which had the same depth and slope except that one of them car- 
ed sediment and the other did not. The bed in both cases was composed of 
avement of sand grains which did not move. According to the then preva- 
t theory, the mean velocity of the clear flow should exceed that of the 
uded flow. However, the opposite result was found in eleven of the sixteen 
rs of runs selected for the comparison. The depth measurements in these 

were apparently not as reliable as desired, and it was thought that this 
y have introduced a systematic error into the result. Therefore, Gilbert 
lected 13 more sets of comparable runs, 11 of which showed higher velocity 
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for the clear flows than for the loaded flows. However, it was also noted that 
the reduction in mean velocity due to sediment transport was greater for 
small load than for large, for gentle slopes than for steep, for low velocity 
than for high, and for large depths than for small. As a result of these data, 
Gilbert was forced to admit that the theory was not supported and that there 
was probably some physical law which escaped his analysis. 

A valuable contribution to the resistance problem was made by Buckley (6) 
in his work on the Nile. He made elaborate continuous observations of slope, 
mean velocity, and cross-sectional area at several stations on the river 
which showed an appreciable fluctuation of the coefficient of roughness. 
Typical values of these are shown in Table 1. 


TABLE 1. Measurements at the Khannag Nile Discharge 
Station by Buckley (6) 


River elev. Cond. of Flow rate Sed.conc. Mean 
Date meters flood meter /sec grams/meter Vel. 


(parts per million) meter/s¢é 


Aug. 16,1920 87.54 rising 7200 1900 1.48 


Sept. 16,1920 87.54 falling 6400 1000 1,32 


Data applying to this station on both dates: 


Cross-sectional area 4850 square meters 
Slope 0.000072 
Hydraulic radius 8.6 meters 


These data show clearly that at the given stage the mean velocity of the river 
is larger on the rising than on the falling stage, while other channel charac- 
teristics remain unchanged, The reduction in mean velocity represents an 
increase in the Manning friction factor of 11 per cent. A further pertinent 
point to be noted is that the sediment concentration, like the velocity, is 
higher on the rising than on the falling flood. 

Fig. 1 shows the measurements made by Buckley(6) in 1921 at Beleida 
Station on the Nile, which is about 28 miles upstream from Cairo and about — 
40 miles upstream from the Delta Barrage, or dam, which diverts water to — 
the canals of the delta system. The variation in slope shown is due to the 
change in level of the water at the Delta Barrage required to feed the canals, 
The value of Kutter’s roughness coefficient is seen to decrease appreciably 
as the discharge and mean velocity increase. For example, the following ex: 
treme values are obtained from Fig. 1: 


‘ 


On. ey a 


April 12, 1921 Sept. 8, 1921 
Discharge, cu. meters/sec 500 6200 


Kutter “n” -063 .027 
Hydraulic radius, meters 3.1 7.8 
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Such great reduction in channel roughness during the flood must have been due 
primarily to smoothing out of dunes and bars on the river bed, with the damp- 
ing effect of the suspended sediment being secondary, as will be seen fromthe 
laboratory experiments described below. 

In discussing Buckley’s paper (6) p. 231, Lacey pointed out that in the re- 
port of the Indus River Commission covering the period 1906 to 1910, it was 
stated that the velocity and discharge for a rising river were higher than for 
a falling river for the same stage, and the reverse was true for the water 
cross section of the stream, 

Kantlack(6) p. 265, cited evidence which disagreed with that obtained by 
Buckley. He called attention to the result of research on rivers by Swiss and 
French engineers which showed that the stream velocity was decreased in 
every case as the sediment concentration increased, According to Kantlack, 
this had been demonstrated prior to 1870. He also stated that in his experi- 
ence in the Punjab Canal system of India there was never an excess of flow 
during the time of high sediment load, as would be expected from Buckley’s 
findings on the Nile and as was reported by other workers in India. 

In the opening paragraph of his aber Buckley also called attention to 
some laboratory work by Duponchel which showed that a sediment load tended 
to increase the mean velocity of flow. 


Recent Studies of the Resistance Problem 


Interest in the resistance problem in the past two or three decades has 
been stimulated by both laboratory and field studies of the sediment problem, 
One of the main objectives of research in sedimentation is to develop a theory 
which will enable one to predict the sediment load of a stream as a function : 
of the water discharge and the stream channel characteristics. Once the dis- 
charge is given it is necessary to calculate the flow velocity and depth, whic 
requires that one know the roughness coefficient of the channel. Thus the 
need for knowledge of resistance of streams arises naturally in dealing with 
sedimentation problems. ; 

Experiments by Vanoni(7,8) in a flume with fixed, artificially roughened ~ 
bed showed that a flow with suspended load had a smaller friction factor than 
a clear water flow of the same depth in the same channel. These experiment 
were made first with clear water and then with various small amounts of 
sand, so that the sediment load was controlled. Experiments by Ismail(9) in 
a rectangular pipe showed a decrease in friction factor with increase in load 
in some cases, but the opposite trend in others, thus duplicating the conflict- 
ing results reported by field observers. The apparent paradox was probably 
due to the fact that in some runs the bed became roughened with small dunes 


making the roughness increase instead of decrease when sediment was adde 
to the system, 


ment, It is clear that the sediment-laden flow which has the higher velocity 
has the lower friction factor. The right-hand graph, in which the relative 
distance up from the bed is plotted on a logarithmic scale, shows that the 
velocity profile follows the von Karman logarithmic law. For two- 
dimensional flow, profiles can be expressed by 
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ahs 


d=0.295 ft d=0.295 ft 
S$=0.0025 = S=0.0025 


Oo 10 2.0 3.0 4.0 1.0 2,0 30 4.0 
Velocity U, feet per sec Velocity U, feet per sec 


Fig. 2. Linear and semi-logarithmic graphs of velocity profiles in 
a flow 0.295 ft deep and 33.5 in. wide with clear water and 
with a heavy suspended load of 0.1 mm sand. 


te a ee Bee ee 
k p k p 810d (1) 


sre u is the point velocity at a distance y up from the bed, u is the mean 
ocity at the profile, 7, is the shear stress at the bed, p is the mass den- 
y of the fluid, d is the flow depth, and k is the von Karman constant. From 


1 one can write, 
. ee ae T) ¥2 
uae Sipe merits, 18 (2) 


‘re uy and ug are the velocities at levels yy and yg, respectively. Since 
‘can be calculated from the measured slope and depth, one can calculate k 
m a pair of points on the velocity profile by Eq. 2. The effect of the sus- 
ded sediment in reducing k, as indicated by Fig. 2, is a general one which 
urs at all times(8) in the laboratory as well as in the field. 

ecent experiments by Brooks(10) demonstrated clearly how the friction 
or is affected by the dunes which develop on the sand bed of a stream. He 
ined friction factors for dune-covered beds which were as much as seven 
s as large as for flat sand beds of the same material. Studies by Ali and 
rtson,(11) and Barton and Lin(12) also showed clearly the importance of 
S on the channel roughness. 

scent studies of rivers have contributed appreciably to the understanding 
e resistance problems, Leopold and Maddock(13) presented detailed ob- 
ations of streams which show that at a given discharge high sediment 

S are associated with high mean velocities and low water depths. This 

ns that high loads are also accompanied by low friction factors, as was 
don the Indus River system by Lacey(6) over 35 years ago. 
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Eden(14) made a study of Manning’s roughness coefficient “n” as a functior 
of river stage for the lower Mississippi River. He found that n varied from 
about 0.040 at low stages to about 0.030 or less at flood stages. Carey and 
Keller(15) have recently investigated the bed configuration of the lower 
Mississippi River by means of a sonic fathometer and found that sand waves 
were present at all times and that the largest waves were formed at the high- 
est flow rates. They concluded that these waves must be a major factor in 
determining the friction coefficient of the river, although they were unable to 
measure it during their investigations. Carey and Keller are of the opinion 
that the loop in the rating curve of the River can be explained, at least par- 
tially, by a lag in the adjustment of the sand waves to the flow, i.e., during the 
rising flood the waves are growing and hence have a smaller resistance than 
they do at the same flow rate during the receding flood when the sand waves 
are diminishing, 

Harrison(16) presented data on the Missouri River near Fort Randall whic! 
showed that the Manning roughness coefficient diminished with increase in 
discharge. The minimum roughness coefficient obtained by him was about 
equal to that which would be expected by a sand roughness equal to the grain 
size of the bed material, so it seemed reasonable to conclude that dunes 
existed on the bed at lower discharges. 


Apparatus and Procedure 


The equipment and techniques used in the experiments are essentially 
those described by Brooks(10) and therefore are described only briefly here. 
The experiments were carried out by Nomicos and are described in his 
thesis(17) as well as in a report to the Corps of Engineers, U.S. Army(18) 
who sponsored the research. 


« 
Flume 4 


The flume used in these studies is 40 ft long, 10-1/2 inches wide, and 10 it 
deep. It is equipped with an axial-flow pump with a variable speed drive, 
venturi meter, transparent section in the return pipe, and electric immersion 
heaters to control the temperature of the water. The open flume section is — 
made of hot rolled steel, and its interior is painted with a smooth bitumastiél 
paint. Along the top of the channel are mounted steel rails and a movable 
instrument carriage which can measure vertical elevations relative to the 
flume to the nearest 0.001 foot. The entire flume is mounted on a truss, the 
slope of which can be easily adjusted even during the course of a run. Al- 
though the channel itself is 10 in. deep, experiments were made with depths 
of flow less than 4 inches to avoid excessive wall effects on the velocity and 
shear distribution of the flow. 


x 


Depth and Slope 


The elevation of the water surface and sand bed relative to the rails on 
flume was read with the point gage. To get the bed elevation once the water 
surface had been measured, the flow was stopped and the water drained out. 
The bed was then leveled in reaches of 3 or 4 ft with a scraper carried on tl 
instrument carriage rails. The elevation of this leveled surface was taken ; 
the mean elevation of the bed. The water depth at any point was taken as t 
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erence between water surface and bed elevations and the local elevation of 
energy grade line was calculated by adding the velocity head based on 

an depth to the water surface elevation. 

Profiles along the flume of the water surface, bed, mean depth, and energy 
de line were plotted for each run to determine if uniform flow obtained. 
ally several trial runs were required before uniform conditions were es- 
ished. The water surface and energy grade line were not always parallel 
he flume itself, even with uniform flow. Therefore, the actual slopes were 
uined by adding the slope of the flume and the slope of the water surface or 
rgy grade line relative to the flume. 


ocity 


The discharge was measured with a venturi meter with a throat diameter 
-00 in. placed in the return line, which was a 4-in. standard pipe. The 

in velocity was calculated as the discharge divided by the mean depth times 
flume width. Point velocities were measured with a 3/16 diam. Prandtl 
t-static tube and a water manometer reading to .001 ft. 


iment Concentration 


the mean sediment discharge concentration, C, was determined by 

pling the flow in the suction well of the pump at the downstream end of 
flume. The sampler was a brass tube 0.302 in. in inside diam. with a 
-loop at the end so that when the tube was held vertically during sampling 
inlet faced upward into the flow which was downward. 

*oint concentrations, c, in the flume were determined from samples taken 
a 3/16 in. outside diam. brass tube flattened at the end so the tip opening 
-040 in. high by .217 in. wide. The tube was bent in a right angle much 
the common pitot tube, and was mounted on the instrument carriage with 
tip pointing upstream. 

the velocity at the inlets to both samplers was kept the same as the local 
velocity by controlling the sampling rate. Each sample consisted of at 

t three one-liter bottles of the sediment-water mixture. The concentra- 

3 were determined by filtering, drying, and weighing the sediment collected 
ach bottle. 


perature 


he water temperature in all runs was maintained at 25° + 1°C. This was 
implished with immersion heaters in the return pipe which could be regu- 
i manually to give power input ranging up to 4,000 watts in convenient 


wall corrections 


nce the width-to-depth ratio of the flows in the experiments were small 
ared with natural streams, the effects of side-wall friction are relatively 
important. To reduce the laboratory results to correspond to those 

1 would obtain for much wider streams a correction is made. This fol- 

a method presented by Johnson(19) using a procedure developed by 
(20,18) This yielded an hydraulic radius rp, shear velocity U,,, and 


7-Weisbach friction factor fp applying to the bed itself relatively free of 
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side-wall influence. The friction factor fp is defined by the Darcy-Weisbach 
equation, 

Pe NST a at 
b ar, 2g (3) 


S = f. 


where § is the slope of the energy grade line, V is the mean flow velocity in 
the cross section, and g is the acceleration of gravity. To calculate the 
mean value, f, of the friction factor, one substitutes the hydraulic radius r 
of the section into Eq. (3) in place of rp. 


Sand Characteristics 


Cumulative size frequency curves for the two sands used in these experi- 
ments are shown in Fig. 3, and the size characteristics are summarized in 
Table 2, These sands were prepared from white foundry sands with well- 
rounded grains. The material was almost entirely silica, so that the mean 
specific gravity of the grains was 2.65. The mean settling velocity, w, of the 
various sands used was determined from sedimentation diameters of sieve 
fractions measured by Vanoni(21) where the fractions were obtained with a 
next of sieves with openings of adjacent sieves in the ratio 4/2. These 
sedimentation diameters are shown in Table 3. 


Procedure for Experiments , 

A run with uniform flow was first established by the usual techniques with 
determinations of depth, water discharge, sediment discharge, and slope. 
After the final determination of the bed profile by the leveling process, the 
flow was again started in order to regenerate the bed configuration which was 
destroyed when the bed profile was taken. When the bed was reestablished, 
the flow was stopped and the water carefully drained off. The bed was then 
solidified by spraying with chemicals, as explained in detail in the following — 
paragraph, After the chemicals had set, experiments with the stable bed were 
made, first using clear water and then with varying small amounts of sedi- 
ment added to the system. 


: 


TABLE 2, Summary of Characteristics of Sand Used in Experiments 


Dy é dete Lis wo 
Sand Experiments Geom.mean Geom.std. Mean sed. Mean fall 
No. used in, diam. deviation diam. velocity 
Set No. mm mm abe 
ee eS eh Se Se aie ee eae ey nce 
6 Ly APE 0.091 16 0.105 0.031 
7 IV 0.148 1.16 0.161 0.062 


Stabilization of sand bed 


The method of stabilizing the sand bed was developed by tests of sample 3 
in trays and in actual flume tests. The most satisfactory method involved 
the use of sodium aluminate and sodium silicate as stabilizing chemicals. 
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Fig. 3. Cumulative size frequency curves of the sands used in 
the experiments. 
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steps followed in the stabilizing process are as follows: 


HY 


a) Drain the water from the flume and allow the sand to dry for 24 to 48 
hours so that the moisture content is about 10 per cent by weight. 


b) Spray the sand with a mixture of the following three materials: 
1. 68% by volume of a 2%solution of sodium aluminate, 
2. 22% by volume of a solution of sodium silicate of specific gravity 1. 


3. 10% by volume of pure water. 


This mixture was applied to the bed with a paint sprayer until the sanc 
showed signs of saturation. This required an amount of the mixture 
about equal to 10% by weight of the sand being treated. 


c) After drying for about 12 hours, the sand was sprayed with a light ap- 
plication of a calcium chloride solution with specific gravity 1.2. 


d) Finally, a thin coat of synthetic varnish manufactured by Krylon, Inc., 


and sold under the trade name ‘Krylon Acrylic,” 


sand to form a waterproof surface. 


was sprayed on the 


One of the most difficult problems in the stabilization process was the 
preservation of the proper grain roughness of the surface of the bed. To dig 
turb the surfacé characteristics as little as possible, care was taken not ap- 
ply an excess of chemicals which would flood the surface, and the solutions 
were applied gently with a paint sprayer aimed parallel to the bed instead of 
at the bed, Although the results appeared satisfactory visually, there was, 
unfortunately, no method of testing whether the grain roughness of the finish 
surface was exactly the same as that of the loose sand itself. 


TABLE 3. Relation Between Sedimentation Diameter and Sieve 


Diameter Based on Vanoni (21) for Nevada White Sand 
D w. 


4 
3 
si i bet 7 
Tyler Sieve Sieve Mean Mean Fall Fall ; 
Sieve No. opening opening sieve sedimen- velocity velocity 
retained next retained diam. tation at 20°C at 25°C . 
on above on diam. 
mm mm mm. mm. fps fps ; 
. 
80 0.208 0-175 OL9R 0.198 0.078 0.083 j 
100 0.175 0.147 0.161 0.168 0.061 0.066 
115 0,147 0.124 0.135 0.153 0.053 0.058 
150 0.124 0.104 0.114 0.128 0.038 0.042 
170 0.104 0.088 0.096 0.107 0.029 0.032 
200 0,088 0.074 0.081 0.096 0.024 05.027. 
250 0.074 0.061 0.0675 0.086 0.020 0.022 


: STREAMS 89 


Results 


ne of Experiments 


he general plan of the experiments may be seen from Table 4, which 

'S some of the principal quantities measured. The experiments are di- 

j into sets and individual experiments. The flow depth, discharge, and 

1 velocity are kept the same in each set of experiments, but the slope is 
ad. In Sets II and IV the depth, discharge, and velocity are kept constant 
ighout, Sets I, Il, and III were made with a .091 mm sand and a .148 mm 
was used in Set IV. 

ach set contains runs with (i) movable bed, (ii) stabilized bed and clear 
r, and (iii) stabilized bed with varying amounts of loose sand added to the 
2m where the material added is the same as that in the bed. In Runs 2, 4, 
id 8, the bed was stabilized and no sediment was being transported by the 
i.e., the flow consisted only of clear water. 

able 4 shows the most important measured and calculated quantities ob- 
.d in the experiments. Where blank spaces occur in the table, the items 

. either not measured or not calculated. Measured values of the exponent 
the suspended load equation (Eq. 7) are listed in Fig. 8 for Runs £5395 
7, which were those with a movable bed. 


Configuration 

bservations of the bed configuration and the movement of the configuration 
made visually and by means of systematic photographs during the ex- 
ments. Fig. 4 shows a side view of the bed for Run 1 and a plan view for 
2, where the bed configuration for Run 2 is the same as that for Run 1 
pt that it has been stabilized by the method described above. Fig. 5 

§ a similar set of pictures for Runs 3 and 4. Fig. 6 shows the side view 
plan view of the bed in Run 5, which has been described as “flat” in 

e 4. It will be seen that the bed is actually not entirely flat but that 

2 are small ripples near the wall. The pictures of Fig. 4, 5, and 6 illus- 
» the meaning of the terms used to describe the bed condition in the 15th 


mn of Table 4. 


city profiles 


elocity profile measurements were made at the centerline of the flume 
ation 24, i.e., 24 ft downstream from the inlet to the flume. These 
surements fitted the logarithmic law (Eq. 1) very well, as may be seen 
.g. '7 which shows the profiles for Runs 1, 3, 5, and 6. The quantity m 
1 on the figure is the slope of the lines fitted to the points in fps per 

2, Eq. 2 can be transposed to give 

| 


u,- Uy Zee ie r 2.3 
mM = ——————_ = bet 
Y2 k § k 
Log === 
Soy 
which one gets 
53 
P ere 2. (4) 
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Fig. 4. Bed configuration, Set I, U = 1.23 fps, d = .284 ft. 
a) Run No, 1, side view, loose sand, during flow 


b) Run No, 2, plan view, stabilized bed, looking upstream, — 
without flow. 


HY 


Fig. 5. 


STREAMS 


Bed configuration, Set II, U = 2.02 fps, d = .244 ft. 

a) Run No. 3, side view, loose sand, during flow 

-b) Run No. 4, plan view, stabilized bed, looking upstream, 
without flow. 
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Run 
No. 


Notes: 


Key for 
(1) 
(2) 


(3) 
(4) 
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SUMMARY OF } 


Discharge 


1. Temperature for all runs: 25,0°C, 
2. ** Not measured 


Bed Condition: 

Dunes 

Dunes, largest at walls tapering down to nearly 
in center of flume. 

Ripples near walls, flat in center. 

Flat except for slight ripples extending about 2 


inches out from walls. : 7 


es 


_* 


STREAMS 


0.135 
0.150 
0.134 
0.135 
0.133 
0.133 
0.134 


for Bed Treatment: 

(a) Bed of loose sand covering bottom of flume. 

(>) Sand bed stabilized in natural configuration with no loose sand. 
(c) Sand bed stabilized in natural configurations with loose sand 

" 


added. l 
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Fig. 6. Bed configuration, Set Ill, U = 2.28 fps, d = .255 ft. 
a) Run No. 5, side view, loose sand, without flow 
b) Run No. 5, plan view, looking upstream, loose sand 
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Velocity, fps 


Fig. 7. Measured velocity profiles at the centerline of the flume 
at Station 24 for Runs 1, 3, 5, 6. 


vonKarman constant ke) for the centerline profile was determined from 
fy aorta 7 6) 


re the centerline shear velocity Uych is given by 
= af f a 6 
u SPAS —— u ( ) 


re Uc Q is the mean velocity calculated from the measured profile at the 
erline and f is the friction factor for the channel obtained from Eq. (3) 

g the hydraulic radius of the section r instead of rp. The values of Kf, 
+h are listed in Table 4, consistently decrease as the concentration C 
eases. 


ment Concentration Profiles 


rhe suspended load equation for distribution of concentration c of parti- 
with settling velocity w over a vertical is 
Z 
Se ey Beaks a) (7) 


a 


: c y 
e y is the distance up from the bottom, d is the flow depth, c, is the 
entration at the level y = a, and z the exponent is given by 

(8) 


Sie i Deve 
=] uP Koay 
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where 8 is the ratio of the diffusion coefficients for sediment and momentun 
Fig. 8 shows logarithmic graphs of (d- y)/y versus c for measurements mad 
on the longitudinal centerline of the flume. It is seen that these sediment pr¢ 
files follow Eq. 7, the well-established equation for distribution of suspended 
load. The values of the exponent z calculated by Eq. 8 using k,j from Tabl 
4, B= 1, and u, = uxc, (the centerline value defined by Eq. 6) are very 
close to those appearing in Fig. 8 which are determined from the slope of th 
lines in the figure. It is of interest to note that the profile for Run 1 agrees 
with the theory as well as the other runs in spite of the fact that it is the onh; 
one of the four for which the bed was covered with dunes. 


Friction Factor 


From Table 4 it is seen that the bed friction factor fp is smaller for 
sediment-laden flow over loose sand beds than for comparable clear flow 
over fixed beds of the same configuration. The fp - values for comparable 
sediment-laden and clear-water flows are summarized in Table 5. This 
shows clearly the effect of moving sediment in reducing the friction factor, — 
which in these experiments amounted to from 5 to 28 per cent. The results” 
shown in Tables 4 and 5 also indicate that the larger decreases in the frickiy 
factor f are associated with larger mean concentrations rea 

Attention is called to the fact that fp - values for the “identical” runs, 
Runs 8 and 8B, with clear water and stabilized bed, are not the same. Run 
was made once and then repeated on the following day, and the results agree 
closely. However, when the water had been in the flume for longer thana ~ 
total time of about 10 hours, some changes occurred, giving the results of L 
Run 8B. To further investigate this problem, Run 8B was repeated three 
times, once on each of three days, drying the sand bed after each experimel 
On the first day the water was run continuously. 


OE 
ee Beane 
Eera pa Pee eee: 1 


Runs |, 3&5: O.1Omm sand 
Run 7: O.16mm_ sand 


a2 
1.0, ia 1.0, ae 3 1.0, nt 5 10, Run 7 
c, Concentration, grams/liter 


Fig. 8. Measured concentration profiles on centerline of the flume 
at Station 24 for Runs 1, 3, 5, 7. 
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Qn the second day the experiment was completed in as short a time as 
sible and the pump stopped, leaving the water in the flume. After about 

se hours the run was repeated. The flow characteristics and friction fac- 
of the flow were changed by approximately the same amount as they would 
e been if the water had been running continuously. The run for the third 
was a repetition of that for the first, with approximately the same results. 
m this it was concluded that the water was in some way affecting the var- 
1 coat applied to the sand when the water remained in the flume for long 
lods. Run 8 is considered to give the valid results, with Run 8B being 
cted by swelling of the varnish between grains. 


Discussion of Results 


»ct of Sediment Load on the Friction Factor 


The effect of the sediment load on the friction factor may be seen in Table 
’ comparing the results of clear and sediment-laden water for each set of 
s. Sets I, Il, and III were made with the same sediment, so the results 
directly comparable. Set IV was made with a coarser material. It is 

ir that by adding a sediment load to a flow the friction factor is decreased, 
‘e€ in all four sets of experiments the sediment-laden flows had smaller 
tion factors than the comparable clear flows. The relative magnitude of 
reduction in friction factor seems to vary with concentration, the values 

| highest concentration producing the largest change. Referring again to 
le 5, it is seen that in Set I the friction factor fp for a flow with a con- 
ration of 3.64 grams per liter was 5 per cent lower than in the compara- 
clear flow. In Series II and III the reductions in fp were 25 per cent and 
er cent with concentrations of 4.60 and 8.08 grams per liter, respectively. 
m these figures it appears that rate of change in fp with concentration de- 
uses at high values of the concentration. 


TABLE 5. Comparison of Friction Factors for Sediment-Laden 
Flows with Those for Clear Water Flows over Stabilized 
Sand Beds of the Same Configuration. 


Sed, Disch, 


Concent, Friction Factors Decrease 
Run Depth Bed Cc Mean For Bed 
No, ft Cond. grms/liter £ i: Af o/o 
I 0.284 Dunes 3.64 0.074 0.106 
| . 006 5 
| 2 0.284 Dunes 0 COTTE 207112 
3 0,244 Small 4.60 0.0198 0.0211 


r dunes ,0072 25 
Del 2hGme0/ 0203.0... 


0.0165 0.0165 


.0064 28 
0.0203 0.0229 az 
0.0207 0.0227 

.0035 13.5 


0.0230 0.0262 
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The reduction in fp may be attributed to alterations of the turbulent struc 
ture of the flow by the suspended load. The mechanism of this process is 
discussed in more detail below. 


Effect of Bed Configuration on the Friction Factor 


The effect of dunes and ripples on the friction factor agrees with the data 
obtained by others.(10,11,12) For example, the large dunes in Set I cause fp 
to attain the high value of 0.106 in Run 1, compared with the values for all 
other runs which were less than 0.030 because of predominantly flat beds. 
From these data one can conclude that in flume experiments changes in fric- 
tion factor due to dunes can greatly exceed changes due directly to the movin 
sediment. 

In a study of the Mississippi River, Eden(14) observed that at Vicksburg 
the largest dunes occurred at the highest stage. At Memphis, he observed 
bed irregularities only on the rising part of the highest stage. The largest of 
these had heights much as 34 ft. Carey and Keller(15) observed bed configu= 
rations of the Mississippi River near Baton Rouge by means of a sonic fatho- 
meter. They also found that the largest dunes occurred at the highest stages 
There were several systems of dunes of varying wavelengths, the largest of - 
which had lengths of as much as two miles and heights of the order of 30 ft. © 
The shorter dunes were superposed on the larger ones. A study of these pro 
files showed that at low stage some of the longer waves disappeared but that | 
short waves of the order of 20 to 50 ft in length and 1 to 5 ft in height re- 
mained. Fig. 9, reproduced from the work of Carey and Keller, shows pro-— 
files at Donaldsonville gage for two stages. For the high stage which occurr 
on April 10, 1956, it is seen that there are waves of considerable length. Fo 
the lower stage which occurred in July, some of the longer waves are no 
longer present. Despite the increase in dune height with stage, Eden(14) 
found that the Manning coefficient n decreased as the stage increased. The 
sediment load also increases with stage, and can account for some reduc 9) 
in the roughness factor but not necessarily for all of it. It is possible that tk 
long. waves actually do not contribute much to the bed roughness and that the 
shorter, steeper ones, which seem to predominate at the low stages, may 
actually result in a higher roughness than the larger ones which have been 
observed at the highest stages. 

Data presented by Leopold and Maddock(13,22), an example of which is _ 
Fig. 10, for the spring flood of the Colorado River at the Grand Canyon (194 
show that for a given discharge the velocity is higher on the rising than ont 
falling flood means that the friction factor is lower on the rising flood, as- 
suming insignificant changes in slope and slow rate of change of stage. This 
follows the same pattern as observed in flume experiments where flows wit 
the higher velocities and sediment concentration tend to have lower friction 
factors, although it is not clear from the field data whether the damping ef: 
fect of suspended sediment or change in the bed configuration is the pre- — 
dominant factor in effecting the change. 

Leopold and Wolman(23) report that in degrading portions of the Rio 
Grande, the river channel characteristics follow a pattern similar to that 0 
the Colorado River shown in Fig. 10. In the lower reaches of the river whe 
the bed was aggrading, the velocity and sediment concentration were lower | 
the rising than on the falling flood, so that the friction factor was not highes 
on the rising flood. The lower friction factor is still associated with high — 
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Fig. 9. Bed profiles for lower Mississippi River (Louisiana) from 
Carey and Keller. (15) 


entration, but for some reason the sequence of change has been reversed. 
e same authors(23) also report that the characteristics of streams in the 
qumid areas, such as in the eastern United States, are much more stable 
those of the arid western part. In these streams the sediment loads are 

1 lower and the beds apparently much more stable, so that the changes in 
ed are probably much smaller than in streams whose beds are easily 

ad. 


t of Sediment Load on the Velocity Profile 


1e effect of suspended sediment on the velocity profile of a stream is il- 
ated by Fig. 2, which shows profiles for a clear flow and sediment-laden 
ith the same depth and slope. It is clear that in addition to having the 
velocity at any level, the sediment-laden flow also has a higher veloci- 
dient. Fig. 2 shows that very near the bed the velocity of the two flows 
ut the same, so that at some distance from the bed the flow with the 
velocity gradient will have the higher velocity. By differentiating the 
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» in the two flows the shear velocity uy is the same, the gradient for any 
} of y is inversely proportional to k, and since the gradient is seen to 
ase with load, k must decrease with the load. Thus it can be seen that 
ffect of sediment on the velocity gradient can be studied by observing k. 
1e fact that suspended sediment increases the velocity gradient was ob- 
ad by Gilbert(1) page 229, although the full significance of this was ap- 
itly not realized. Flume experiments(7,8) have shown that this is a 
matic effect, and observations on the Missouri River by the Corps of 
1eers(24) showed that this effect was also present in natural streams. 
mechanism by which the sediment increases the gradient, and hence the 

| velocity, has been explained in terms of the damping effect of suspended 
rial on the turbulence.(7) 

) keep sediment in suspension, i.e., to prevent it from settling due to 
tational force, work must be done on the sediment grains. This work 
come from the vertical components of turbulence fluctuations and must 
tt in damping of the turbulent motion. To investigate this problem, 

ein and Chien(25) correlated k with the ratio of the power, P,, to sus- 
sediment to the power, Pr, to overcome hydraulic resistance to the 
using data obtained from measurements on the Missouri River and from 
atory experiments by Ismail.(9) A similar graph was presented(8) for 
experiments made with 0.1 mm sand. Buckley(6) and Chien, (26) in dis- 
ng the effect of sediment on the flow, postulated that the main damping 

t of the sediment on turbulence occurred near the bed where the concen- 
In was highest. This appears reasonable since it is in this zone that the 
- rate is highest and most of the turbulence is produced. With this idea 
nd, k was correlated against the ratio of the power pi to suspend the 


nent in a thin layer near the bed to the power Py required to overcome 
rictional resistance. The power Pg to support sediment in a water 
n with unit horizontal area and depth d equal to the stream depth is 
Bice | ws 
Pp. = —2+———__ ¢t wd 
s es 
ey and Yg are the specific weights of the water and sediment respec- 
in c is the mean concentration over the depth in weight per unit volume, 
y is the settling velocity of the sediment. The power P¢ to overcome 
*iction of the prism of water is 


P= yaus 
2 § is the slope of the energy grade line and u is the mean velocity over 


pth. The ratio of the two powers is 


| Pp oa 

| a eiiea ae (12) 
ie! f Mie yus 

e power P - to support the sediment suspended between levels y = yj 

“=y2 ina prism of unit horizontal area is, 5 


P! (1 yy) Swly y,) 
= ——_—_—_— iS = 
8 Ye 1 ary 
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where cy is the mean sediment concentration in weight per unit volume be- 
! : 
tween the levels y; and yg. The ratio of P, to Pf is 
es Cc, WwW > Pete 
ES ES ae ) | a eS eee (13 


P, ¥s yus d 


Fig. 11a shows a graph P'./ Ps plotted against k for flume data obtained 
with uniform sand of 0.1 mm size, where y, and yg have been taken as .001 
and .01 times the flow depth respectively. Fig. 11b shows a graph of the 
same data(8) plotted according to the method of Einstein and Chien(25) i.e., 

k is plotted against P,/ Pr. It is seen that the correlation of the data in Fig. 
11a is better than in 11b. This tends to support the ideas of Buckley and 

Chien and may also explain why the velocity remains logarithmic in a 
sediment-laden flow, as seen in Fig. 2, even though the turbulence has been 
modified considerably. If significant damping of the turbulence occurred . 
over the entire depth, one would expect a distortion in the distribution of the 
turbulence, and hence in the shape of the velocity profile. However, if the t 
principal interference with the turbulence occurs near the bed, the turbulence 
energy will be reduced but it may still follow the same laws in diffusing 
through the section, and one might expect the same shape of velocity aan 


tion. 

Referring again to Fig. 2, it is seen that the two flows have the same dept 
and hence, hydraulic radius, and the same slope. The one with the higher 
mean velocity, i.e., the sediment-laden flow, will have the lower friction fac- 
tor. The decrease in f, shown in Table 4, for the sediment-laden flows, re 
sults from the damping phenomenon discussed above. A comparison of the 
k-values listed in Table 3 will show that the turbid flows have much the lowe! 
values, thus indicating that damping has occurred. 


General Discussion 


The present experiments have shown that a sediment-laden stream acts il 
two ways to change its roughness coefficient: (1) a smoothing-out of the bed 
configurations tends to reduce the coefficient as the velocity increases; dam 
ing of the turbulence by suspended sediment tends to increase the velocity — 
gradient so that the difference between the velocity near the bed and near thé 
surface increases along with the mean velocity. In the present laboratory 
experiments with 0.091 mm sand, the change in the bed friction factor fp re 
sulting from change in bed configuration alone was as much as five-fold in 
going from a dune-covered bed to a flat bed. On the other hand, the maximt 
reduction in fp due to the sediment load was only 28 per cent (Run 5A). Sint 
the Manning roughness coefficient n is approximately proportional to Vip? 


the corresponding changes in n are about twofold and 13 per cent, respecti 
ly. It is significant that the run which showed the maximum reduction in 
friction factor, i.e., Run 5A, also had the lowest value of k as determined 
from the centerline velocity profile. From these few data, it appears that 
change in bed configuration has a much larger influence on the roughness ¢ 
efficient than the damping effect. From this one would judge that the large 
reductions in roughness coefficients reported by Buckley for the Nile, 
for the Mississippi, and Lindley for the Indus River, could not be the res 
only of damping, but must have been also affected by significant changes 
dunes on the bed. Judging from the flume experiments one would expect th 
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‘Rig. 11. Reduction of the von Karman constant in sediment-laden 
flow. 


to diminish as the velocity increased, while, at the same time the load 
Ls Since both of these effects tend to reduce the roughness coeffi- 
the field observations are in agreement with flume studies and the 

s can be explained at least qualitatively by the combination of the two 


Perr iion of the observations that for the same discharge the rough- 
Oefficient of a stream on a rising flood is smaller than on a falling 
oresents some difficulties. Buckley and Lindley noted that for a given 
e concentration was higher during the rising flood than during the fall- 
od, and explained the reduction in roughness by the effect of the higher 
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load. Judging from the results of the present experiments as shown by 
Table 5, the changes in roughness factor observed in the field are too large 
to be accounted for by damping effect alone, and some readjustment of the be 
must be involved in order to achieve the changes observed. During the risin 
stage the stream bed probably becomes smoother because the transportation 
of a-high sediment load requires a lower depth and higher velocity compared 
with the falling-stage flow for which the transport rate is usually less; how- 
ever, since the damping effect also reduces the roughness, it is not possible 
to determine the relative magnitude of the two effects from the field data 
alone. Bottom configuration data of the kind obtained by Carey and Keller, 
along with stage, slope, discharge and load data, would do much to clarify 
this problem. 

The field data of Buckley, Lindley, and Leopold and Maddock on the varia: 
tion of stream characteristics, report only total sediment concentration and 
do not give size distribution of the load or division between wash load and be 
material load. As seen from Fig. 11, the effect of concentration on k, and 
hence upon friction factor, depends on the settling velocity w as well as the 
concentration, so that for the same C, coarse material with high settling 
velocity has a greater effect on k and on the reduction in friction factor. By 
the same token, fine materials such as silt and clay have low settling veloci- 
ties, and hence small damping effects even for large concentrations. The 
loads in streams such as the Indus, Colorado and Rio Grande, which show a 
correlation between mean concentration and velocity (or friction factor) (6, 13 
are predominantly fine sediments which should have small damping effects, — 
yet changes in load are accompanied by substantial changes in friction factor 
for the same discharge. According to laboratory results it is unlikely that a 
of this is caused by damping effects and changes in bed configuration are 
probably involved. ‘ 

The discussion thus far has dealt with those observations that show a re- 
duction in friction coefficient as the load increases; however, it will be re- 
called that Kantlack(6) reported that Swiss engineers found that the velocity 
was reduced by an increase in the load, and that Latham(3) reported a simi . 
lar behavior in a storm sewer, The latter might be explained by postulating 
that dunes were formed when sediment was being transported, thus causing 
high friction factor, and that during clear flows sediment was cleaned out of 
the sewer, leaving the water in contact with the smooth walls, which resulte 
in low friction at low flow. The explanation of the behavior of the Swiss 
rivers may be associated with a phenomenon observed by Brooks(27) ina _ 
flume at low transportation rates. He observed that there was a tendency a 
low velocities (i.e., low transportation rates) for the friction factor to in- 
crease with the velocity to a maximum, and then to decrease as the veloci ty 
increased further. The experiments which showed these results were mad 
in a flume 33.5 in. wide. Two similar sets of experiments made by Brooks 
in the 10.5-in. flume used in the present experiments did not yield a maxi- 
mum for fp but showed f, decreasing as the velocity increased for all valu 
of the velocity. Therefore, the explanation offered above, based on this on 
set of experiments, must be considered as speculative. 


CONCLUSIONS - 


Four sets of runs were made in a flume 10.5 inches wide in an effort to 
separate the direct effect of suspended sediment on the friction factor fror 
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f the bed configuration. First, a natural run was established over a 
bed of fine sand, and then, after draining the water, the bed was chemi- 
solidified in its naturally formed configuration. Loose sand was re- 

d from the system then and a flow of clear water was run over the 

ized bed, with the usual measurements being taken. Following that, 
increments of sand were added, noting the change in the friction factor 
resulted. From these experiments the principal conclusions are as 

1S: 


Suspended sediment reduces the friction factor of a flow. In a limited 
series of laboratory experiments, the amount of this reduction varied 
from 5 to 28 per cent. 


The changes in the friction factor due to variation in bed configuration 
are much larger than those due to the suspended load itself. A friction 
factor may change by a factor of 3 to 6 in the laboratory flume between 
runs of differing velocity at constant depths. 


Evidence indicates that the suspended load reduces the friction factor 
by damping the turbulence or interfering with its production near the 
bed where the sediment concentration is the highest and the rate of 
turbulence production is the greatest. 


As observed by previous investigators, the von Karman constant k is 
also substantially reduced by the suspended sediment. 


In the field it is impossible to distinguish between the two above- 
mentioned factors affecting alluvial channel roughness. However, it 
has been established that the reduction due to the suspended sediment 
is of importance only for streams carrying a very high suspended load 
over a flat bed, and is of minor importance when there are dunes on the 
bed. 
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GRAVEL BLANKET REQUIRED TO PREVENT WAVE EROSION 


Enos J. Carlson, M. ASCE 


SYNOPSIS 


rface wave erosion tests were performed in a hydraulic wave flume on 
rials shipped from Yakima Project, Washington. The tests were made to 
mine the cover blanket needed to prevent leaching of fine base material. 
lar talus and rounded gravel materials were shipped to the laboratory 
the field, together with fine base material to make the tests. 
e most stable cover blanket was formed by screening the talus material 
3/4-inch screen. The minus 3/4-inch material being placed over the 
material with the plus 3/4-inch material on the surface. Motion pic- 
were taken for a record of the 14 tests. 


INTRODUCTION 


me factors which must be considered in the design of unlined and earth 
canals are those which may cause erosion and, consequently, develop 

in extensive maintenance problem. They include: excessive tractive 

s which are caused by high velocities and, consequently, result in ero- 
over a large part of the wetted perimeter; and those forces caused by 
ce waves which produce erosion on the banks at the water surface. The 
t and rate of erosion produced by these forces which exist in many un- 
and earth lined canals designed and constructed by the Bureau of Recla- 
m, depend upon the type and gradation of the soil material through which 
anals are constructed. 

1e Kinnewick Main Canal, feature of the Bureau of Reclamation’s Yakima 
ct in the state of Washington, was constructed in an area where the ma- 
_is fine and silty and has little cohesion. It is described as a highly 
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compressible micaceous, nonplastic silt which had its origin in old lake bedi 
When dry, it is easily crumbled into a powdery dust. Over 70 percent of the 
silty material passes the No. 200 United States standard screen. This type ¢ 
material, combined with the prevalent wind in the area, provided the condi- 
tions for possible erosion on the canal banks due to waves created by the wil 
The prevailing direction of the wind is generally west to northwest and is in 
the same general direction as a large part of the canal. The many curves 
which were necessary on the canal, constructed in hilly topography, are es- 
pecially vulnerable to erosion caused by surface waves. It was decided that 
some type of protective cover over the fine material for the Kennewick Mair 
Canal was essential, particularly on the curves. A hydraulic model study w 
made to help determine the adequacy of the available materials to be used a 
a cover blanket, and determine the method of placement that would give the 
best protection for the least amount of material and cost. 

The Kennewick Main Canal is located in Benton County, Washington, near 
the Yakima River, between the towns of Prosser and Kennewick (Figure 1). 
The canal is approximately 42.3 miles in length and is designed to carry a 
discharge of from 500 to 275 cfs. The unlined and earth lined sections have 
a 14-foot bottom width and side slopes of 2:1. The depth varies from 7.61 
feet to 5.66 feet, corresponding to the discharges cited above. Construction 
of the canal provides irrigation for approximately 14,500 irrigable acres of 
new land and additional irrigation for approximately 4,600 acres of land in 
the area that had been under irrigation. The design capacity of the canal wi 
permit eventual irrigation of approximately 6,100 acres of additional land. 


Hydraulic Model 


As shown in Figure 2, a 2:1 slope of plywood was built in an existing flun 
which was 2 feet wide, 6 feet deep, and approximately 24 feet in length. To 
make a maximum number of tests with the material which had been shippe 
from the field to the Bureau’s Hydraulic Laboratory in Denver, a plywood — 
divider was placed on the slope, dividing it into two sections. The 1-foot- é 
wide area on the side of the plywood divider next to the glass observation © 
window was used for placing the test section. On the other side of the divic¢ 
a gravel wave absorber was placed to reduce reflections from this area. — 

A wave machine, Figure 3, consisted of a bulkhead hinged at the bottom 
the flume and moved back and forth about the hinge by an arm connected to 
large pulley wheel. The pulley wheel had holes drilled at various radii fro 
the center, thus giving a stroke to the wave-making bulkhead which could be 
varied from 2-1/8 to 11 inches. A 3/4-hp motor connected to a set of vari 
ble speed pulleys completed the wave-making apparatus. With this setup, i 
was possible to change the frequency of waves from 6 to 315 waves per 
minute. 

In these model tests, the waves were made to move at right angles to tht 
canal slope, and resulting beaching action was more severe than the same 
height of wave would have been in a prototype canal where the waves would 
move longitudinally along the side slopes. However, it was possible to con 
pare the various test setups with this condition and thereby select the best 
test arrangement. With the various setups on the wave-producing machine 
it was possible to obtain wave heights up to 0.75 foot from trough to crest. 
Prior to operation of the wave machine, the model was filled slowly to the 
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desired depth which was maintained constant during a test. Water surface 
elevations were determined by using a hook gage in a well outside the wave 
area of the model. 

For each test, both still and motion pictures were taken as an aid in re- 
cording the data from the model. The test procedures were designed to de- 
termine which cover material available on the project would give the best 
protection for the fine silty base material and to determine which would be 
the most efficient way to place the material to give the best resistance to 
erosion and least leaching of the base material and still at the minimum cost. 
In the early tests, a capacitance-type water surface indicator, together with 
a recorder, was used to determine the type of wave which should be produced 
as a standard for all tests. After recording the various waves that were pro- 
duced by adjusting the wave machine to various strokes and frequencies, 
standard wave patterns were set, and the recording apparatus was no longer 
used. 


Test Slope 


Filter Criteria 


On the project there were two materials available which could be used as 
a cover blanket for the Kennewick Main Canal. Both were’equally available 
and in adequate supply for the canal job. Both of these materials fell largely 
outside the Bureau of Reclamation filter criteria for the base material, 
as shown on Figures 4 and 5. The placing of a cover blanket over finer base 
material on a canal bank is similar to the filter problem for the foundation of 
engineering structures to allow drainage and prevent leaching of the base ma: 
terial. The following criteria are used by the Bureau of Reclamation for fil- 
ters in foundations work. F 


5 
’ 
5 


a. The filter material should be more pervious than the base material to 
assure there will be no hydraulic pressure built up to disrupt the filter and — 
adjacent structures. 


b. The voids of the in-place filter material must be small enough to pre- 
vent base material particles from penetrating the filter which causes cloggi 
and failure of the protective filter system. 


c. The layer of the protective filter must be sufficiently thick to provide 
good distribution of all particle sizes throughout the filter and also to provi 
adequate insulation for the base material where frost action is involved. 


Results of a large series of tests resulted in the following requirements | 
for a graded filter. A graded filter material is defined in this paper as a 
material having a comparatively broad range of particle sizes. Graded ma- 
terial may have concave, convex, S-shaped or straight line gradation curve 
and may be defined as “poorly” and “well-graded” material, depending on 
their gradation curve shape. Uniform grain-size filters have a similar set 
of requirements. 

For graded filters, it is necessary that the ratio of the 50 percent grain 


size of the filter material to the 50 percent grain size of the base mater 
be between the limits of 12 and 58 
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_ 50 percent size F.M. 


R' = 
50 percent size B.M. 


= 12 to 58 
ratio of the 15 percent grain size of filter material to the 15 percent 
n size of the base material should be between the limits of 12 to 40 


RM" = 15 percent size F.M. 


~ 15 percent size B.M. | ie 


n addition, the following requirements for a graded filter should be met: 


. The filter material should pass the 3-inch screen for minimizing parti- 
segregation and bridging during placement. Also, filters must not have 

e than 5 percent minus No. 200 particles to prevent excessive movement 
nes in the filter and into drainage pipes causing clogging. 


. The gradation curves of the filter and the base material should be ap- 
imately parallel in the range of finer sizes, because the stability and 
yer function of protective filters depend upon skewness of the gradation 
ye of the filter toward the fines, giving a support to the fines in the base 
erial. 


. The filter material adjacent to the drainage pipe should have sufficient 
"se Sizes to prevent movement of filter material into the drainage pipe. 
maximum size of the perforations or joint openings of the drainage pipe 
selected as 1/2 of the 85 percent grain size of the filter material. 


. In designing filters for base materials containing particles larger than 
4 size, the base material should be analyzed on the basis of the gradation 
aterial smaller than No. 4 size. 


serials Tested 


he gravel, Figure 4, was rounded river-deposited material. Twenty 

ent of this material was smaller than the No. 30 United States standard 
e, and 80 percent was smaller than the 1-1/2-inch size. Between the 

30 and the No. 4 sieves, there was only 4 percent of the material. Be- 
e of the dearth of material in this range, the gradation fell largely out- 
the limits of the graded filter criteria for the base material. The other 
rial available for a blanket cover was an angular talus found at the base 
cal mountain slopes. This material was obtained from two locations. 
combined material gives a gradation of 10 percent smaller than the No. 
mited States standard sieve, and 80 percent smaller than the 1-1/2-inch 
ing. The talus material was better graded than the rounded gravel ma- 
i; however, it also fell outside the limits for a graded filter criteria for 
vase material, as shown on Figure 5. 
gradation of the base material, Figure 6, showed 6 percent finer than 
mm, 27 percent finer than 0.037 mm, 74 percent finer than 0.074 mm, 
small amount of material as large as 3 to 5 inches. Several tests to 
rmine the liquid limit and plastic limit showed the base material to be 
astic. 


Hydraulic Wave Tests 


1—Gravel Moistened and Well Mixed ‘ 
e first test was conducted on the gravel material, Figure 7. Six inches 
‘material moistened and well mixed was placed uncompacted on the 2:1 
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WAVE EROSION 


)) Slope failure after 45 minutes 
of operation at 32 waves per 
minute - wave height 0. 32 ft. 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
y TEST 1 (GRAVEL MOISTENED AND WELL MIXED) 
, ; Model study 


Figure 7 
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slope. The moist sand particles adhered to the gravel and were well dis- 
tributed throughout the 6-inch layer. The depth of water was raised to 2.7 
feet, and the wave machine was operated at 32 waves per minute. This pro- 
duced a wave having a height from trough to crest of 0.32 foot. Immediately 
after the wave machine was started, the gravel and sand material began to 
slip and roll down the slope. The sand which had adhered to the gravel par- 
ticles began washing from between the particles, reducing the friction throug 
out the layer, and the force of the water flowing out of the blanket when each 
wave receded from the slope was enough to start the particles to slide and 
roll. Approximately 20 minutes after the beginning of the test, the board sloj 
was exposed. The photograph in Figure 7b shows Test 1 after the test had 
been in operation for 45 minutes. Complete failure near the water surface 
had occurred. It seems possible that if the water surface were raised and 
lowered a few times, allowing the sand particles to wash from between the 
gravel particles, that a greater coefficient of friction could be obtained 
throughout the gravel. This cover layer then might have been more stable. 
The fine sand particles acted as ball bearings, allowing the larger gravel 
particles to roll and slide, causing failure of the cover blanket. 


Test 2—Gravel Separated on 3/4-inch Screen 


The gravelly material used in Test 1 was dried and separated on a 3/4- 
inch screen. The resulting analysis curve and Bureau of Reclamation filter 
criteria are shown in Figure 8. The analysis curve for the fine base materi- 
al, minus 3/4-inch material, and plus 3/4-inch material are shown. Filter — 
criteria limits for the base material and the minus 3/4-inch material are 
shown on the same graph. The graph shows that the minus 3/4-inch material 
falls roughly within the limits of the graded filter criteria for the base ma-_ 
terial, except that there is very little material between the No. 4 and No. 30 ; 
sieve sizes. The 3/4-inch material falls within the limits of the uniform fil- 
ter criteria for the minus 3/4-inch material. 

A 4-inch layer of the material smaller than 3/4-inch size was placed on — 
the 2:1 slope and covered by a 4-inch layer of material larger than 3/4-inch 
in size. The water was raised to a depth of 2.7 feet, and a wave frequency — 
was set at 32 waves per minute. The resulting wave height was 0.32 foot. TI 
model, ready for operation, may be seen in Figure 9a. After 1 hour’s opera- 
tion at these conditions, there was practically no observable damage on the 7 
cover blanket. The wave frequency was then increased to 42 per minute, wit 
the water depth maintained the same. This increase in frequency caused a 
wave height to increase from 0.32 to 0.43 foot. a 

After operating the model for 23 hours at 42 waves per minute, there was 
only a small amount of beaching. It was believed the resulting damage, shoy 
on Figure 9b, was not sufficient to cause failure in a canal installation. The! 
was some movement of sand in the gravel below the 3/4-inch layer. The sa: 
appeared to move down the slope through the gravel in the area where the ty 
layers were in contact. As may be seen by studying the photographs (Figure 
9), the sand moved in the top 1-1/2 inches of the 4-inch layer of minus 3/4- 
inch material. Some of the sand in the material smaller than 3/4-inch in si 
concentrated at the bottom of the layer, thus, probably improving the action 
material and preventing beaching. ‘ 
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Figure 9 


(b) After 23 hours of operation at 42 


waves per minute--wave height 
0. 43 ft. 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 2 (GRAVEL SEPARATED ON 3/4-INCH SCREEN) 
Model study 
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3—Gravel Placed Dry and Raked 


-attempt was made to find a less expensive method of placing the cover 
‘ial and stabilizing it with a raking or harrowing process. It was thought 
ith harrowing the gravel material, the finer material would tend to set- 
the lower portion of the cover and the coarser material migrate toward 
irface; in other words, a separation with the coarse material on the sur- 
was attempted to be produced by raking, Figure 10, rather than screen- 
sS was done in Test 2. 

© wave machine was started and operated at 32 waves per minute, with 
e height of 0.32 foot for a short time. This created only a small amount 
vement of the finer material. The wave machine was then stepped up to 
ves per minute, with a resulting wave height of 0.43 foot. The depth of 
was again held constant at 2.7 feet. As can be seen on Figure 10b, the 
wed material withstood the wave action much better than the unraked 
10ist material placed in Test 1. However, beaching did start near the 
-surface as the machine was first turned on, and a small amount of 

ing continued until the wave machine was turned off 24 hours later. 
photograph showing the model after it had been in operation for 19 hours 
ye seen in Figure 10b. The final beach slope, after the model had been in 
tion for 24 hours, was approximately 5.14:1. This slope seemed to be 

r stable at the end of the 24-hour period. 


4—Gravel Separated on i-inch Screen 


Test 4, the gravel was separated on a 1-inch screen, and 4 inches of 
‘ial smaller than 1-inch in size was placed on the 2:1 slope and covered 
| inches of the plus 1-inch material. The resulting analysis curve 
d with the base material, and filter criteria are shown in Figure 11. 
nodel was filled to a depth of 2.7 feet, turned on, and operated at 42 
3 per minute for approximately 45 minutes. Although the surface ma- 
‘withstood the 0.43-foot waves very well, some sand moved from the 
layer quite freely. 
sure 12a shows a photograph of the model after 45 minutes of operation. 
rea from which the sand moved in the lower layer is seen in the photo- 
. At this condition, the surface material was very stable. To create a 
severe erosion condition, the water surface was then raised to a depth 
feet. The wave machine was operated at 42 waves per minute, and the 
ting waves increased from 0.43 foot to 0.5 foot from trough to crest. 
this larger wave, there was some movement of the surface material; 
rer, it was not considered serious. Increased movement of the sand 
the material smaller than 1-inch in size was noted. 
study of the photograph in Figure 12b reveals that the sand had moved in 
wer layer to the full depth of 4 inches. It was obvious that the compara- 
| larger voids in the plus 1-inch material allowed a flow of water through 
avel great enough to cause the fine sand to leach out through these voids. 
ger voids would not reduce the force of the water enough by the time it 
the fine sand material in the lower layer. The fine base material 
under a protective blanket similar to that which was placed for Test 4 
apparently leach out, as shown by the test. 
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Figure 10 


(a) Rake used in Test 3 


ee ee 


(b) After 19 hours of operation at 42 
waves per minute--wave height 
0. 43 ft. 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 3 (GRAVEL PLACED DRY AND RAKED) 
Model study 
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Figure 12 
(a) After 45 minutes of operation-- 

2.7 ft. depth--42 waves per 
minute--wave height 0. 43 ft. 
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(b) After 20 more minutes of 4 


operation--3. 2 ft. depth-- 
42 waves per minute-- 
wave height 0.5 ft. 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 4 (GRAVEL SEPARATED ON 1-INCH SCREEN) 
Model study 
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sis 5 and 6—Gravel Separated on 1/2-inch Screen 


To find the optimum screening size, tests were continued, and in this test 
gravel was separated on a 1/2-inch screen. For Test 5, the resulting 
lysis curve with filter criteria and the bed material curve are shown in 
ure 13. Four inches of material smaller than 1/2-inch in size was placed 
the 2:1 slope and covered with 4 inches of material larger than 1/2-inch in 
e, Figure 14a. The material in this test was screened while damp, and it 

; noticed that the sand adhered to the coarse gravel particles even after 
siderable shaking in the screening machine. This same condition occurred 
Test 1. Water was again raised to a depth of 2.7 feet, and the wave machine 
; operated at 42 waves per minute, creating waves with a height of 0.43 foot. 
er a few waves had washed up the slope, a slide occurred. This slide ap- 
red to start in the area where the water from the waves flowed out of the 
vel material and continued to the toe of the slope. A photograph showing 
condition of the slope after the slide may be seen in Figure 14b. 

Motion pictures taken of this test, including the slide, are included in the 
ord of the hydraulic model tests. From viewing the motion pictures, it 
eared that the sliding surface started at the surface of the plus 1/2-inch 
terial, progressed through the plus 1/2-inch material, then through the 

nus 1/2-inch material, and ended at the board surface at the toe of the 

pe. As the material had been placed in the same manner as for the other 
t, except that it was damp and the sand remained coated to the rock, itis 
ieved the sand reduced the friction between rocks in the top layer and the 
ce of the water flowing out of the blanket when a wave receded down the 

pe was enough to start the slide. 

After the slippage, the material was removed, rescreened in a dry condi- 
n, and replaced in the same two layers, as for Test 5. The material placed 
this dry condition was called Test 6. The model was then backfilled to a 
th of 2.7 feet, the wave machine was set as before and operated for 24 

irs. No slippage occurred, and the fine material appeared stable; there 

5 little movement of sand. Some of the surface material (larger than 1/2- 
h) moved, forming a beaching slope. This 1/2-inch and larger material 
1a tendency to move more readily than the 3/4-inch and larger or the 1- 

h and larger, which had been tested previously. The photograph, showing 
condition of the model after 23 hours of operation, is shown in Figure 15a. 


sts 7 and 8 


sts to determine effect of sand on slippage failure 

Tests 7 and 8 were conducted in an attempt to determine the effect of the 
id when it remained coated to the rock, causing failure. For Test 7, the 
vel material was screened while damp on the 3/4-inch screen and placed, 
in Test 2. A 4-inch layer of material smaller than 3/4-inch in size was 
vered by a 4-inch layer of material larger than 3/4-inch in size. The only 
erence between Tests 2 and 7 was that the material in Test 7 was screened 
p, thereby leaving a coating of sand particles on the larger gravel parti- 
s. The flume was filled to a depth of 2.7 feet, and the wave machine was 
rted at 42 waves per minute, resulting in a wave height of 0.43 foot. 
Immediately on starting the wave action, considerable slippage occurred. 
slippage was less pronounced, however, than when the material was 
arated on the 1/2-inch screen. The coarse ‘material slipped near the sur- 
e until it had decreased in thickness by approximately 1 inch. After 
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Figure 14 


(b) Slippage that occured immediately “ 
after test started 


Yakima Project - Washington 
, KENNEWICK MAIN CANAL 
TEST 5 (GRAVEL SEPARATED ON 1/2-INCH SCREEN) 
Model study 
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(a) Test 6--After 23 hours of operation 
2.7 ft. depth--42 waves per minute-- 
wave height 0. 43 ft. 
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(b) Test 9--Prepared for operation-- 
Recommended cover - talus 
separated on 3/4-inch screen 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TESTS 6 and 9 
Model study 
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roximately 1 minute of operation, the material then became stable. Fur- 
r operation produced little change. The test appeared to confirm the fact 

t the sand coated to large material contributed largely to the slippage 

ure in Test 5. 

For Test 8, the gravel was placed, as in Test 1, well mixed, but instead of 
ng in a damp condition, it was dry. An 8-inch layer was placed on the 2:1 
pe, the model was filled to a depth of 2.7 feet, and the wave machine 

rted at 42 waves per minute, creating 0.43-foot wave height. Approximately 
aches of beaching at the water surface occurred in the first 2 minutes, and 
2r about 3-1/2 inches of beaching had occurred, the gravel appeared to be 
te stable. A photograph showing the condition of the cover after 5 hours of 
ration is shown in Figure 19. The dry gravel resulted in a more stable 
dition than the damp gravel of Test 1. 


st 10—Talus Material Moistened and Well Mixed 


A 6-inch layer of talus material moistened and well mixed was placed on a 
slope in the testing flume. The method of placing was similar to Test L, 

which gravel material was used. The flume was filled to a depth of 2.7 

t, and the wave machine started at 32 waves per minute. There was move- 
nt of the finer talus material, but not as much as that observed in the grav- 

material of Test 1. After 20 minutes of operation, the frequency of the 

ves was increased to 42 per minute, and after an additional 20 minutes to 
per minute, with a resulting wave height of 0.32, 0.43, and 0.5 foot, respec- 

ely. At the higher frequency and wave height, there was considerable ero- 
n of the talus material which was washed down to form a beaching slope of 
6:1. The model was operated for 2-1/3 hours at 62 waves per minute. 

A photograph showing the condition of the talus material after 3 hours of 

al operation, described above, may be seen in Figure 16a. After the beach- 

‘slope had been recorded, the model was left running at 62 waves per 

nute, and the water was raised slowly to a depth of 3.2 feet, then lowered 

a depth of 2.2 feet. Raising and lowering the water surface caused various 

ve heights for 62 waves per minute. The scour observed showed consider- 
e material had moved down the slope, but at no point was the 2:1 slope ex- 

sed. A photograph showing the condition of the model after raising and 

rering the water surface is shown in Figure 16b. The model was operated 
a total of 5 hours. Test 10 indicated that the angular talus material was 

\siderably more stable than the rounded gravel material placed in the 


née manner. 


st 11—Talus Separated on 1/2-inch Screen 


The talus material was dried and screened on the 1/2-inch screen, result- 
in the analysis curve shown in Figure 17. Four inches of the material 
er than 1/2-inch in size was placed on the 2:1 slope and covered by 4 
es of material larger than 1/2-inch in size. The model was filled to a 
th of 2.7 feet, and the wave machine was started at 32 waves per minute, 
ng a wave height of 0.32 foot. At 32 waves per ‘minute, there was very 
e movement of the surface material. The material in the bottom layers 
er than 1/2-inch in size also appeared very stable. 
As no damage occurred at 32 waves per minute, the frequency of the waves 
increased to 42 and then to 62 waves per minute, and the model was 
rated for a 2-hour period. The resulting wave height was 0.43 and 0.5 foot, 
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Figure 16 


(a) After 3 hours operation at various 
wave heights and frequencies 


(b) After 5 hours operation with 
various waves and depths 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 10 (TALUS MOISTENED AND WELL MIXED) 
Model study 
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respectively. At the higher wave frequencies, the surface material moved 
more freely, and approximately 1-inch was eroded in forming the beaching 
slope. However, the material in the bottom layer below the 1/2-inch in size 
still appeared very stable and very little movement occurred. 


Test 9—Recommended Cover--Talus Separated on 3/4-inch Screen 


The mixed talus material was screened on the 3/4-inch screen, together 
with the base material analysis curve and the Bureau of Reclamation filter 
criteria. The resulting analysis curve is shown in Figure 18. Before placir 
the talus material on the 2:1 slope, 3.6 percent moisture, by dry weight, was 
added to the material smaller than 3/4-inch in size to simulate field condi- 
tions. Four inches of the material smaller than 3/4-inch in size was placed 
on the 2:1 slope and covered by 4 inches of material larger than the 3/4-incl 
size, Figure 15b. The model was filled to a depth of 2.7 feet, and the wave 
machine was started. The model was operated at a frequency of 32 waves 
per minute, creating 0.32-foot waves for 15 minutes. As the material ap- 
peared to be very stable, the frequency of the wave machine was increased t 
42 waves per minute. The model at this frequency and wave height (0.43 foo 
was operated for 30 minutes. There was a small amount of movement of the 
fine material. The wave frequency was then increased to 68 waves per min- 
ute, with a resulting wave height of 0.5 foot. 

The 0.5-foot wave approaching the slope may be seen in the photograph, 
Figure 20a. At 68 waves per minute, the surface material was still relative 
stable. There was a little movement of the fine material in the lower layer; 
however, this movement in the lower layer appeared to stop after about 2 
hours of operation. The model was operated at 68 waves per minute for ap- 
proximately 5 hours. The photograph, shown in Figure 20b, gives the condi- 
tion of the slope at the end of the 5-hour period. A comparison of Figures 1 
and 20b shows very little change in the appearance of the material. ‘ 


Test 12 (Fine Base Material) 


A 6-inch layer of the fine base material (analysis curve shown in Figure 
was subjected to wave action. In preparing the material, enough water was” 
added to bring the moisture content up to 18 percent of the dry weight of the 
material. The 18 percent moisture is the optimum moisture content for 
maximum compaction. The material was placed on the 2:1 slope in layers ¢ 
approximately 2 inches, and each layer was tamped to near maximum com= 
paction with a wooden Caper Before placing the base material on the 2: 1 
board slope, it was roughened with a layer of metal lath to prevent any ten- 
dency for the base material to slide on the plywood board slope. 

The water was raised slowly in the flume to a depth of 2.7 feet, Figure 2 
There was some settlement of the material on the slope, and near the wate1 
surface, cracks could be observed. The wave machine was started and 
operated at 32 waves per minute. The wave height varied from 0.25 to 0.4 
foot as the material washed down the slope. During this operation, comp] 
failure of the base material slope occurred. Immediately on starting the ’ 
wave machine, the earth material began to fail. The photograph taken 2 
minutes after operation started may be seen in Figure 21b. The density c 
rent of fine material flowed down the slope to the bottom of the flume. Aft 
12 minutes of operation, the 2:1 slope was completely exposed. The model 
was operated for 1 hour, after which it was turned off and slowly drained. 
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Figure 19 
Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 8 - AFTER 5 HOURS OF OPERATION 
Model study 


," 


The fine base material assumed an approximately level position, as shown 
Figure 22a. . 


Test 13—Base Material and Plus 1/2-inch Talus Material . 


A 4-inch layer of the base material was compacted on the 2:1 slope and 
covered with 4 inches of talus material larger than 1/2-inch in size, Figui 
17. The fine base material was placed at optimum moisture content (18 p 
cent by dry weight) and compacted in thin layers with a wooden tamper. 
talus material was dumped on the compacted base material, simulating th 
way that it would possibly be placed in the field, Figure 22b. The model w 
slowly filled to a depth of 2.7 feet, and the wave machine was started at 3 
waves per minute, giving a wave height of 0.32 foot. 

Immediately after the model was started, the fine base material besa” 
leaching through the talus cover, the material being leached out formed a 
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Figure 20 


by 


») Slope after 5 hours of operation at 
_ varying wave heights and frequencies 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 9 (RECOMMENDED COVER - 
TALUS SEPARATED ON 3/4-INCH SCREEN) 
Model study a 


138 


May, 1959 HY 


rekon) Bs EAE Figure 21 


(a) Ready for operation 


(b) Two minutes after operation started 
at 2.7 ft. depth and 32 waves per 
minute--wave height 0. 26 ft. 


Yakima Project - Washington . 
KENNEWICK MAIN CANAL 
TEST 12 (FINE BASE MATERIAL) 
Model study 
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Figure 22 


(a) Test 12--after 1 hour of operation 
at 2.7 ft. depth--32 waves per 
minute--wave height 0. 26 ft. 


wie — 


“NW35-S7_ 


(b) Test 13--ready for operation. = 
Base material covered with 
plus 1/2-inch talus 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TESTS 12 and 13 _~ 
Model study a 
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density current which moved down the slope in the water immediately above 
the talus cover and continued on to the flume floor. The density current ma 
be seen in the photograph of Figure 23a which was taken 7 minutes after the 
model was turned on. After 30 minutes of operation at 32 waves per minute 
the wave machine frequency was increased to 42 waves per minute. The in- 
crease in frequency increased the wave height to 0.43 foot. At this increase 
frequency, the leaching of base material increased. 

A photograph taken 90 minutes after operation started, shows the base 
material had washed within approximately 1/2 inch of the board slope (Figu: 
23b). The model was operated an additional 17-1/2 hours at 42 waves per 
minute. During this time, the upper part of the slope failed completely. Th 
finer material was completely washed from beneath the talus cover. The 
photograph showing this failure may be seen in Figure 24a. 


Test 14—Recommended Cover and Base Material 


A final test was run using the talus material as protective cover for the 
fine base material. The moisture in the base material was increased to op: 
timum moisture content, and a layer 4 inches thick was placed by the tampi 
operation explained previously. The talus material was screened on the 3/: 
inch screen. Four inches of talus material smaller than 3/4-inch in size w 
placed by dumping over the fine base material, and on top was placed 4 inck 
of talus material larger than 3/4-inch in size, Figure 24b. The analysis 
curves are shown in Figure 18. The total thickness of the test section was 
12 inches. The model was filled to a depth of 2.7 feet and operated for 28 
hours at a frequency of 42 waves a minute, giving a wave height of 0.43 foot 

As shown by Figure 25, the test section withstood the wave very well, a1 
very little movement occurred in the talus material. No movement occurr¢ 
in the fine base material. To check the stability of the section, the model vy 
operated for a total of 128 hours. After the first 28-hour period, the opera 
tion was as follows: 24 hours at 3.2 foot depth, 42 waves per minute, 0.43 
foot wave height; 24 hours at 2.2 foot depth, 42 waves per minute, 0.42 foot 
wave height; 4 hours at 2.7 foot depth, 22 waves per minute, 0.6 foot wave 
height; 24 hours at 2.7 foot depth, 51 waves per minute, 0.62 foot wave heig 
and 24 hours at 3.03 foot depth, 22 waves per minute, 0.6 foot wave height. 

As may be seen from a study of Figure 26a, the material withstood the — 
operation very well. At the frequency of 51 waves per minute, there was © 
some movement of the surface material. During the complete test, the talt 
material below 3/4-inch in size was very stable, and no leaching of the fine 
base material was observed. The waves of 22 per minute frequency ran ts 
and down the slope from maximum to minimum, a distance of 2.2 feet. Thi 
0.6-foot high waves at 22 waves per minute resulted in the most damaging 
condition where leaching action was concerned. However, no leaching of th 
base material occurred. 4 


The section was operated to failure after the above test, the water surfé 
was raised to a depth of 3.03 feet, and the wave machine set for 45 waves } 
minute, with a resulting wave height of 0.75 foot. For this condition, the n 
terial failed. Failure started on the surface and progressed to the bed ma 
terial. Total failure occurred in 12 minutes, as may be seen from the pho 
graph of Figure 26b. The large waves were large enough to wash the larg 
gravel particles off the slope, exposing the finer gravel and eventually the 
fine base material. 
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(a) After 7 minutes of operation at 
32 waves per minute--depth 
2.7 ft.--wave height 0. 26 ft. 


4 


q 
/ 
; 
4 
: 
| 
4 


(b) After 90 minutes of operation 

: at 32 and 42 waves per minute 
| leaching of base material 
progressed as shown by lines 
on the photograph 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 13 (BASE MATERIAL AND 
PLUS 1/2-INCH TALUS MATERIAL) 
Model study 


Figure 23 
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(a) Test 13--after the model was left 
running all night at 42 waves per 
minute. Base material completely 
leached out in area near water 
surface 


(b) Test 14--ready for operation. 
Base material covered with 
talus separated on 3/4-inch 
screen - recommended cover 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TESTS 13 and 14 
Model study 
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Figure 24 
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Figure 25 
Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 14 (RECOMMENDED COVER) AFTER 28 HOURS OF OPERATION 
AT 42 WAVES PER MINUTE--2.7 FT, DEPTH--WAVE HEIGHT 0.43 FT, 
Model study 


Results and Conclusions of Hydraulic Model Tests 


The gravel used in the first test was well rounded, and sand adhered to the 
irgest particles when the material was moistened and mixed. When sand ad- 
ered to the larger particles, failure of the 2:1 slope by a slip occurred. The 
ounded particles of the gravel moved more readily when exposed to wave 
ction than did the same mean size angular particles of the talus material 
js in Tests 9, 10, 11, 13, and 14. - 

The tests indicated it was desirable to screen the cover material and place 
_in two layers, the larger material protecting the finer material. The larger 

e average sized material on the surface layer, the more stable was the sur- 
ee to wave action and the more susceptible were the fines to leaching. The 
maller the average size of the surface material, the more susceptible the 


os was to movement by wave action. 
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Figure 26 


(a) After 128 hours of operation at 
various depths, wave heights 
and frequencies 


(b) Failure of slope from waves 0.75 ft. high and 


45 per minute--approximately 12 minutes 
after final test started 


Yakima Project - Washington 
KENNEWICK MAIN CANAL 
TEST 14 (RECOMMENDED COVER) 
Model study 
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The material that fell within the Bureau of Reclamation filter criteria 
mits appeared to provide an effective filter against leaching of fine materi- 
S by wave action. The study indicated that a protective material must be 
ne enough to prevent leaching of the base material through it and that par- 
cles on the exposed surface of the protective material must be of sufficient 
ze and weight to resist rolling down the slope from wave action, 

The study further indicated that the angular talus particles form a more 
able protection blanket than the rounded gravel particles. 

The results of Tests 9 and 14 showed the best protective cover was the 
lus material separated on a 3/4-inch screen and placed in 2 layers: (1) a 
yer of talus material smaller than 3/4 inch placed over the fine base ma- 
rial, and (2) a layer of talus material larger than 3/4 inch placed on the 
irface. The resulting analysis curves in the Bureau of Reclamation filter 
‘iteria are shown in Figure 18. The test showed that the layers of cover 
aterial should not be less than 4 inches in thickness for adequate resistance 

erosion by small waves. 
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HYDRAULIC MODELS OF THE ST. LAWRENCE POWER PROJECT 


John B. Bryce! 


ABSTRACT 


Hydraulic models were used extensively in the design of many hydraulic 
atures of the St. Lawrence Power Project and in the development of a plan 
yr controlling the river during construction. A description is given of the 
roject and the models, and some of the more important findings are dis- 
issed. 


INTRODUCTION 


In the design and construction of the hydraulic features of the St. Lawrence 
ower Project, a fundamental part was played by hydraulic models. In all, 
even models were constructed for the Power Project, some being compre- 
snsive river models and others reproducing complete structures or parts of 
ructures. The models were constructed and tested at the Hydraulic Model 
aboratory of the Hydro-Electric Power Commission of Ontario for the use 
‘the Commission and the Power Authority of the State of New York, who 
ere jointly responsible for the design and construction of the Power Project. 
the models, the final detailed design of all channel enlargements in the 
ternational Rapids section was developed. Also evolved and demonstrated 
as a satisfactory plan of river control which would permit the construction 
‘the channel enlargements and the structures without adversely affecting 
ver levels and flows, and without interfering with existing navigation. In ad- 
tion, the optimum location of temporary and permanent structures was indi- 
ited, and assistance was given in predicting flow and velocity conditions dur- 
g construction. The hydraulic design of certain features of the structures 
as also developed. In the Power Project Model investigation, many interest- 
g and complex problems were encountered, It is the purpose of this paper 
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to describe the models, and to outline the major design problems and the 
methods used in their solution, 


Brief Description of Project 


As shown in Fig, 1, the St. Lawrence Power Project extends throughout the 
International Rapids section of the St. Lawrence river from Ogdensburg, New 
York, and Prescott, Ontario, downstream to Cornwall, Ontario, a distance of 
about 45 miles. In nature a fall of approximately 92 feet occurred in this 
distance, the main concentrations being at: the Galop Rapids, 9 feet; the Plat 
Rapids at Ogden Island, 12 feet; and the Long Sault Rapids, 30 feet; the re- 
mainder of the fall occurring at lesser concentrations and as river slope. The 
flow in the International Rapids Section is remarkably uniform as a result of 
the natural regulating effect of the Great Lakes. The long term average flow 
is slightly greater than 240,000 cfs, while the maximum and minimum monthly 
mean recorded flows between 1860 and 1954 were 315,000 cfs and 166,000 cfs 
respectively. 

The principal features of the Power Project may be briefly described and 
summarized as follows: 


The Powerhouse 


Located at the outlet of the channel around the north side of Barnhart 
Island, a short distance upstream from Cornwall, the powerhouse is some 
3,120 feet long, about 160 feet high and contains 32 generating units with a 
total rated capacity of about 2,500,000 hp at 81.0 feet head. Six ice sluices 
also are provided. Together with the Long Sault Dam and 13 miles of dykes, 
this power dam enables the level to be raised over 80 feet at this point. Down- 
stream from the powerhouse to Polly’s Gut, channel enlargements are being 
undertaken to lower tailwater levels. 


The Long Sault Dam 


The Long Sault Dam is required to contain the power pool by closing the 
main river channel south of Barnhart Island. This dam provides the spillway © 
for the project and must be capable of discharging the full flow of the river. 
Thirty sluices, each 50 feet wide, are provided, and each with a gate 52 feet 
wide and 30 feet high. The dam is about 2,400 feet long and is about 110 feet — 
high at the deepest point. : 


The Massena Intake 


An intake structure is provided near the entrance to the Massena Power _ 
Canal to permit continuance of the flow to the Aluminum Company plant during 
construction of the project and domestic supply thereafter. The intake con- 4 
tains four sluices or ports each controlled by a gate. 

¢ : 

' 

Iroquois Control Dam a 
Located just upstream of Iroquois, Ontario, the Iroquois Control Dam en- ¥ 
ables the St. Lawrence River flows to be controlled and regulated. The dam ~ 
_ is some 2,340 feet long, 70 feet high, and contains 32 sluices, each 50 feet 
wide, and each provided with a gate 52 feet wide and 48 feet high. 4 
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Channel Enlargements Morrisburg to Sparrowhawk Point 


Extensive channel enlargements are required in this reach of the river to 
provide velocities not exceeding those at which an ice cover may form by 
packing during winter operation. These enlargements also must provide an 
adequate course for 27-foot navigation. 


Channel Enlargements Cardinal to Chimney Point 


In this reach, very extensive channel enlargements are required to provide 
a deep-draught navigation channel and velocities not exceeding those stipulated 
for navigation. 


Navigation Features 


While not directly connected with the Power Project, certain navigation 
features in the International Rapids section are interrelated and will be 
mentioned here. The Canadian St. Lawrence Seaway Authority has improved 
the channel north of Cornwall Island and portions of the channel south of 
Cornwall Island. The United States Seaway Development Corporation has en- 
larged the remainder of the south channel around Cornwall Island and has 
constructed the Snell Lock, the Eisenhower Lock and the Wiley-Dondero Canal 
upstream from the locks. Twenty-seven-foot navigation will enter the head- 
water pool by this means and after leaving the Wiley-Dondero Canal will sail 
up the flooded river channel to Morrisburg and the enlarged river channel to 
Iroquois. At Iroquois, the St. Lawrence Seaway Authority has provided a lock 
to by-pass the Iroquois Dam. After passing through the lock, shipping will 
proceed up the enlarged channel to Chimney Point and the natural river there- 
from. The location of these navigation features may be noted on Fig. 1. 


Scope of Model Investigation 


The hydraulic problems associated with the St. Lawrence Power Project 
were many and varied, and all were accentuated by the large scale of the 
Project and the continuously large flow in the river. For many of these 
problems it was considered that the construction and testing of hydraulic 
models would yield the best, and in some cases the only reliable method for | 
their satisfactory solution. Accordingly, it was decided to carry outa | 


comprehensive program of model testing to provide the maximum assurance _ 


‘ 


that the Project would be constructed and would operate as planned. Also % 
from past experience, it was considered probable that substantial economies : 
might be realized in view of the large scale of the project. Construction of } 
the first St. Lawrence models commenced in 1953 and testing began later that 
year. From that time throughout the construction of the project building and 
testing of the models proceeded continuously. As final designs were complet- 
ed, the emphasis shifted to the steps in construction and the provision of ; 
information useful to construction. 


In general, the problems to be investigated fell into four categories: : 


(1) The design of channel enlargements, ; 

(2) Problems at structures, such as their optimum location, the design of 
their hydraulic elements, and the provision of adequate discharge 
capacity both during and after construction. 
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(3) The development of a plan of river control during the construction of 
the Project. 

(4) The determination of velocities and currents in the navigation channels 
throughout the Project. 


For the investigation of the problems eleven hydraulic models were construc 
ed. Four of these models were distorted-scale river models and reproduced 
altogether over 40 miles of the St. Lawrence River. The remainder were 
natural scale models, three of these being large comprehensive models of 
principal structures. With regard to the large river models, it might be ob- 
served that one of the principal reasons for their construction was that in the 
International Rapids section many islands divided the flow and abrupt change 
in direction occurred at some critical points. It was considered that by no 
other means could reliable information be obtained under flooded conditions 
of the flow distribution through the multiple channels, the velocity distributio 
across the sections, the actual velocity in the navigation channels, and the de 
tailed effects of modification in channel enlargements, 


Collection of River Information 


To reproduce over 40 miles of the St. Lawrence River in the models, it 
was essential that accurate river bed contours and topographic data be ob- 
tained for their construction, and that river levels, flow distributions and 
currents be available for model verification. The collection of such data ove 
the length of such a large river posed a major problem. 

Fortunately, earlier surveys of the International Rapids section provided 
a general coverage of subaqueous contours and other topographic data and 
these were used as a base. To verify the information, echo soundings were 
carried out at sections about one mile apart over the length of the river to be 
reproduced in the models. Where discrepancies occurred, or in critical 
areas, extensive programs of additional echo soundings were carried out. In 
the Long Sault Rapids section, soundings were made from helicopters as the 
most feasible method. Aerial photographs of the whole reach were made to 
provide accurate mapping. 

With regard to hydrometric information, records from several long term 
lock gauges were available, and these were augmented by over 30 automatic: 
recording gauges which provided a more general coverage, However, for — 
verification purposes a much denser coverage was required and this was ob- 
tained by establishing over 100 temporary gauges. Intensive reading of these 
gauges at selected flows provided the necessary information. Fortunately, — 
these gauges were established at the time of the peak flow in 1952, and as the 
flow subsequently subsided a considerable range in levels was obtained. 3 
Current meterings to determine velocities and flow distribution around the _ 
islands were made, and float observations indicated the current pattern. 


Islington Model Testing Laboratory 


a PND een ey ae 


All the models for the St. Lawrence Power Project were constructed and 
tested at Ontario Hydro’s Hydraulic Model Laboratory at Islington, just west 
of Toronto, Ontario. This laboratory, constructed specifically for model te 
ing, was begun originally for Niagara redevelopment model tests and was 
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‘ended for the St. Lawrence investigation. The building is about 500 feet 

i and 70 feet wide over above three-quarters of its length. Over the re- 
‘ining one-quarter, the width is doubled to 140 feet, as shown on Fig. 2; 

ich also indicates the plan of the laboratory and the general arrangement 
the models. The laboratory is serviced by a travelling gantry crane over 
yst of its length, and accurately levelled rails provide vertical control from 
welling point gauges. Two constant-level head tanks provide the flow for 

: models through distribution mains; and sumps generally encircling the 
ilding provide the return for recirculation. A special use of Bailey bridging 
yterial provides uninterrupted floor space throughout. 


The River Models 


The four river models constituted a large and important part of the St. 
wrence model test investigation. As these models have many features in 
mmon, they will be described and discussed collectively rather than indi- 


ually. 


scription of River Models 


In Fig. 1 is shown a plan of the International Rapids section on which is 
arked the limits of the individual models. Model No. I extends from 
-densburg to Leishman’s Point and reproduces about 16 miles of the river. 
sis model contains the Galop rapids reach and the Sparrowhawk Point — 
oquois reach where very large channel enlargements were required, and in- 
ides also the Iroquois Dam and the Iroquois 27-foot navigation lock. Model 
». 2 reproduces approximately 8 miles of the river from above Point Three 
»ints to below Canada Island and contains the Ogden Island reach where im- 
rtant channel enlargements were necessary. Model No. 3 extends from 
ove Croil Island to Cornwall, a distance of about 14 miles, and contains the 

ng Sault Dam and Power House. Model No. 3 was used for the dewatering 
udies of the Long Sault Dam, and it was considered necessary to extend the 

del upstream beyond Croil Island where the river is ina single channel in 
der to obtain the right flow distribution in the multiple channels in the Long 
ult area. In Fig. 3 is shown a photograph of Model No. 3 which is typical of 

various river models constructed. Model No. 9, reproduces the area from 
‘ove the Powerhouse to the upper end of Cornwall Island and includes the 
ilrace area where channel enlargements are necessary to obtain an eco- 
mic head regain at the plant. It will be noted that these four river models — 
iproduce the entire International Rapids reach from Ogdensburg to just above 
rnwall, with the exception of an 8 mile reach in the Goose Neck Island area 


nsidered desirable to construct several models rather than one long con- 
uous model, as work could proceed simultaneously on all models without 
tests on one reach interfering with tests on other reaches. 


Oice of Model Scales 

The scales chosen for river models 1, 2 and 3 were 1:500 horizontally and 
00 vertically. The choice of the horizontal scale ratio was largely dictated 
laboratory space requirements consistent with other requirements. A 
reful study was made of the optimum placement of the models in the 


= 
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Fig. 3. River Model No. 3 showing Powerhouse and Dykes in Foreground E 


laboratory, and as will be seen from Fig. 2, a larger horizontal scale would 
have required a very considerable increase in laboratory space. This scale 
was also considered satisfactory from accuracy and workability standpoints. 

Due to the horizontal scale selected and the large river width relative to 
its depth, distortion of the vertical scale was considered necessary for simi 
lar performance of model and prototype. The measurement of model veloci: 
ties and depths was of great importance in these models, and it was desired 
that the vertical scale be sufficiently large to enable these quantities to be — 
measured with the required precision. However, as the vertical scale in- S 
creased, the difficulty of producing the required model roughness also in- E 
creased. A distortion ratio of 5:1 was selected, thus producing a vertical 
scale of 1:100, which was considered the largest feasible in view of the roug 
ness requirements, and adequate from the measurement standpoint. From 
past experience this distortion ratio had proven successful in other models 
a similar nature. The scale distortion in these models required also that tt 
Iroquois Dam and the Long Sault dam structures be distorted as well. In te 
where distortion of these structures affected the results, their true head- — 
discharge relationship was determined from the larger natural scale model 
of these structures, and the correct levels produced in the distorted-scale. 
models. : 

The scales chosen for Model No. 9 of the power house tailrace were 1:16 
horizontally and 1:80 vertically. It was considered that for this model the 
scales should be enlarged and the distortion reduced, as a detailed econo r 
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dy was desired and as a considerable change in direction of the flow oc- 
‘red on entering the tailrace. 


istruction of River Models 


The river models were all constructed in the same manner using a method 
ich has been found to be satisfactory. Over the entire model area, thin 
wood templates were erected on a two-foot grid running both laterally and 
gitudinally. The tops of the templates were accurately cut to the surface 
vation of the ground or river bed. In the areas where no channel im- 
ements were expected, sand was placed between the templates to within 

) inches of their tops, and then a layer of concrete was placed on top of the 
id to produce the final surface. Between templates, a “router” on a travel- 
¥ point gauge traced the contours to the correct elevation. In areas where 
nnel changes were expected a different procedure was used. A concrete 
ple” was poured just below the minimum grade expected. A concrete block 
s then poured on top of the table with the surface finished as described 

ve. Provision was made for lifting out these blocks so that they could be 
laced with others containing desired improvements. In Fig. 4 is shown a 
‘tograph of the models under construction. 


del Roughness 


One of the major problems in the river models was the production of the 
juired bed roughness to reproduce accurately prototype levels and veloci- 
3. Artificial roughness beyond that provided by the concrete channel beds 
3 necessary due to the scale distortion of the models. Anticipating this 
jbolem, a “Roughness Model” was constructed and tested prior to the com- 
tion of the river models themselves. In Fig. 6 is shown a photograph of 

3 model which represented, to the same scales and distortion as the river 
dels, a channel 20,000 feet long, 2,000 feet wide and with a normal depth of 
feet. These dimensions were chosen as they were of the same order of 

th and width as the river model channels. In this model various means of 
ducing roughness were tested, including treatment of the concrete bed it- 
J the addition of various thicknesses of expanded metal screening, and the 
ition of one-half-inch wide vertical metal strips of various heights and at 
ious spacings. The resulting roughness was measured in terms of the 
fficient “n” in the Manning formula, and it is felt that a few of the results 
nined may be of interest. At the normal depth, the mean roughness coef- 
ents obtained for the various treatments were as follows: 


1) Brushed concrete; n = .015 

(2) Single thickness expanded metal screening; n = 0.028 

3) Double thickness expanded metal screening; n = 0.031 

‘4) Metal strips, 8 inch centres, height = 1/2 depth; n = 0.027 
5) Metal strips, 6 inch centres, height = 1/2 depth; n = 0.033 
6) Metal strips, 4 inch centres, height = 3/4 depth; n = 0.060 
'7) Metal strips, 2 inch centres, height = 3/4 depth; n= 0.116 


m an analysis of these measurements and others obtained at various depths 
as concluded that the metal strips would provide the best type of artificial 
hness for models 1, 2.and 3. It was found also that the roughness pro- 

d by the strips was a function of the strip spacing and the ratio of strip 

ht to water depth. As it was desired to have an unobstructed water surface 
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Fig. 5. River Models under Construction 


float measurement, it was decided that the basic roughness pattern would 
produced by one-half-inch wide strips spaced on a six-inch staggered grid, 
a height above the bed of one-half the depth of mean flow conditions. Actu- 
y, slots were provided in the model bed on a three-inch staggered grid to 
»vide flexibility in increasing or decreasing the roughness. 
When the models were verified it was found that the basic roughness pattern 
uld only be considered as a general guide, as the “form-factor resistance”, 
. that caused by the shape and alignment of the river channel and the 
»sence of islands, constituted an appreciable portion of the total resistance. 
such points, and where rapids occurred, it was necessary to appreciably 
the strip spacing to provide correct levels. In such cases the strips and 
slots on the three-inch grid proved to be a very satisfactory and flexible 
ans of adjusting roughness. For Model No. 9, with less distortion and 
er scales, it was found that metal screening over the concrete bed was 
‘icient to produce the desired roughness. 


lature Current Meter 
n all the river model tests, an accurate measurement of velocities in the 
el channels was of the utmost importance as the criterion for the channel 
rgement design was the production of a specified mean velocity. Similar- 
the determination of flow distribution through multiple channels and around 
nds depended upon the ability to measure velocities. To enable these ve- 
ties to be measured with sufficient accuracy an improved miniature 
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Fig. 6. Roughness Model Showing Metal Strip Test Pattern 


propeller-type current meter was developed. The requirements for the met 
were that it be as small as possible, offer minimum resistance to the water, 
register accurately model velocities between 0.1 feet per second and 2 feet 
per second, and be remote indicating. In Fig. 7 is shown a photograph of the 
current meter developed. Propellers of various pitches and having different 
numbers of blades enabled velocities as low as 0.07 feet per second and as 
high as 4.5 feet per second to be measured reliably. A feature of the meter 
is the remote electronic velocity indicator. Previous electrical indication o 
velocity required a commutator which inevitably produced a frictional drag. 
This was eliminated in this meter by devising a means of utilizing changes i 
electrical resistance of the water itself around the propeller. Basically the 
instrument consists of an electronic meter which counts the propeller ro- 
tations by means of the changes in resistance of the water between the blade 
and a small electrode, marked X in Fig. 7, located within 0.05 inch of the o7 
of the blade tips. A second electrode, marked Y and remote from the blade: 
was incorporated to balance any changes in water conductivity. An electron 
circuit counts resistance changes, and the number of changes in a given tim 
period is a multiple of the rate of rotation of the meter and a function of the 
water velocity. In operation, the meter is set up in the model channel and.a 
electric cable carried the impulses to the counter which is set up in a con- 
venient location. The operator determines the velocity from the electronic 
counter and an automatic timing device. A towing channel is provided in ‘i 
laboratory for periodic calibration of the meters. 
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Fig. 7. Miniature Current Meter with Electrodes for Frictionless Counting 


: 
Measurement of Surface Stream Lines and Velocities 


A technique has been developed at Islington for the production of surface 
tream lines from which surface velocities may be determined. The techniqu 

brief is that a camera is mounted on the laboratory roof and a vertical 
faerial” photograph is taken of the model below. Next, in semi-darkness, a 
umber of uniformly-spaced surface floats each containing a lighted candle 
.re released in the model channel at the upstream limit of the photograph. 
The camera shutter is simultaneously opened and the lighted floats trace 
urface stream lines on the previously exposed photograph of the channel. An 
pparatus, attached to the camera, provides a rotating segmented disc which 
mterrupts the field of view once every second. This produces a series of 
ashes instead of a continuous stream line, the length of the dashes being pro- 
ortional to the surface velocity. In Fig. 8 may be observed typical results 


f such a technique. 
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River Model Accomplishments 


The work on the river models extended over a period of several years, and 
in this paper it is practicable only to record the highlights of the results ob- 
tained. 


Design of Channel Enlargements 


The design of the required channel enlargements in the International Rapid 
reach was the largest single problem of the model investigation. From the 
model tests, detailed designs were developed which entailed the removal of 
some 64,000,000 cubic yards of material by land-based and dredging 
equipment. The criteria for the design of these enlargements were stipulated 
in the International Joint Commission’s Supplementary Order of Approval for 
the Project and are as follows: 


“Channel enlargements will be undertaken from above Chimney Point 

to below Lotus Island, designed to give a maximum mean velocity in 
any cross-section of the channel which will be used for navigation not 
exceeding four feet per second at any time, also between Lotus Island 
and Iroquois Point and from above Point Three Points to below Ogden 
Island, designed to give a maximum mean velocity in any cross-section 
not exceeding two and one-quarter feet per second with the flow and at 
the stage to be permitted on the first of January of any year, under 
regulation of outflow and levels of Lake Ontario in accordance with Plan 
of Regulation No. 12-A-9, as prepared by the International Lake Ontario 
Board of Engineers, dated 5 May, 1955. Downstream from the power- 
houses channel enlargements will be carried out for the purpose of re- 
ducing the tailwater level at the powerhouses.” 


In addition to these criteria, navigation interests required a channel having a 
minimum width of 600 feet with special widening treatment at bends, and with 
minimum depths of 29.5 feet and 28.5 feet upstream and downstream from 
Iroquois Dam respectively. 

The general nature of the required channel enlargements was first de- 
termined by computation, and these preliminary designs were then placed in 
Model No. 1 and Model No. 2. The initial controlling levels were estimated 
from backwater computations. These designs were then tested in the models 
and the velocities compared with the criteria. A very intensive series of 
tests then were undertaken to refine these designs and test alternatives, to 
obtain the most economical channels which would satisfy the governing cri- | 
teria. As the final design was approached, the critical water surface profiles 
corresponding to the design criteria were also obtained from the model. The 
design that emerged from these tests is shown in Fig. 1, where hatching indi- 
cates the areas of channel enlargement. It is of interest to note that in the 
Galop Island reach, Model No. 1, improvement of the South Galop channel was 
found to be of considerable economic value in reducing velocities and exca- 4 
vation in the main Galop enlargement which contained the navigation cheaaee 
Similarly at Ogden Island, Model No. 2 proved that economies amounting to — 
many millions of dollars could be effected by diverting more flow to the oa 
of the island by dry excavation, and thus reduce costly and difficult wet exca- 
vation in the Plat rapid channel to the north. In the Sparrowhawk Point - __ 
Toussaints Island area it was found that the alignment of the river was such 
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lat excavation somewhat in excess of that required by the criteria was neces- 
ary to provide satisfactory ice-packing velocities. The navigation channel 
as located in the model to conform with the requirements, and velocity 
1easurements were made throughout to ensure compliance with the criteria. 
loat pictures indicated the direction of the currents. Comprehensive water 
urface profile tests indicated the water levels and depths that would occur 
roughout the Project. In Model No. 9 a very complete model analysis indi- 
ated the economic tailrace enlargement, shown also in Fig. 1. 


roquois Dam Relocation 


In conjunction with an investigation of the foundation conditions for the 
roquois Dam, a series of tests was run in Model No. 1 to establish the opti- 
1um location and design of the dam from the hydraulic standpoint. Tests were 
1ade first on the dam at the location and to the design originally shown on the 
lans approved by the International Joint Commission in its Order of Approval. 
‘he dam as shown in these plans was located at Iroquois Point, was curved in 
lan, and had 40 sluices each 50 feet wide. The tests indicated unsatisfactory 
onditions, as an abrupt change in direction of the flow occurred at the dam 
1aking a considerable portion of the sluice area ineffective. After testing six 
ther plans, it was concluded that much more satisfactory conditions would 
esult with the dam moved downstream to Point Rockway and with a straight 
ather than curved alignment. At the new locations, channel enlargements 
rere necessary upstream to meet the ice-packing velocity requirement. 
fuch-improved hydraulic conditions were found with the dam at the down- 
tream location and with the associated channel enlargements. In Fig. 8 are 
hown stream lines with the dam in the original and final locations. The 
brupt change in direction of the flow at the dam has been eliminated, thus 
ermitting a reduction in the number of sluices from 40 to 32 while providing 
n even greater discharge capacity. Many other advantages also resulted; 
alf the dam could be built in the “dry” at Point Rockway thereby greatly 
implifying the considerable dewatering problem and the effect on upstream 
avels at the original location; head losses in the Iroquois reach were sub- 

ntially reduced resulting in a benefit of many millions of dollars in ultimate 
ower generation; the foundation conditions were much more favourable at the 
ew location, which, together with the fewer number of sluices, resulted in a 
reatly reduced dam cost and an overall saving of close to a million dollars 
n total construction cost of dam and channel enlargements. 


/ 
iversion Channels for Long Sault Dam 


_ During the first stage of construction of the Long Sault Dam, it was planned 
the southern portion extending from the United States mainland to Long 

ult Island would be constructed, and the South Long Sault channel diverted 

rough “Cut C” shown in Fig. 1. During the second stage, the North Long 

ult channel would be cofferdammed just upstream from the Long Sault 

ds, and the flow diverted through “Cut F” diversion channel and through 

ep diversion sluices left in the first stage of the dam. The design of these 

version channels was carried out in Model No. 3. The criteria for their de- 
n was that sufficient discharge capacity must be provided at all times by 

e channels so that levels upstream would be satisfactory for the continuance 

14-foot navigation and for the diversion of water through the Massena powe1 
nal, In the design of “Cut F” a very thorough study was made to develop its 
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Fig. 8. Surface Stream Lines with Original and Final Locations for 
Iroquois Dam 
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19st economic location and dimensions. This channel, which had to carry 
imost the whole St. Lawrence River flow, involved the removal of about 
000,000 cubic yards of material. In Model No. 3 studies were made of the 
equence and timing of the complete series of construction operations re- 
uired in the transitions from one stage into another, and a procedure was de- 
eloped that proved most satisfactory. 


iver Control 


By the terms of the Order of Approval by the International Joint Com- 
1ission it was required that the Project be constructed without adversely af- 
scting either upstream or downstream interests and that existing 14-foot 
avigation continue without interruption throughout the construction period. 

1 effect this required the general maintenance of pre-project levels and out- 
ows from Lake Ontario and pre-project levels in the International Rapids 
ection at the locks and in the river used by navigation. These requirements 
resented a considerable problem in developing an overall plan which would 
oordinate the construction of the channel enlargements and structures to pro- 
uce these conditions. In the development of this plan, the three large river 
1odels, Nos. 1, 2 and 3, were used. The large natural-scale models of the 
-oquois and Long Sault Dams, Models 8 and 5, were also used in determining 
1 detail the effects of construction at these dams. Utilizing all these models, 
step-by-step chronological program of construction was carried out in the 
1odels and the effects of the various operations noted. In the earlier stages, 
1e timing of various operations was determined such that the effect of one 
peration would be balanced by the effect of another. In the later stages, the 
ontrol required at structures was detailed. A coordinated plan was developed 
) ensure satisfactory level and flow conditions throughout, and this plan 
roved most satisfactory. 


The Long Sault Dam Models 


The Long Sault Dam was constructed in three stages. During the first 
‘age, that portion containing the southern 13 sluices was built with openings 
ft generally to river bed elevation in each sluice to pass the full river flow 
1 the second stage. The permanent sluice gates were adapted to act as di- 
srsion gates in these openings for control during the second stage. In the 
scond stage, the remainder of the dam was built, leaving 34 tunnel ports, 
ach controlled by a gate, to pass the full river flow during the third stage. 

1 the third stage the ogee sections were completed in the thirteen southern 
uices. Following this operation, the port gates were closed to raise the 
ater to the operating level, and the flow controlled by the permanent sluice 
ates. ' 

In the hydraulic design of the dam, several problems arose which were 
msidered to require model investigation; these were: 


(1) The need for energy dissipation downstream from the dam during both 
construction and ultimate operation periods, and operation periods, and 
if required, the design of an adequate energy dissipator. 
_ (2) The hydraulic loading on all the gates to ensure adequate hoisting ca- 
, city. 
(3) The Hitec capacity of the sluices and ports during the various 
construction phases. ee = 


164 May, 1959 HY 5 


Description of Models 


The Models were built to investigate these various problems. The first, 
Model No. 10a, was a sectional model of the dam consisting of two piers locat 
ed to form a central full sluice with a half sluice on each side. The ogee roll- 
way was constructed in such a manner that it could be converted to all stages 
in construction, i.e. diversion sluices, closure ports, or final ogee profile. 
This model, shown in Fig. 9, was constructed of plexiglass to a natural 1:50 
scale and tested in a plexiglass sided flume, 30 inches square and 26 feet long 
In this model the energy dissipation problem was investigated and the dissi-~ 
pator design developed, and exact models of the various gates were made and 
tested. The results of these gate tests are reported in a later section. 

The second model, No. 5, was a comprehensive model of the Long Sault 
Dam and reproduced, at a natural scale of 1:80, the entire Long Sault Dam, 
parts of Long Sault and Barnhart Islands, and a portion of the United States 
mainland. The river channel reproduced extended some 1,000 feet upstream 
and 2,500 feet downstream from the dam. The model was constructed entirel) 
in concrete and the dam was made with removable sections such that all 
stages in construction could be reproduced. The energy dissipator developed 
in the sectional model was also constructed downstream from the dam. A 
photograph of the model is shown in Fig. 10. The primary purpose of this 
model was to check the efficacy of the energy dissipator developed in Model 1{ 
over the full range of operating conditions in a comprehensive setting. From 


Fig. 9. Sectional Model of Long Sault Dam Showing Gate Load 7 
Measuring Equipment 7 
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st experience, it had been found that in many cases tailwater similarity 
mnot be achieved in sectional models due to the restricted channel repro- 
iced. This model was also used to determine the overall discharge ratings 
‘the diversion sluices and ports which were affected by high velocity of ap- 
‘oach conditions. 


nergy Dissipating Works 


From the model investigation a very severe energy dissipation problem 
2S found to exist at the Long Sault dam both during construction and under 
nal operating conditions. Because of the nature of the river bed, which 
opes downward rather steeply downstream from the dam, the tailwater depth 
as found to be low under all operating conditions, and insufficient to form a 
draulic jump. Preliminary tests indicated that, without energy dissipating 
orks, the lack of adequate tailwater depths resulted in downstream velocities 
hich would endanger the safety of the structure even with the relatively low 
‘adwater levels that would exist during the construction periods. The design 
‘oblem was unique in that the dissipating works must provide protection for 
ree entirely different flow conditions which would be encountered during 
ystruction and under final operation of the dam. 

After many trials in the plexiglass-sided flume, where over thirty different 
‘Signs were tested, the Recommended Design was developed, and was veri- 
2d in the comprehensive Model No. 5. It was found that this dissipator main- 
ined a forced hydraulic jump under all flow conditions, and that velocities 

the order of 75 feet per second with no dissipator were reduced to 


— Fig. 10. Comprehensive Model of Long Sault Dam 


166 May, 1959 HY! 


approximately 27 feet per second. In Fig. 11 are shown sketches illustrating 
the three flow conditions under which the dissipator had to operate. 

Although the Recommended Design was the best developed, rock condition: 
downstream from some of the sluices made it advisable to raise the bottom o 
the bucket by at least 5 feet. To meet this requirement the Revised Recom- 
mended Design was developed, and was shown to be almost as effective as the 
Recommended Design. Fig. 12 shows details of both energy dissipators. 


The Iroquois Dam Models 


In the final location selected for the Iroquois Dam, it was possible to build 
almost half of the dam in the “dry” on the United States shore. During this 
period, the channel enlargements on that shore were also completed in the 
“dry”, a rim dyke extending a short distance out into the river contained the 
construction area. When this phase was completed, the rim dyke was pro- 
gressively removed and the river diverted through the completed portion of 
the dam. Simultaneously, cofferdams were progressively extended out into 
the river to dewater the area for the second half of the dam. This work had 
to be carefully planned to maintain pre-project levels upstream. 

It was decided that the following problems during construction and oper- 
ation required model investigation: 


(1) The development of a plan of construction to provide sufficient dischar 
capacity during each stage of construction. Also required was a 
knowledge of the velocities and currents that would exist at the coffer- 
dams, particularly during the construction of the second stage coffer- 
dam. 

(2) The determination of energy dissipation requirements just downstream 
from the dam during the later stages of construction, when headwater 
levels would be raised and tailwater levels would be at or below natura 
due to the effect of downstream channel enlargements, and the design ¢ 
necessary energy dissipators. ; 

(3) The determination of the hydraulic forces acting on the sluice gates, te 
assist in selecting the required hoisting capacity. : 

(4) The development of a discharge calibration of the gates for constructic 
and operation conditions. The pattern of gate openings to provide the - 
most satisfactory river conditions was also desired. ’ 


Description of Models $ 


4 


5, 


q 
y 
5 
: 
4 
: 


dissipation requirements and to design the energy dissipator. Exact models 
of the sluice gates were also built and tested in the model to determine the - 
hydraulic loading on the gates. The results of this gate loading test are re- 
ported in a later section. : 
The second model, Model No. 8, was a large comprehensive 1:80 natural 
scale model of the Iroquois Dam, which included the whole dam and its as- 
sociated channel enlargements, and the downstream approach to the Iroquoi; 
Lock. The model reproduced the river upstream for some 5,300 feet and 
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mstream about 3,300 feet. The river portion of the model was constructed 
irely in concrete, while the dam itself was built of plexiglass. The coffer- 
n cells for the second stage cofferdam were also built for construction 

ts. In Fig. 13 is a photograph of the model showing the complete first stage 
onstruction and the cofferdams for the second stage. 

As in the case of the Long Sault Dam, it was desired to check the per- 
mance of the energy dissipator in a comprehensive setting. Also, the large 
del was required to determine in detail the velocities and currents to be 
ected during construction of the cofferdams during the dewatering of the 
ond stage area. The discharge calibration of the dam had also to be ob- 
1¢d from this model; as under final operating conditions the head across 
dam is very low and approach conditions affect the results. 


srgy Dissipating Works for the Iroquois Dam 


Under final operating conditions, the water-level differential across the 
quois Dam is quite small, and downstream velocities are low. However, 
ing one period prior to raising the levels in the forebay of the Barnhart 

ind powerhouse, the river was discharged through the sluices of the first 

[ of the dam with headwater levels raised and with tailwater levels at or 
shtly below natural. From model tests of these conditions, downstream ve- 
ities as high as 45 feet per second were measured, and the need for energy 
sipation for this period was demonstrated. Several schemes were investi- 
2d in the flume and the Revised Recommended Design was developed. By 
action, velocities downstream from the apron of the dam were reduced to 
roximately 17 feet per second. Additional tests indicated that the pro- 

ion of ten sluices would be adequate for this period, and the dissipator was 


Fig. 13. Comprehensive Model of Iroquois Dam 
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constructed to protect this number of sluices. Details of the design develop 
are shown in Fig. 12. 


Performance of Dam 


From Model No. 8 a complete discharge calibration was obtained for the 
dam during all stages of construction. This model, however, proved almost 
indispensable in predicting the optimum gate opening pattern to be used to 
produce satisfactory conditions downstream. This was particularly so in th 
period of construction when the headwater level was raised, but tailwater 
levels were still natural. Under these high head conditions, up to 15 feet, tk 
issuing velocities were necessarily high and unavoidably affected the veloci 
ties and currents downstream in the approaches to the Iroquois Lock. Fron 
the model, however, a pattern of part-gate openings was devised which pro- 
duced optimum possible conditions and which enabled navigation to proceed 
satisfactorily. 


The Massena Intake Model 


It was planned that diversion of flow through the Massena power canal 
would be continued during the construction period and also for a limited per 
after the raising of the St. Lawrence power pool. In the design of the intake 
it was necessary to ensure that there was sufficient gate capacity at normal 
river levels, and to determine the performance of the intake structure unde 
raised water level conditions. It was desired to investigate energy dissipat 
and to determine the hydraulic loadings on the intake gates under high-head 
conditions. 


Description of Model 


To study these hydraulic problems, Model No. 6, a comprehensive 1:60 
natural-scale model of the Massena Intake was built. This model reproduc: 
the complete intake itself, the excavated channel in which the intake is local 
and a portion of the existing power canal to ensure correct approach and ex 
conditions. The intake was constructed of plexiglass and the model bed of — 
concrete. Models of the emergency and service gates were also constructe 


Results of Model Study 


From the model, the best alignment of the structure in relation to the 
channels was determined to minimize scour in the downstream channel. 
Limited energy dissipation was found necessary, and this was achieved by 
ramping up the floor of the sluices near their outlet to direct the jet towar¢ 
the surface, and directing the jet by an extension of one pier. Interesting r 
sults were obtained from the gate tests and these are discussed in a later — 
section. 


The Powerhouse Models 


At the Powerhouse, three models were constructed to investigate the 
following problems: 


—"-=™ J esind en ee 
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(i) The hydraulic loading on the Headgates 

(2) The best location for the piezometer taps for the Gibson method of flow 
measurement 

(3) The hydraulic design of the ice sluices and the calibration of the ice 
sluices control gates. 


Jel of Powerhouse Water Passages 


Model No. 7, the Water Passages Model, reproduced the complete power- 
se water passage from headwater to tailwater. This model, to a 1:36 

iral scale, was constructed entirely of plexiglass to enable observations to 
made of internal flow conditions. Exact models of the headgates were made 
heet metal and installed in the model, and were used to determine the 

tical hydraulic loading on the gates. The results of the headgate tests are 
on in a later section. Fig. 14 is a photograph of the model which shows the 
siderable number of Gibson test piezometer connections tested. At this 
erhouse three large intakes supply each turbine, and the objective of the 

S was to determine the best cross-sections to be used for the Gibson tests, 
the best locations for the piezometer connections in each cross-section 
ecord the true average pressure. The model was found to be most useful 
his regard. 
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Model of Powerhouse Ice Sluices 


At the powerhouse, six ice sluices are provided to enable ice to be dis- 
charged downstream. Two adjacent sluices, 75-foot wide, are located at eac 
end of the powerhouse; while two adjacent sluices, 50-foot wide, are located 
at the middle of the plant. Each sluice is controlled by a drum-type gate. 
When operating, these sluices discharge ice and water through the powerhou 
from headwater to tailwater in a closed water passage which does not run fu 
To study the behaviour of ice being discharged through these sluices, and to 
develop a design which would operate safely, Model No. 4 was constructed 0: 
plexiglass to a natural scale of 1:36. Paraffin wax was used to represent ic 
and Fig. 15 shows a photograph of the model in action. Earlier designs tend 
to trap the ice inside the passage at the foot of the slope with the ice churnil 
heavily against the structure. However, by adjusting the sluice profile and 
selecting the right floor and ceiling elevation in the lower portion, it was 
found possible under all conditions to discharge the ice and water at high ve 
locity through the passage and into the tailwater without the hydraulic jump 
entering the passage. 

The general requirements for the drum-type gates controlling the ice 
sluice discharge were determined in Model No. 4 and further model tests 
were made by the gate manufacturer. To provide a discharge calibration fo 


Fig. 15. Powerhouse Ice Sluice Model 
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gates, Model No. 11 was constructed of the upper portions of the ice chute 
a natural scale of 1:36. Two 75-foot wide sluices and two 50-foot wide 
ices were reproduced and the discharge calibration obtained. 


Hydraulic Loading on Gates 


In a number of the model studies described earlier it was mentioned that 
ts were run to determine the vertical hydraulic loading on the water- 
trolling gates of the various structures. It is considered desirable to re- 
t on all these results in this section, as it is of interest to compare the 
ults obtained. In all seven different project gates were tested to determine 
vertical hydraulic loads which could be expected under the most severe 
‘rating conditions. Four of the gates were to be completely submerged in 
ration while the remainder were not. The submerged gates were the 
yerhouse headgate, the Long Sault Dam tunnel port gate, and the Massena 
ike service and emergency gates. Unsubmerged gates tested, were for the 
quois Dam and for the diversion sluices and spillways of the Long Sault 

n. Because of the interesting nature of the results, a summary of the test 
lings has been included in this paper. It should be noted that in all the test 
ults the effects of buoyancy, gate weight and friction in the checks have 

n eliminated, and the loads shown are therefore due to vertical hydraulic 
ion only. The apparatus used in the gate tests is shown in Fig. 9, where 
gate is suspended from a hoist through a spring balance used to measure 
vertical loading. Friction between the gate and the checks is eliminated 
sulling the gate upstream just free of the checks. This is accomplished by 
9ecial harness which assures a horizontal pull and also that the gate is 
ting just free of the checks. 


ts on Tunnel Port Gates at Long Sault Dam 


Perhaps the most interesting tests were those on the temporary gates, 

d to close the tunnel ports in the Long Sault Dam. While these gates were 
smallest of those tested, their weight being about 25 tons each, it was 

id that they were subject to the most severe loading conditions. As shown 
"ig. 16, investigations on the preliminary design indicated that a maximum 
1 of 154 tons could be expected if the holes required for eventual concreting 
1e ports were plugged, and 124 tons if the holes were left open. By alter- 
the shape of the bottom section to 45° and then to 60° as shown in Fig. 16, 
as possible to reduce the maximum downpull to 73 tons and 54 tons re- 
tively. 


ts on Other Project Gates 


‘he tests on the other project gates produced a variety of results which 

Shown in Fig. 17, where the grouping depends upon whether the gates were 

merged or not. While it is not practicable to analyse the results in detail, 

ain observations may be of interest. 

‘or unsubmerged gates conventionally installed with upstream skin plates, 
as at the Iroquois Dam and Long Sault Dam spillways, it would appear 

neither downpull or uplift is a significant force. The long Sault Dam di- 

sion sluices, being an adaption of the final sluice gates, were not con- 

ional in their installation because the sluice width was only 80 per cent 

: 
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that of the gate, with the result that water collected on the unsheeted down- 
stream members and thus increased the indicated downpull. 

The three submerged gates were generally similar in size and illustrate 
the effect of some of the varying factors. At Massena the skin plate was on 
the downstream side of the emergency gate and on the upstream side of the 
service gate. The large hydraulic downpull on the emergency gate, and the 
tendency for uplift to occur on the service gate, indicates that the position 0: 
the skin plate influences greatly the hydraulic loading. These loads may hax 
a considerable bearing on the required hoisting capacity or the ballast re- 
quired to ensure that the gate will close. The powerhouse headgate was test 
under turbine runaway conditions with the other two gates in the supply 
passages closed, and the large resulting downpull shows the effect of the in- 
crease in head across the gate. 
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CONCLUSIONS 


From the foregoing review, it is apparent that a very considerable effort 
has been made by the Power Authority of the State of New York and the Hydro: 
Electric Power Commission of Ontario to test thoroughly all aspects in the 
hydraulic design of the St. Lawrence Power Project. It is considered that the 
models have proven to be indispensable in providing designs which carry the 
maximum assurance of successful performance and at the same time make 
possible maximum economies. In the case of this Project, economies were 
effected which exceeded by many times the cost of the models. Information 
from the models enabled each step in construction to be taken with confidence. 
in the knowledge that its effect had been tested in advance and found to be 
within the criteria. Another considerable advantage of the models, was the 
ability to demonstrate in advance proposed plans and procedures to super- 
visory boards and other interested agencies. In all, it is felt that the models 
have played an essential part in the Project. 
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MECHANICAL ANALOGS AID GRAPHICAL FLOOD ROUTING? 
Closure by Max A. Kohler 


MAX A, KOHLER,! M. ASCE.—The author is indeed flattered by Mr. 
Clark’s discussion of his paper. The analogies drawn by Mr. Clark, as for 
sxample between the instantaneous hydrograph and the template routing, are 
realistic and accordingly lead to a better understanding of flood-wave phe- 
1omena. His caution that lag is in reality translatory storage is well taken. 


[alae 

. Proc. Paper 1585, April, 1958, by Max A. Kohler. 

. Chf. Research Hydrologist, U. s. Weather Bureau, Washington, D. C. 
¢t , ’ 
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NORTHEASTERN FLOODS OF 1955: RAINFALL AND RUNOFF2 


Closure by Tate Dalrymple 


TATE DALRYMPLE,! A. M. ASCE.—The author wishes to thank Mr. 
urner for his careful review of the paper. The discussion presents con- 
iderable material that adds to the subject of the paper, and thereby improves 


The author would hesitate to assign a frequency to a given flood that ne- 
2ssitated a long extrapolation of knowledge. Generally, a 50 year recurrence 
terval is about as far as is justified, except in a few instances where long 
rm knowledge is available, as in New England. A project now being com- 
eted by the Geological Survey has uncovered considerable historical flood 
formation in New England, and promises to form a basis for extending flood 
‘equency curves to 300 years. An extrapolation beyond this point seems un- 
ise, regardless of the need for evaluating greater floods. 


Proc. Paper 1662, June, 1958, by Tate Dalrymple. 
-Chf., Floods Section, Surface Water Branch, Water Resources Div., U. S. 
Geological Survey, Washington, D. C. 
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WAVE FORCES ON SUBMERGED STRUCTURES 


Discussions by I. Alterman and N. Hamlin 


I. ALTERMAN.!—This paper is a welcome addition to the vast but 

attered literature on the subject, as evidenced by the large number of Refer- 
ces following the article. 

The addition of semi-empirical data of horizontal forces on barge models 
a timely and very important contribution to problems arising in the ever- 
‘reasing number of sea-structures executed. The article is partly of inter- 
-for those who prefer a theoretical approach and partly for those who prefer 
yractical one. 

For compensation, it might be of interest that the writer has made use of 
theoretical formulae, such as (6), (7) of the article, or rather an inte- 
ition of those, and the results have shown a dramatical application: 

A series of 320 steel piles, 16" diameter, driven 33 feet into the sea-bed 
fe served as a temporary jetty. Once the superstructures had to be dis- 
ntled, the writer calculated the maximum waves which these piles could 
hstand freely, as winter storms approached. 

The computation showed that 16-ft. waves, with a period of 15 seconds 

ild cause breaking stresses to the piles. Orders were given accordingly 
cut-off the piles as quickly as feasible, before the stormy period would 
roach. 

Before all piles could be cut off, a first storm came, with waves 12 feet 

h and 15 seconds period striking the piles. An inspection of the piles 

wed no detrimental effect. The Resident-Engineer thus questioned the 
uracy of the formulae and slowed down the work. A check on the formulae 
this case proved that these waves could not cause over-stressing of the 
2S, because the bending-moments were proportional to the square of wave 
shts, 

The next storm actually brought 1 feet waves, while only 4 out of the 320 
*S were left. All four piles were bent at the sea-bed level beyond repair. 
The above full scale experiment may add evidence to the applicability of 
theory. 

Jn the other hand, it would be desirable to complete the article by intro- 
ing data relating to horizontal cylindrical submerged elements. 

axamples of such elements are: 


a) Horizontal connections of pile bend. 
b) Submarine pipes. 


Proc. Paper 1893, November, 1958, by Ernest F. Brater, John S. McNown 
and Leslie D. Stair. 

Head, Div. of Concrete Pipes, Water Planning for Israel Ltd. (TAHAL), 
Tel Aviv, Israel. 


184 May, 1959 HY 


(c) Floating submarine vessels. 
(d) Horizontal elements of “Jumbo” structures etc. 


A useful formula for the forces exerted by waves on an immersed cylinde 
may be found in literature. (1,2) 
Thus, the resistance to be taken up by the cylinder is: 


R= 477° g 9 rt o-4 -28f/e° ( 
where: 
g = acceleration due to gravity 


= sea-water density 


radius of cylinder 


p 
r 
c = wave celerity 
f 


= depth of centre of cylinder below mean water level. 


It is assumed that f = r, i.e. the cylinder is not too close to the free wate: 
level. 

In the special but common case that the cylinder is at the sea-bed(3) leve 
(e.g. a pipe), the above formula becomes very simple: 


2 st wih oth . 


Yar 
while 


we simply get: 


This simple formula shows the main difference between the action of 
waves on a plate or a vertical cylinder to a horizontal cylinder. In the 1 
case the wave action is shown to increase with the 4th power of the radius ; 
inversely with the square of the depth. This explains why a small diamet 
pipe may be laid without any support on the sea-bottom, while a larger one 
will suffer more than an equivalent vertical pile.(4) Also a submarine ve ; 
will “bounce” in shallow water while in deep water, since the action decrez 
with the square of the depth, it will be quite steady. _ 

Also, a flat plate will be affected by forces equivalent to the exposed ar 
(see e.g. formula (8) of the article), while a submerged horizontal cylinde 
may have forces proportional to a higher degree, e.g. the 4th power. _ 
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Y. HAMLIN.1—The authors’ tests are a welcome addition to the store of 
wledge available to the designer of offshore structures. To the writer’s 
wledge these tests were the first to measure wave forces methodically on 
arate fixed bodies representative of submerged barges in the absence of 
sorting piles. 

fhe significant conclusions from the tests and the authors’ analysis thereof 
first, the predominant wave forces on the barge-like models were inertial 
second, values of mass coefficient needed to determine these forces from 
retical calculations compared favorably with values estimated from purely 
retical considerations or determined experimentally by oscillating models 
till water. 

n order to estimate wave forces on a rectangular barge, the writer’s 

-tice has been to use a mass coefficient based upon data reported by Yu in 
authors’ Reference 19. An overall comparison of wave force data from 

el tests of a complete drill rig, together with preliminary data from the 
ors’ tests which were made available to the writer’s company as one of 
sponsoring organizations at the time of the tests, indicated that the Yu 
qula (modified to make it dimensionless) brought about reasonable corre- 
mn of the measured forces. However, mass coefficients from the authors’ 
iminary measurements of horizontal forces on the barge models were 
ificantly lower than those from the Yu formula. The comparison referred 
; reported in the writer’s discussion of a 1957 paper before the 

A.M.E.2 The final horizontal barge forces reported in the present paper 
considerably higher than the preliminary data indicated and lead to mass 
ficients which now agree favorably with the Yu formula. (The final values 
about 15% higher than formula values.) This tends further to confirm the 
ata as providing means for estimating mass coefficients for wave forces 
ubmerged rectangular barges with acceptable engineering accuracy. 
Owever, the results of tests with a flat plate are surprising in that the 

ial mass derived from the tests was high compared to theory. In addition, 
2 appeared to be a substantial drag force acting on the plate which could 
be explained by a drag coefficient higher than has been found from tests 
eady flow. Regarding the flat plate tests, a possible source of scatter in 
lata and of high transient forces might be the tendency of shed vortices 

€ vicinity of the model to retain their strength and thereby to continue to 
ence the forces developed on the plate. This would seem to be more of a 
effect rather than inertial. Have the authors any evidence in support of 
an explanation? In this connection, it is of interest that the flat plate 


ethlehem Steel Co., Shipbuilding Div., Central Technical Dept., Quincy, 
lass. 
echtin, Steele and Scales, P. 633-681, 1957 S.N.A.M.E. Transactions. 
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test data reported by Yu tends to confirm the theoretical value of C,,V fora 
infinitely long plate which, as the authors state, is equal to the volume of a 
circumscribing cylinder. Incidentally, with regard to the authors’ reser- 
vations concerning Cy, values obtained from high frequency oscillations, the 
Yu tests in water are reported for frequencies ranging from about .29 to pp | 
c.p.s., which are lower than the wave frequencies used in the authors’ tests. 
There are presently two interpretations of Eqs. (4) and (5) in use. Many 
who have worked with wave theory, including the authors, take the symbol 
(d +z) to mean the instantaneous particle height above the bottom. The writ 
is of the opinion that (d + z) should mean the average particle height, in othe 
words, the height of the center of the particle’s orbit. The latter interpreta 
enables one to depict orbital motion and other features of the wave in a mor 
realistic manner. An apparent confusion is indicated, but is of little conse- 
quence in analysing the authors’ tests, inasmuch as both interpretations leat 
to equal but opposite upstream and downstream horizontal inertial forces. ’ 
However, assuming (d + z) to represent the average particle height, theoreti 
cal vertical inertial forces are slightly larger when directed upward than wl 
downward for bodies fixed in the wave but near the surface. The tests with 
the barge models appear to confirm this tendency, but also show an unexpec 
trend for upward forces to exceed downward when the models were near the 
bottom. 
The authors are to be congratulated for a clear presentation of test resul 
which are of fundamental importance. . 7 
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SNOWMELT RUNOFF? 


Discussion by H. C. Riggs 


iac. RIGGS,! A. M. ASCE.—The authors propose a method of computing 
rate of runoff from melting snow on the premise that, for practical use, 
off from melt must be related to generally available meteorological data. 
yputing snowmelt runoff on a daily basis requires first a relation between 
actual melt water produced on a unit area and the causative factors. The 
relation probably would be one derived from data obtained on a small 

on which melt water and the meteorological factors were accurately 
sured. The second part of the problem is that of routing the melt water 
stream channel and thence downstream to the gage site. 

‘he subject paper treats these two parts as one. A disadvantage of this 
-oach is that the reliabilities of the several procedures (transferrence of 
eorological data, distribution of runoff from one day’s melt to streamflow 
several days, altering of discharge data on the beginning day of melt to 

ct a trial-and-error solution) cannot be evaluated separately. 

‘he authors’ results are incomplete for the purpose of verifying their pro- 
ires. Stream discharge is not related to melt potential alone; the amount 
vailable water in the snow pack and the precipitation during the melting 
od must be included in order to define a usable relation. The wide di- 
‘ence of the curves of the authors’ Fig. 8 indicates the need for additional 
ables. The equation 


Y = 27.4X - 1.88 x2 


ely puts the divergent curves of Fig. 8 into register at a melt potential of 
nches. This equation was an acceptable fit in only 7 of the 15 years tried 
therefore, does not constitute verification of the method. No equation of 
type, however reliable, could be used for forecasting because results are 
1, not in quantitative forms, but as a percentage of.cumulative runoff at 
nches melt potential. This cumulative runoff at 7.0 inches melt potential 
known at the time of forecast. 

he most reliable data available should be used in developing and testing a 
method. It is unfortunate that the streamflow records in the selected 

1 are not of highest accuracy. However, the use of stream discharges re- 
i by the authors does not lend confidence to their results. The “ice index” 
od is not new. It is a simplification of a graphic method described by 

(1) in 1913. Present practice is described by Moore. 2) 


TOC. Paper 1834, November, 1958, by J. Harold Zoller and Arno T. Lenz. 
ydr. Engr., U. S. Geological Survey, Washington 25, D. C. 
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A method of forecasting the rate of runoff from snowmelt is of practica 
use only when the causative factors can be predicted in advance or when thi 
is appreciable lag between occurrence of the causative factors and the resu 


ing runoff. Compliance with one of these necessary conditions needs to be 
demonstrated. 


‘ 
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ia), 1825 (SM4), 1826(SM4), 1827(ST6)°, 1828(SM4)°, 1829(HW3)°, 1830(PO5)°, 1831(EM4)°, 1832(HY5) 


MBER: 1833(HY6), 1834(HY6), 1835(SA6), 1836(ST7), 1837(ST7), 1838(ST7), 1839/ST7), 1840(ST7), 
41(ST7), 1842(SU3), 1843(SU3), 1844(SU3), 1845(SU3), 1846(SU3), 1847(SA6), 1848(SA6), ares 
st 1851‘SA6), 1852(SA6), 1853(SA6), 1854(ST7), 1855(SA6)°, 1856(HY6)°, 1857(ST7)°, 1858 


BER: 1859(HY7), 1860(IR4), 1861 (IR4), 1862(1R4), 1863(SM5), 1864(SM5), 1865(ST8), 1866(ST8), 1867 
}), 1868(PP1), 1869(PP1), 1870(PP1), 1871(PP1), 1872(PP1), 1873(WWws), 1874(WW5),1875(WW5), 1876 
5), 1877(CP2), 1878(STS), 1879(ST8), 1880(HY7), 1881(SM5)°, 1882(ST8)°, 1883(PP1)°, 1884(Wws5)°, 
CP2)¢, 1886(PO6), 1887(PO6), 1888(PO6), 1889(PO6), 1890(HY7), 1891(PP1), 


VOLUME 85 (1959) 


RY: 1892(AT1), 1893(AT1), 1894(EM1), 1895(EM1), 1896(EM1), 1897(EM1), 1898(EM1), 1899(HW1), 
), 1901(HY1), 1902(HY1), 1903(HY1), 1904(HY1), 1905{PL1), 1906(PL1), 1907(PL1), 1908(PL1), 
1), 1910(ST1), 1911(ST1), 1912(ST1), 1913(ST1), 1914(ST1), 1915(ST1), 1916(AT1)°, 1917(EM1)°, 
iA 1919(HY1)°, 1920(PL1)°, 1921(SA1)©, 1922(ST1)©, 1923(EM1), 1924(HW1), 1925(HW1), 1926 
1927(HW1), 1928(HW1), 1929(SA1), 1930(SA1), 1931(SA1), 1932(SA1). ¥ 
RY: 1933(HY2), 1934(HY2), 1935(HY2), 1936(SM1), 1937(SM1), 1938(ST2), 1939(ST2), 1940(ST2), 
T 2), oanerae 1943(ST2), 1944(ST2), 1945(HY2), 1946(PO1), 1947(POL), 1948(PO1), 1949(PO1), 
2)°, 1951(SM1)©, 1952(ST2)°, 1953(PO1)¢, 1954(CO1), 1955(CO1), 1956(CO1), 1957(CO1), 1958 
1), 1959(CO1). 
ie a  1961(HY3), 1962(HY3), 1963(IR1), 1964(IR1), 1965(IR1), 1966(IR1), 1967(SA2), 1968 
fee Sanaa, ition 1972(ST3), 1973(ST3), 1974(ST3), 1975(ST3), 1976(WWw1), 1977 
; sero), 1979(WW1), 1980(Wwi), 1981(WW1), 1982(WW1), 1983(WW1), 1984(SA2), 1985(SA2)°, 


4), 1995(HY4), 1996(HY4), 1997(HY¥4). 1998 
arin 2004(ST4) case 2006(PO2), 2007 
3(P02)°. Hee 
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Me 
Paw. 4 
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